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ABSTRACT 


The  CMARC  panel-code  is  evaluated  for  the  development  of  an  aerodynamic 
model  of  the  Naval  Postgraduate  School  FROG  Unmanned  Air  Vehicle  (UAV). 
CMARC  is  a  personal  computer  hosted  panel-code  software  suite  for  solving  inviscid, 
incompressible  flow  over  complex  three-dimensional  bodies.  A  panel  model  of  the  NPS 
FROG  UAV  is  developed  to  obtain  stability  derivative  data  at  the  cruise  flight  condition. 
Emphasis  is  placed  on  comparing  the  CMARC  data  to  aerodynamic  models  obtained 
from  classical  design  techniques  and  parameter  estimation.  Linearized  longitudinal  and 
lateral-directional  state-equation  models  are  used  to  compare  the  dynamic  response  of 
each  data  set.  In  addition,  CMARC  is  used  to  generate  static-source  and  angle-of-attack 
sensor  position  corrections.  Position  corrections  are  provided  in  look-up  table  and  curve- 
fit  formats.  The  aerodynamic  model  obtained  with  CMARC  demonstrated  higher  fidelity 
dynamic  longitudinal  response  than  the  classical  design  model.  Dynamic  lateral- 
directional  response  is  similar  to  that  obtained  from  classical  design  techniques. 
Adjustment  through  comparison  with  flight-test  data  is  still  required  to  optimize  the 
CMARC  model.  Future  studies  should  concentrate  on  improving  CMARC  modeling  of 
fuselage  side  force  through  the  addition  of  wake  separation  lines.  Additionally,  the 
propeller  disk  should  be  modeled  in  an  attempt  to  capture  the  effects  of  increased 
dynamic  pressure  over  the  horizontal  and  vertical  tail  surfaces. 
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I.  INTRODUCTION 

A.         BACKGROUND 

Computational  fluid  dynamics  (CFD)  is  increasingly  used  as  a  design  and  analysis 
tool.  As  the  price  of  computer  hardware  drops  and  computational  power  increases,  CFD 
becomes  more  attractive  to  a  larger  audience.  CFD  tools  range  from  the  high-end  three- 
dimensional  (3D)  Navier-Stokes  solvers  for  compressible,  viscous  fluids  to  potential  flow 
solvers  for  incompressible,  inviscid  flows.  This  paper  discusses  the  development  of  an 
aerodynamic  model  of  the  Naval  Postgraduate  School  (NPS)  Fiber  Optic  Guided  (FROG) 
Unmanned  Air  Vehicle  (UAV)  using  a  personal  computer  hosted  panel  code. 

Validation  of  the  Personal  Simulation  Works  panel  code  software  and  initial 
FROG  UAV  modeling  was  completed  in  June  1997.  Results  are  reported  in  Ref.  [1]. 
This  study  expands  on  the  modeling  effort  to  include  the  development  of  a  complete 
aerodynamic  model  at  the  trim  flight  condition.  Stability  derivative  data  are  obtained 
from  the  panel  code  and  then  compared  to  data  from  classical  design  calculations  and 
parameter  estimation.  A  linearized  state  space  simulation  is  used  to  extract  the  air  vehicle 
dynamic  modes  and  to  model  control  response.  Results  are  compared  for  all  three  data 
sets. 

The  Personal  Simulation  Works  software  suite,  consisting  of  LOFTSMAN, 
CMARC  and  POSTMARC,  is  used  for  all  aspects  of  the  study.  The  software  provides 
for  panel  model  development,  input  file  processing  and  the  visualization  of  results. 
Emphasis  is  placed  on  verifying  both  the  accuracy  and  suitability  of  the  CFD  programs 
for  aerodynamic  modeling. 

Until  recently,  personal  computers  (PC)  did  not  have  the  computational  power  or 
memory  to  be  practical  for  panel  code  CFD  programs.  Things  have  changed  with  the 
introduction  of  the  Pentium  class  PC  and  low  cost  RAM.  AeroLogic  capitalized  on  the 
power  of  the  Pentium  class  PC  and  developed  Personal  Simulation  Works  (PSW).  PS W 
is  based  on  the  3D,  low  order,  inviscid  potential  flow  solver  named  CMARC.  CMARC  is 
a  re-hosted  version  of  NASA's  Panel  Method  Ames  Research  Center  (PMARC). 
PMARC  was  re-written  in  the  C  language  and  compiled  for  IBM  compatible  PCs. 
CMARC  runs  under  the  DOS  operating  system.  It  will  also  run  in  a  DOS  window  under 
the  WINDOWS  3.x,  95  or  NT  operating  systems.  Enhanced  capabilities  include; 
improved  memory  management,   an  expanded  set  of  command   line   switches   and 


provisions  for  expanded  boundary  layer  post-processing  capabilities.  However,  the  core 
processing  algorithms  remain  the  same  as  implemented  in  PMARC. 

LOFTSMAN,  the  PSW  pre-processing  program,  is  used  to  mesh  complex  3D 
bodies  and  create  input  file  patches.  The  program  runs  under  the  Windows  operating 
system  and  allows  the  user  to  loft  conies  based  3D  surfaces.  The  program  automatically 
creates  CMARC,  PMARC  or  VSAERO  input  patches  based  on  desired  panel  densities 
and  distribution. 

POSTMARC  is  used  for  flow  visualization  and  integration  of  resultant  forces.  It 
runs  under  the  Windows  3.x,  95  and  NT  operating  systems.  POSTMARC  reads  CMARC 
or  PMARC  output  files  and  provides  for  the  visualization  of  model  geometry,  wake 
stepping,  on  and  off-body  streamlines  and  surface  phenomena. 

B.  REQUIREMENTS 

The  Naval  Postgraduate  School  Aeronautics  Department  is  integrating  UAV 
hardware  and  software  to  demonstrate  autonomous  flight,  trajectory  tracking  and 
automatic  landing.  Closed-loop  flight  control  development  requires  a  valid  aerodynamic 
truth  model  for  the  UAV  airframe.  The  introduction  of  each  new  airframe  requires  the 
development  of  a  new  aerodynamic  truth  model.  Most  recently,  Papageorgiou  [Ref.  2] 
developed  and  tested  an  aerodynamic  model  for  the  NPS  FROG  UAV  based  on  classical 
methods.  His  model  produced  a  close  match  to  flight  test  results  in  the  longitudinal  axis. 
However,  the  lateral-directional  axis  required  modifications  based  on  measured  flight  test 
data  to  produce  acceptable  results.  With  the  availability  of  low  cost  panel  code  CFD 
capabilities,  it  is  suggested  that  a  panel  code  model  of  the  FROG  UAV  will  provide  more 
refined  stability  derivative  data  for  an  initial  aerodynamic  model. 

In  addition  to  a  valid  truth  model,  accurate  pitot-static  and  angle-of-attack  sensors 
are  required  for  highly  augmented  flight  control  systems.  CMARC  is  well  suited  for 
solving  on-body  static  pressure  distributions  and  off-body  flow  velocities  over  the 
predominately  attached  flow  fields  of  a  fuselage  forebody.  This  proves  particularly 
useful  for  generating  pitot-static  and  angle  of  attack  correction  curves  and  look-up  tables. 

C.  STATEMENT  OF  OBJECTIVES 

The  Naval  Postgraduate  School  Aeronautics  Department  has  both  active  CFD 
research  and  avionics  development  programs.   The  primary  purpose  of  this  investigation 


is  to  verify  the  accuracy  and  suitability  of  the  PSW  software  suite  for  developing  an 
aerodynamic  model  of  the  NPS  FROG  UAV.  Specific  objectives  are  as  follows: 

•  Demonstrate  panel  code  modeling,  processing  and  visualization  on  a  Pentium 
PC  using  the  PSW  software  package. 

•  Develop  and  analyze  a  panel  code  model  for  the  NPS  FROG  UAV  using  PSW 
to  obtain  a  complete  set  of  stability  derivatives  at  the  cruise  trim  condition. 

•  Verify  the  aerodynamic  model  through  comparison  to  data  obtained  from 
classical  design  calculations  and  parameter  estimation. 

•  Demonstrate  techniques  for  producing  angle-of-attack  vane  and  pitot-static 
correction  curves. 

•  Develop  a  user  guide  that  outlines  PSW  panel  code  modeling  for  the 
extraction  of  stability  derivatives. 


II.  OVERVIEW  OF  PERSONAL  SIMULATION  WORKS 

A         GENERAL 

Personal  Simulation  Works  is  a  PC  based  software  suite  that  provides  for  the  three 
primary  CFD  requirements;  3D  modeling  of  an  aircraft  (LOFTSMAN),  panel  code  flow 
solver  (CMARC),  and  post-processing  of  the  computed  flow  field  (POSTMARC).  The 
software  package  contains  three  applications  hosted  on  the  IBM  compatible  personal 
computer.  Each  software  program  is  discussed  separately. 

B.         LOFTSMAN 

LOFTSMAN  is  a  Windows  hosted  3D  modeling  tool  that  generates  surface  panel 
distributions  for  CMARC  or  PMARC  input  files.  The  program  is  based  on  conies,  which 
allows  rapid  lofting  of  streamlined  bodies  such  as  aircraft  fuselages  and  engine  nacelles. 
In  addition,  wing  and  control  surfaces  can  be  designed  with  the  extensive  library  of  airfoil 
templates  or  with  user  specified  coordinates.  The  software  is  well  documented,  including 
a  tutorial,  in  the  Personal  Simulation  Works  User  Guide  [Ref.  3].  LOFTSMAN  is 
primarily  designed  for  creating  new  objects,  but  an  existing  airframe  can  be  matched 
quite  closely  with  just  a  detailed  three-view  drawing  that  includes  frame  cross  sections. 
Appendix  A  outlines  the  development  of  the  FROG  UAV  model. 

1.  Streamlined  Bodies 

LOFTSMAN  functionality  is  divided  into  Body  Objects  and  Wing  Objects.  In 
general,  they  remain  separate  unless  the  intersection  between  a  wing  and  body  is  required. 

Body  Objects  are  created  using  a  family  of  curves  called  second-degree  conies. 
Circles,  ellipses,  parabolas  and  hyperbolas  are  among  this  group.  An  entire  fuselage  is 
described  by  specifying  just  four  lines.  These  are  the  top  waterline  (TW),  bottom 
waterline  (BW),  the  maximum  breadth  line  (MB)  and  the  waterline  of  the  maximum 
breadth  line  (WW).  For  each  line,  the  beginning,  ending  and  a  few  points  along  the  line 
are  specified.  Control  points  are  also  specified  with  a  curvature  factor  that  allows 
LOFTSMAN  to  generate  a  smooth  conic  between  the  points.  The  power  of  conic  lofting 
will  become  evident  when  discussing  the  modeling  of  the  complex  FROG  UAV  fuselage 
in  Chapter  V. 


2.  Wings  and  Control  Surfaces 

Wings  and  control  surfaces  are  easily  specified  in  LOFTSMAN  using  a  short 
input  file  created  with  any  text  editor.  The  file  specifies  root,  intermediate  and  tip  rib 
section,  location,  axis,  chord  and  incidence.  LOFTSMAN  then  fairs  a  smooth  surface 
through  the  rib  sections.  Washout  is  specified  by  varying  the  incidence  of  the  root  and  tip 
ribs.  Sweep-back  is  controlled  by  staggering  the  tip  rib  location  with  respect  to  the  root 
rib.  Once  the  general  wing  surface  is  specified,  control  surfaces  such  as  ailerons,  flaps 
and  elevators  can  be  deflected  and  meshed. 

3.  Patches 

LOFTSMAN  automatically  meshes  3D  surfaces  and  creates  patches  for 
CMARC/PMARC  input  files.  The  distinction  between  a  mesh  and  a  patch  is  important. 
A  mesh  is  a  set  of  quadrilateral  and  triangular  panels  that  represent  the  surface  of  a  wing 
or  body.  When  the  set  of  panels  is  organized  and  formatted  to  create  a  sub-component 
portion  of  a  CMARC  or  PMARC  input  file,  it  is  called  a  patch. 

A  body  or  wing  surface  is  first  meshed  at  a  density  specified  by  the  user.  Panel 
compression  options  include  cosine  and  half-cosine  spacing.  After  meshing  the  object, 
one  saves  it  to  a  text  file  as  a  formatted  patch.  One  then  opens  the  patch  file  with  any  text 
editor  and  copies/pastes  the  patch  text  into  the  appropriate  location  in  the  CMARC  input 
file. 

Each  control  surface  deflection  requires  a  separate  mesh  and  formatted  patch.  For 
instance,  to  evaluate  roll  performance  one  needs  to  separately  mesh  an  upward  aileron 
deflection  on  the  right  wing  and  a  downward  deflection  on  the  left  wing.  If  multiple 
deflections  of  a  single  control  surface  are  required,  each  deflection  must  be  meshed 
separately. 

C.         CMARC 

CMARC  is  the  C  version  of  PMARC  low-order,  3D  panel  code.  Inviscid, 
irrotational,  incompressible,  potential  flow  is  assumed.  Low-order  means  that  source  and 
doublet  strength  distribution  is  constant  across  a  panel.  There  is  no  attempt  to  match  the 
source  or  doublet  strength  of  an  adjacent  panel  at  a  common  edge.  Advanced  features 
include  internal  flow  modeling  and  time  stepping  wake  models. 


PMARC  version  12.19  was  released  as  FORTRAN  77  source  code  in  1992. 
CMARC  was  rewritten  in  the  C  language  and  compiled  for  hosting  on  IBM  compatible 
personal  computers  by  AeroLogic,  Inc.  The  program  runs  under  the  DOS  operating 
system.  It  will  also  run  in  a  DOS  window  from  Windows  3.1,  95  or  NT.  Enhanced 
features  include  command-line  options  and  flexible  memory  management.  Command 
line  options  simplify  batch  processing  by  adding  an  extensive  set  of  switches  that  can  be 
set  external  to  the  CMARC  input  file.  Flexible  memory  management  provides  for  the 
automatic  sizing  of  arrays  without  having  to  recompile  the  source  code. 

D.         POSTMARC 

POSTMARC  is  a  Windows  post-processing  program  for  the  visualization  of 
CMARC  and  PMARC  output  files.  Capabilities  include  body  geometry,  wake  stepping, 
surface  pressure  and  streamline  visualization.  POSTMARC  also  provides  the  capability 
to  integrate  pressure  and  skin  friction  forces  over  the  model  geometry.  This  proves 
particularly  useful  when  one  desires  to  recalculate  loads  around  a  different  center  of 
gravity. 

An  interesting  feature  for  design  work  is  the  integration  of  panel  surface  area  to 
obtain  total  wetted  area.  After  lofting  a  new  geometry  in  LOFTSMAN,  a  quick  check  of 
geometry  is  made  by  running  CMARC  with  the  -g  command  line  toggle.  The  total 
wetted  area  is  then  checked  in  POSTMARC.  This  function  is  particularly  useful  when 
working  to  reduce  skin  friction  drag. 

Versions  1.17.3  and  later  of  POSTMARC  include  the  capability  to  integrate  skin 
friction  drag  coefficient  over  the  model  geometry.  It  is  important  to  note  that  a  key  piece 
of  the  drag  equation  is  missing  from  a  POSTMARC  solution.  CMARC  provides  induced 
drag  from  the  surface  pressure  distribution  and  skin  friction  drag  from  the  2D  boundary 
layer  code.  Skin  friction  is  only  calculated  up  to  the  point  of  boundary  layer  separation. 
Pressure  drag  due  to  separation,  a  major  portion  of  the  drag  equation,  is  missing  from  a 
CMARC/POSTMARC  solution. 

In  fact,  if  one  isn't  careful,  POSTMARC  drag  calculations  can  be  misleading. 
Take  for  instance  two  similar  model  configurations  with  only  minor  geometry  differences 
that  do  not  affect  wetted  area.  It  is  possible  for  the  model  with  more  flow  separation  to 
have  less  skin  friction  drag  because  there  is  no  CMARC  output  for  skin  friction 
coefficient  after  the  boundary  layer  code  predicts  separation.    During  iterative  design 


work,  this  could  lead  to  the  incorrect  conclusion  that  the  design  team  is  reducing  overall 
drag.  Perhaps  a  better  function  for  LOFTSMAN  than  integrated  skin  friction  drag  would 
be  a  function  that  predicts  the  percentage  of  attached  flow  and  laminar  flow.  Iterative 
design  changes  could  be  made  that  maximize  laminar  flow  and  minimize  separated  flow. 


III.  CMARC  PANEL  CODE  THEORY 

A.         POTENTIAL  FLOW  PANEL  CODE  THEORY  (CMARC/PMARC) 

Potential  flow  theory  involves  the  superposition  of  sources  and  doublets  to 
generate  the  desired  flow  field  around  a  3D  body.  It  assumes  inviscid,  irrotational  and 
incompressible  flow.  As  such,  valid  solutions  are  only  obtained  at  low  Mach  numbers 
and  for  flow  fields  without  large  areas  of  separation. 

The  basic  concept  of  panel  methods,  as  outlined  by  Bertin  and  Smith  [Ref.  4], 
requires  the  modeling  of  the  desired  3D  configuration  with  a  large  number  of 
quadrilateral  and  triangular  panels  representing  the  surface  of  the  aircraft.  A  series  of 
sources,  doublets  and  vortices  is  then  distributed  on  each  panel.  Superposition  allows  the 
simultaneous  computation  of  the  singularity  strengths  required  to  satisfy  flow  tangency 
on  the  surface.  The  inviscid,  irrotational  and  incompressible  flow  field  represented  by  the 
superposition  of  sources  and  doublets  satisfies  the  Laplace  equation: 

V2(D  =  0  3.1 

Using  Green's  Theorem,  the  potential  at  any  point  P  in  the  flow  is  represented  by: 

•«•  =  i  JJfa-O'H^-i  Jl(;)«  (vo- v<D,ys  3.2 


s+w  s+w 


Where     (o-O,)     represents    the    potential     from    the    doublet    distribution    and 
n  ■  (VO  -  VOj  represents  the  potential  from  the  source  distributions. 

CMARC  is  a  low  order  panel  code  that  assumes  constant  source  and  doublet 
strength  distributions  across  each  panel.  Figure  3.1  shows  a  panel  layout  for  a  generic  3D 
wing  fuselage  configuration.  It  is  important  to  note  that  for  a  3D  solution,  there  is  an 
equivalence  to  surface  doublet  and  surface  vortex  distributions.  CMARC  implements 
source  and  doublet  distributions. 


Figure  3.1  Typical  Wing/Body  Panel  Code  Configuration,  from  Ref.  [5]. 

As  mentioned  previously,  the  general  boundary  condition  imposed  is  tangential 
flow  at  the  surface.  CMARC,  as  outlined  in  Ref.  [3],  allows  the  modification  of  this 
boundary  condition  on  individual  panels  or  groups  of  panels.  A  normal  surface  velocity 
distribution  may  be  specified  to  simulate  flow  into  or  out  of  ducts. 

In  order  to  produce  lift,  a  potential  flow  panel  code  requires  a  method  to 
implement  the  Kutta  condition.  As  noted  in  Anderson  [Ref.  6],  the  Kutta  condition  at  the 
trailing  edge  implies  that  the  circulation,  T,  around  an  airfoil  is  such  that  the  flow  exits 
the  trailing  edge  smoothly.  In  addition,  the  velocities  leaving  the  top  and  bottom  surfaces 
are  finite  and  equal  in  magnitude  and  direction. 

Panel  codes  impose  the  Kutta  condition  by  the  shedding  of  wake  panels  along  the 
trailing  edges  or  separation  lines.  Wake  panels  are  similar  to  a  surface  panel  with  only  a 
doublet  distribution.  The  doublet  strength  of  the  attached  wake  panel  equals  the 
difference  in  doublet  strengths  of  the  two  adjacent  surface  panels. 


K) 


The  CMARC  core  panel  code  processing  engine  is  functionally  equivalent  to  the 
PMARC  panel  code  module.  The  implemented  equations  are  well  documented  by  Ashby 
et  al.  [Ref.  5].  The  PMARC  documentation  includes  a  wing-body  combination,  shown  in 
Figure  3.1,  evaluated  by  PMARC  with  good  correlation  to  experimental  data.  The 
results  are  shown  in  Figures  3.2  and  3.3.  In  addition,  Lambert  [Ref.  7]  compared 
PMARC  panel  code  results  to  several  theoretical  and  experimental  test  cases  with  good 
correlation  at  low  angle-of-attack.  Sensitivity  to  wake  placement  is  highlighted  by  his 
studies. 

Wake  positioning  can  have  a  large  influence  on  potential  flow  solutions.  A  wake 
is  obviously  attached  to  the  trailing  edge  of  wings  and  control  surfaces  with  sharp,  thin 
trailing  edges  to  produce  the  Kutta  condition.  However,  wake  positioning  on  streamlined 
fuselages,  missile  airframes  and  nacelles  is  more  of  an  art  than  science.  Recently,  Tuncer 
and  Platzer  [Ref.  8]  investigated  generalized  wake  placement  techniques  for  cylindrical 
bodies  of  revolution  with  good  correlation  to  experimental  data  at  up  to  20  degrees  angle- 
of-attack.  The  techniques  were  used  with  success  by  this  author  in  Ref.  [1]  for  the 
verification  of  CMARC  calculations  for  flow  over  an  inclined  6:1  prolate  spheroid. 
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Figure  3.2     Comparison  of  Experimental  Data  and  PMARC  Results  for  Two 
Spanwise  Stations  of  the  Wing/Body  (a  =  4°),  from  Ref.  [5]. 
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Figure  3.3     Comparison  of  Experimental  Data  and  PMARC  Along  the  Fuselage 
Centerline  of  the  Wing/Body  Configuration  (a  =  4°),  from  Ref.  [5]. 


B. 


CMARC  BOUNDARY  LAYER  ANALYSIS  THEORY 


CMARC  and  PMARC  use  the  same  two-dimensional  integral  method  to  calculate 
boundary  layer  characteristics  along  a  surface  streamline.  A  transition  model 
automatically  switches  from  laminar  to  turbulent  calculations.  The  developers  of  the 
PMARC  code  chose  a  2D  integral  routine  over  a  3D  finite  difference  grid  method  due  to 
speed  and  robustness  of  the  calculations  [Ref.  5].  Building  a  finite  difference  grid  is  a 
difficult  and  time  consuming  process  requiring  the  user  to  develop  grids  over  complex  3D 
surfaces.  In  addition,  boundary  layer  calculation  times  can  easily  exceed  that  required  for 
the  basic  potential  flow  solution.  Reference  [9]  gives  a  good  outline  of  three-dimensional 
finite  difference  methods. 


12 


The  CMARC  and  PMARC  User's  Guides  [Refs.  3  and  5]  contain  detailed 
discussions  on  the  development  of  the  CMARC/PMARC  boundary  layer  code  starting 
from  the  two-dimensional  momentum  equation: 

dO     ,  ,    6  dU     1 

—  +2  +  //)* =  -C, 

dri  U  dr/      2    f 

3.3 
The  momentum  integral  equation  is  numerically  integrated  along  a  surface  streamline. 

The  laminar  region  of  the  boundary  layer  is  modeled  by  numerically  integrating 
the  following  exact  differential  equation.  The  equation  is  solved  iteratively  through 
numerical  integration  along  a  streamline  starting  at  a  stagnation  point  [Ref.  5]: 


0.45v 


<K*Y  =  ^  J(l  +  2.222g(K,  ^(rjf  dr?  +  0(0) 


3.4 


1/(0) 

where:  U  -  velocity  at  outer  edge  of  boundary  layer 

9  -  momentum  thickness 

v  drj 
r)  -  generalized  coordinate  along  a  streamline 


The  value  g(K,u)  is  based  on  exact  solutions  for  a  number  of  pressure 
distributions.  Initial  work  was  conducted  by  Thwaites  with  improvements  by  Curie 
[Ref.  10]: 

g(K,ju)  =  F0(K)-MG0(K)-0A5  +  6K  3.5 

CMARC  uses  an  empirical  transition  model  based  on  the  average  pressure 
gradient,  K  ,  for  predicting  laminar  to  turbulent  transition.  The  following  relations  are 
used  to  calculate  the  transition  point  [Ref.  5]: 
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\Kdrj 
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K  =  ^ 3.6 

where    r/ms   is  the  streamline  coordinate  at  instability.     And,  K  is  the  local  pressure 
gradient  at  boundary  layer  instability  [Ref.  5]: 


K  =  -0.4709  +  0.1 1066*ln(Re,)  +  0.0058591  *ln2(Re5)  (0<Re,  <650) 

3.7 

K  =  0.69412-  0.23992 *ln(Re5)  +  0.0205* ln2(Re e)  (650  <  Re^  <10000) 


The  local  Reynolds  number  at  transition  is  correlated  to  K  with  the  following  expressions 
[Ref.  5]: 

K  =  -0.0925  +  0.00007  *  Re^  (0  <  Re^  <  750) 

K  =  -0.12571  +  0.0001 14286 *Re„  (750  <  Ree  <1100)      3.8 

Z  =  1.59381-0.45543*ln(Re<?)  +  0.032534*ln2(Re(?)  (llOO  <  Ree  <3000) 

At  transition,  the  initial  turbulent  shape  factor,  H,    is  given  by  the  following  empirical 
formula  that  is  a  fit  to  data  developed  by  Coles  [Ref.  10]: 

1.4754 
H  = +  0.9698  3.9 

logio(Re*) 

Provisions  are  made  to  check  for  turbulent  reattachment  if  laminar  separation  is 
encountered.  At  laminar  separation,  a  point  calculation  is  made  to  determine  if  the 
boundary  layer  will  reattach.  If  reattachment  is  predicted,  the  boundary  layer  code 
immediately  switches  to  turbulent  calculations.  No  attempt  is  made  to  model  the  laminar 
separation  bubble  or  provide  a  transition  length.  After  laminar  separation  is  predicted, 
the  following  empirical  relations  are  used  to  determine  if  reattachment  occurs  [Ref.  5]: 
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K  =  0.0227  -  0.007575  *  Ree  -  0.000001 157  *  Re/  (Re^  >  125) 

3.10 

A"  =  -0.09  (Re^<125) 


The  boundary  layer  code  in  CMARC  uses  a  point  transition  model.  No  attempt  is 
made  to  model  a  more  representative  transition  length.  Turbulent  calculations  begin  at 
transition  using  the  Nash-Hicks  model  [Ref.  5].  Calculations  continue  along  the 
streamline  until  turbulent  separation  is  predicted  or  the  end  of  the  streamline  is  reached. 
No  boundary  layer  data  is  available  after  separation. 

The  authors  of  PMARC  caveat  that  their  boundary  layer  calculations  are  quite 
accurate  for  predominately  2D  flow  but  break  down  in  regions  of  large  cross  flow  near 
separation.  This  premise  was  tested  in  Reference  [1]  by  comparing  the  predominately  2D 
flow  over  the  inboard  region  of  a  high  aspect  ratio  wing  to  the  finite  difference 
calculations  performed  by  the  Naval  Postgraduate  School  Unsteady  Potential  Flow  Code 
(UPOT).  In  general,  CMARC  provided  correct  trends  for  both  the  transition  and 
separation  points.  However,  in  all  cases,  CMARC  predicted  early  transition  and  late  flow 
separation.  As  expected,  the  differences  were  greatest  at  lower  Reynolds  numbers  where 
boundary  layer  thickness  is  larger. 

Reference  [1]  also  modeled  flow  over  an  inclined  prolate  spheroid.  The  6:1 
prolate  spheroid  was  chosen  because  extensive  experimental  data  is  available.  In 
addition,  the  three-dimensional  flow  around  a  prolate  shperoid  is  similar  to  a  streamlined 
slender  fuselage.  With  proper  wake  placement,  CMARC  was  found  to  produce  accurate 
normal  force  and  pitching  moment  coefficients.  Over  the  three-dimensional  body, 
CMARC  boundary  layer  calculations  also  predicted  early  transition  and  late  flow 
separation.  Despite  inaccuracies,  CMARC  boundary  layer  calculations  remained  useful 
when  used  as  a  design  tool  for  visualizing  the  trend  in  transition  and  separation  points 
with  configuration  changes. 
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IV.  AERODYNAMIC  MODEL  OF  THE  FROG  UAV 

A.  BACKGROUND 

The  Naval  Postgraduate  School  Aeronautics  Department  is  integrating  UAV 
hardware  and  software  to  demonstrate  autonomous  flight,  trajectory  tracking  and 
automatic  landing.  A  core  requirement  for  flight  control  law  development  is  a  valid 
aerodynamic  truth  model  for  the  UAV  airframe.  A  panel  code  model  of  the  FROG  UAV 
is  one  method  for  estimating  many  of  the  stability  derivatives  required  for  an 
aerodynamic  truth  model.  This  development  effort  concentrates  on  finding  the  static,  rate 
damping  and  control  power  derivatives  for  both  the  longitudinal  and  lateral-directional 
axes. 

Panel  code  modeling  usefulness  goes  beyond  the  development  of  aerodynamic 
coefficients.  Flight  control  systems  require  accurate  pitot-static  and  angle-of-attack 
sensor  inputs.  CMARC  accurately  solves  on-body  static  pressure  distributions  and  off- 
body  flow  velocities  over  the  predominately  attached  flow  fields  of  fuselage  fore  bodies. 
In  this  study,  correction  curves  are  generated  for  static  source  and  angle-of-attack  probe 
position  errors. 

B.  FROG  UAV  DESCRIPTION 

The  FROG  UAV  is  a  small  single  engine  flight  test  vehicle  used  for  autonomous 
flight  research  by  the  Naval  Postgraduate  School  Aeronautics  Department.  The  aircraft 
was  originally  designated  the  FOG-R  by  the  U.  S.  Army.  It  was  designed  as  a  small 
lightweight,  battlefield  observation  platform  that  could  be  guided  by  a  fiber  optic  data 
link.  Table  4.1  presents  the  basic  aircraft  specifications. 

The  aircraft  is  somewhat  unconventional  in  that  the  engine  is  mounted  in  a  nacelle 
tractor  style  above  the  fuselage  and  wing.  The  aft  fuselage  consists  of  a  1.75  inch 
diameter  aluminum  tube  which  connects  the  tail  surfaces  to  the  main  fuselage.  Figure  4.1 
displays  a  three  view  drawing  of  the  FROG  UAV. 
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PARAMETER 

MEASUREMENT/UNITS 

Length 

8.125  ft                 97.5  in 

Height 

1.75  ft                    21  in 

Weight 

67.7  lbs 

Power  Plant 

12  Hp  /  2  Cycle 

Wing  Airfoil 

NACA2415 

Horiz.  Stab.  Airfoil 

NACA  0006  (Approx.) 

^w  (^ref) 

17.57  ft2               2530  in2 

st 

3.174  ft2             457.1  in2 

Sv 

0.9818  ft2              141.4  in2 

c 

1 .66  ft                   20  in 

c, 

0.958  ft                 11.5  in 

bw 

10.54  ft               126.5  in 

b, 

3.313  ft               39.75  in 

bv 

1.25  ft                 15.0  in 

It 

4.44  ft               53.25  in 

lv 

4.44  ft               53.25  in 

ARW 

6.32 

ARt 

3.46 

ARV 

1.59 

VH 

0.49 

vv 

0.16 

Table  4.1   FROG  UAV  Characteristics,  after  Ref.  [2]. 
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Figure  4.1  FROG  UAV  Three- View  Drawing. 
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The  FROG  UAV,  as  operated  by  NPS,  is  equipped  with  airspeed,  angle-of-attack, 
altitude  and  control  surface  sensors.  In  addition,  a  miniature  Inertial  Measurement  Unit 
(IMU)  captures  aircraft  attitude,  acceleration  and  body  rates.  Data  is  down  linked  to  a 
mobile  SGI  workstation  through  a  spread  spectrum  modem.  Onboard  GPS  provides 
differential  GPS  navigation  capability  with  the  ground  station  used  as  a  reference.  The 
aircraft  can  be  flown  by  conventional  radio  control  or  by  up-linking  flight  control 
commands  from  the  SGI  workstation. 

Current  flight  control  development  requires  a  linearized  aerodynamic  model 
around  the  cruise  trim  point  of  60  m.p.h.  or  88  ft/s.  This  flight  condition  is  selected  for 
the  development  of  stability  derivative  data  with  CMARC.  Table  4.2  lists  the  aircraft 
parameters  for  the  trim  flight  condition. 


PARAMETER 

MEASUREMENT 

UNITS 

Weight 

67.73 

lbs 

IXX 

12.52 

slug-ft 

IYY 

8.43 

slug- ft 

IZZ 

18.55 

slug-ft 

Airspeed 

60/88 

mph  and  ft/s 

Altitude 

800 

ftMSL 

Air  Density 

0.002327 

slug/ft3 

Center  of  Gravity 

34.5% 

M.A.C 

*-  L  trim 

0.4295 

n/a 

*-*•  trim  (est) 

-0.04 

degrees 

°  Etnm 

5.1 

degrees 

Table  4.2  FROG  UAV  Trim  Flight  Condition,  after  Ref.[2]. 
C.         FROG  UAV  MODELING 

1.    General 

LOFTSMAN  is  utilized  for  the  creation  of  all  CMARC  input  file  patches  except 
for  wing  tips.  In  some  cases,  CMARC's  more  efficient  built-in  capability  to  model 
standard  NACA  4-digit  wing  surfaces  could  be  used.  However,  with  the  requirement  to 
mesh  deflected  control  surfaces,  all  patches  are  created  with  LOFTSMAN  from  the  start. 
The  complete  FROG  UAV  model  with  all  patches  activated  is  displayed  in  Figure  4.2. 
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Appendix  A  contains  detailed  descriptions  of  obtaining  stability  derivative  data  in  a  less 
formal  user  2uide  format. 


Figure  4.2  FROG  UAV  Panel  Code  Model  with  All  Patches  Activated. 

Some  assumptions  are  made  to  simplify  the  modeling  process.  First,  the 
horizontal  and  vertical  stabilizers  are  modeled  with  a  NACA  0006  section.  The  actual 
surfaces  are  constructed  with  a  flat  section,  rounded  at  the  nose  and  tapered  starting  at  the 
control  surface  hinge  line  to  a  sharp  trailing  edge.  The  NACA  0006  provides  a  close 
approximation  and  allows  the  use  of  LOFTMAN's  built-in  wing  lofting  capability.  For  a 
potential  flow  solution,  this  simplification  is  considered  minor. 

A  second  simplification  is  made  regarding  the  vertical  stabilizer's  tip  rib 
orientation.  The  actual  rib  is  canted  down  5°  with  respect  to  the  longitudinal  waterline. 
LOFTSMAN  will  only  model  a  chord  line  that  is  parallel  to  the  waterline  (constant  BL). 
Thus,  for  this  study,  the  vertical  tail  tip  rib  is  modeled  with  a  constant  BL,  but  the  span  is 
adjusted  to  maintain  the  same  overall  surface  area. 

A  third  assumption  is  made  regarding  aileron  control  surface  hinge  lines.  The 
actual  FROG  UAV  ailerons  are  piano  hinged  along  the  upper  surface.  For  simplification 
of  the  LOFTSMAN  model,  all  surfaces  are  hinged  at  the  centerline. 
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Finally,  there  is  no  attempt  to  model  the  tricycle  landing  gear  struts  or  wheel 
assemblies.  The  landing  gear  components  do  not  contribute  significantly  to  the 
aerodynamic  stability  derivatives.  However,  they  certainly  need  to  be  taken  into  account 
when  measuring  moments  of  inertia  for  a  dynamic  model. 

2.  Modeling  Coordinate  System 

The  model  is  developed  using  a  coordinate  system  selected  to  simplify  fuselage 
measurements.  The  +x-axis  starts  even  with  the  nose  and  runs  aft  along  the  bottom  of  the 
fuselage,  parallel  with  the  tail  boom.  The  bottom  of  the  fuselage  is  used  as  the  waterline 
with  +z-axis  in  the  up  direction.  This  allows  for  easy  vertical  measurements  when  the 
aircraft  is  placed  flat  on  a  horizontal  surface.  The  +y-axis  runs  from  centerline  outboard 
parallel  to  the  right  wing.  Figure  4.3,  which  displays  static-pressure  source  and  alpha 
vane  locations,  also  shows  the  location  and  origin  of  the  modeling  coordinate  system. 

3.  LOFTSMAN  Patches 

LOFTSMAN  is  used  to  generate  all  the  model  patches  except  for  wing  tips. 
CMARC's  built-in  capability  is  used  to  create  wing  tip  patches.  Appendix  B  contains 
listings  of  all  the  LOFTSMAN  input  files.  Once  a  surface  is  meshed,  the  mesh  is  saved 
to  a  file  as  a  CMARC/PMARC  patch.  The  resulting  text  file  is  then  opened,  and  the  text 
is  copied  and  pasted  with  any  text  editor  into  the  patch  definition  section  of  the  CMARC 
input  file.  LOFTSMAN  patch  files  are  not  listed  because  they  are  redundant  with  the 
patches  in  the  final  CMARC  input  file  listed  in  Appendix  C.  Multiple  patches  for  the 
same  control  surface  may  be  stored  in  the  input  file.  For  instance,  two  horizontal 
stabilizer  patches  with  different  elevator  deflection  angles  may  be  kept  in  the  CMARC 
input  file.  However,  only  one  may  reside  in  the  active  patch  input  section.  All  other 
patches  for  the  same  aircraft  component  must  reside  after  the  patch  with  the  "last  patch" 
toggle  set  (TNODS=5). 

When  saving  a  patch,  LOFTSMAN  automatically  takes  care  of  all  CMARC  input 
formatting.  A  patch,  as  formatted  by  LOFTSMAN,  assumes  additional  patches  will 
follow  in  the  CMARC  input  file.  Therefore,  the  last  segment's  TNODS  variable  is  set 
TNODS=3.  When  the  patch  is  the  last  patch  in  the  input  file,  the  TNODS  variable  must 
be  manually  set  to  TNODS=5.  If  CMARC  hangs  up  while  reading  in  geometry 
information,  most  likely  TNODS=5  is  missing  on  the  last  patch. 
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a.  Fuselage  Model 

The  fuselage  is  lofted  as  a  B-type  body.  A  B-type  body  is  used  when 
major  portions  of  the  fuselage  have  a  circular  or  oval  cross  section.  The  input  file  is 
listed  in  Appendix  B.  Only  the  right  side  is  meshed,  with  a  symmetric  left  side  created 
by  toggling  the  EPATSYM  variable  to  IPATSYM=1.  LOFTSMAN  assumes  that  B-type 
bodies  converge  to  a  specific  point  at  the  fore  and  aft  ends.  The  flat  aft  fuselage  face  does 
not  provide  this  single  point.  A  slight  modification  was  made  to  the  aft  face  to  allow 
automatic  meshing  as  a  B-type  body.  The  center  of  the  aft  face  is  extended  very  slightly, 
approximately  1/8  inch,  to  provide  a  convergence  point  for  the  final  rear  triangular 
panels.  This  small  deviation  is  assumed  not  affect  the  aerodynamic  fidelity  of  the  model 
for  a  potential  flow  solution. 

The  right  side  was  originally  meshed  separately  from  the  wing  as  a  20  x 
20  panel  patch.  This  created  a  low  order  fit  when  the  wing  patch  was  butted  to  the  side  of 
the  fuselage,  resulting  in  overlapping  panels.  A  final  mesh  was  created  that  flowed 
around  the  wing  root  and  fuselage  intersection  for  a  high  order  fit.  All  the  fuselage 
panels  at  the  wing  root  join  with  the  adjacent  wing  panels.  This  mesh  requires  that  the 
fuselage  be  broken  up  into  six  separate  panels  per  side.  They  are  the  nose  patch,  the 
forward  transition  patch,  the  top  and  bottom  wing  root  patches,  the  aft  transition  patch 
and  finally  the  rear  fuselage  patch.  Some  manual  editing  is  required  to  straighten  out 
panels  on  the  upper  fuselage  patch.  When  the  six  patches  are  added  together,  the  final 
configuration  is  modeled  with  a  44x15  panel  patch. 

b.  Main  Wing  Patch 

The  NACA  2415  wing  is  created  with  four  separate  patches  to  allow  the 
addition  of  a  deflected  aileron  mesh.  CMARC  comes  with  a  broad  selection  of  "*.SD" 
airfoil  template  files  that  are  automatically  loaded  during  installation.  The 
"NACA2415.SD"  file  is  used  for  this  model.  The  10x15  inboard  wing  patch  runs  from 
the  wing  root,  past  the  flaps,  to  the  start  of  the  aileron.  The  10x15  mid- wing  patch  covers 
the  portion  of  the  wing  spanned  by  the  aileron.  The  tapered  outboard  wing  extension  is 
made  with  a  6x15  patch.  Finally,  a  4x15  semi-circular  wing  tip  patch  is  created  using 
CMARC's  built-in  wing  tip  functionality.  The  wing  is  set  to  a  4.5°  incidence  in  the 
LOFTSMAN  input  file.  Alternatively,  the  patch  could  be  created  with  no  incidence  and 
then  the  patch  coordinate  system  could  be  rotated  in  the  CMARC  input  file.  Together, 
the  three  wing  patches  add  to  make  a  30x15  panel  wing  model. 
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The  ailerons  are  meshed  at  zero  degree  deflection  and  five  degrees  of  right 
rolling  moment  deflection.  For  the  no  deflection  case,  a  single  right  wing  patch  is 
meshed  with  CMARC's  symmetric  patch  toggle  (IPATSYM=1)  creating  the  left  patch. 
For  deflected  ailerons,  separate  patches  are  meshed  for  each  wing.  The  right  wing  aileron 
deflection  is  five  degrees  trailing  edge  up  (T.E.U.)  and  the  left  wing  five  degrees  trailing 
edge  down  (T.E.D.). 

c.  Horizontal  Stabilizer  Patch 

The  horizontal  stabilizer  and  elevator  patch  is  created  with  a  single  10x16 
mesh  using  the  "NACA0006.SD"  airfoil  template.  As  stated  earlier,  the  actual  UAV 
horizontal  stabilizer  has  flat  upper  and  lower  surfaces  and  a  rounded  leading  edge.  The 
NACA  0006  airfoil  is  substituted  because  it  provides  a  close  approximation  to  the  leading 
edged  radius  and  thickness  ratio.  No  other  special  modifications  are  required.  A  tip 
patch  is  not  added  because  some  of  the  resulting  panels  would  be  too  small.  In  particular, 
the  triangular  panels  closing  out  the  aft  end  of  the  tip  are  too  small  in  proportion  to  the 
other  panels,  creating  singularities  due  to  machine  resolution.  Initially,  an  attempt  was 
made  to  model  horizontal  and  vertical  stabilizer  wing  tips,  but  the  model  would  not 
converge  with  them.  Leaving  off  tip  patches  does  not  significantly  influence  results 
according  to  the  CMARC  User's  Guide  [Ref.  3]. 

The  elevator  spans  the  entire  horizontal  stabilizer.  Two  patches  are 
meshed,  one  with  zero  deflection  and  the  second  with  positive  five  degrees  (+5°  T.E.D.) 
of  deflection.  The  deflected  patch  is  used  to  obtain  the  elevator  control  power  derivatives 
(Clso  CM8e  and  CoSe)- 

d.  Vertical  Stabilizer  Patch 

Two  different  versions  of  the  vertical  stabilizer  are  used  for  FROG 
modeling.  An  extended  vertical  stabilizer,  which  includes  the  "effective"  area  of  the 
empennage  tail  boom,  is  used  to  measure  the  static  sideslip  and  yaw  rate  damping 
derivatives.  It  is  a  10x14  panel  patch.  The  second  patch  models  the  actual  area  of  the 
vertical  stabilizer.  It  is  a  10x12  patch  and  is  meshed  with  a  five-degree  rudder  deflection 
for  measuring  control  power  derivatives. 

The  vertical  stabilizer  and  rudder  patches  are  created  with  a  single  mesh 
using  the  "NACA0006.SD"  airfoil  template.  As  with  the  horizontal  tail  stabilizer,  the 
NACA  0006  airfoil  closely  approximates  the  vertical  surface  leading  edge  radius  and 
thickness  ratio.    The  LOFTSMAN  input  file  is  different  in  that  a  vertical  wing  surface 
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requires  a  modification  to  the  rib  axis.  The  rib  axis  must  be  specified  with  an  x-axis 
rotation  of  90°,  a  y-axis  rotation  of  0°  and  an  unspecified  (999.0)  z-axis  rotation.  No 
symmetry  is  selected  for  the  vertical  stabilizer  because  the  patch  is  already  symmetric 
about  the  y=0  plane.  As  with  the  horizontal  stabilizer,  a  tip  patch  is  not  added  because 
some  of  the  resulting  panels  would  be  too  small. 

The  rudder  spans  the  entire  vertical  stabilizer.  Two  patches  are  meshed, 
one  with  zero  deflection  and  the  second  with  positive  five  degrees  (+5°  T.E.L.)  of 
deflection.    The  deflected  patch  is  used  to  obtain  the  elevator  control  power  derivatives 

(ClSc  CM6e  and  Co6e)- 

e.     Tail  Boom  Patch 

The  tail  boom  patch  is  created  as  a  single  12  x  10  mesh  using  a  B-type 
body.  Again,  only  the  right  side  is  meshed  due  to  symmetry.  The  LOFTSMAN  input  file 
requires  modifications  at  both  ends  in  a  similar  fashion  to  the  aft  fuselage.  A  single  point 
is  added  to  allow  convergence  of  the  triangular  panels  at  either  end.  With  this  point,  the 
tail  boom  has  the  appearance  of  being  tapered  at  both  ends.  The  point  is  then  manually 
edited  out  in  the  CMARC  input  file  by  replacing  the  "x"  coordinate  of  the  beginning  and 
ending  section  panels  with  the  correct  value.  In  most  cases,  the  tail  boom  is  left  out  of 
solution  to  aid  in  convergence.  This  is  due  to  the  small  overlapping  panels  at  the  fuselage 
tail  boom  junction.  Being  a  slender,  round  tube  directly  in  the  fuselage  slipstream,  the 
tail  boom  is  assumed  to  have  a  small  influence  on  the  stability  derivatives. 

/    Engine  Pod  Patch 

The  engine  pod  patch,  or  nacelle,  is  created  as  a  single  15x10  mesh  using 
a  B-type  body.  Only  the  right  side  is  meshed  due  to  symmetry.  The  prop  spinner  is  an 
integral  part  of  the  patch.  No  attempt  is  made  to  model  the  prop,  engine  heads  or  exhaust 
system.  As  with  the  tail  boom,  in  most  cases  the  engine  pod  is  left  off  the  model  used  to 
establish  stability  derivatives. 

g.    Engine  Pylon  Patch 

The  engine  pylon  patch  is  modeled  with  a  single  15x10  mesh  using  an  A- 
type  body.  A-type  bodies  are  used  to  model  surfaces  similar  to  boat  hulls  with  cornered 
surfaces  or  sharp  chines.  In  addition,  A-type  bodies  do  not  require  the  body  to  be 
completely  enclosed.  As  a  result,  an  A-body  was  selected  to  model  just  the  sides  of  the 
pylon.  Only  the  right  side  is  meshed  due  to  symmetry.  A  low  order  fit  is  achieved  with 
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the  adjacent  fuselage  and  engine  pod  panels.  This  results  in  questionable  pressure 
distributions.  As  a  result,  the  pylon  patch  was  turned  off  for  most  configurations.  A 
future  attempt  could  be  made  to  create  a  high  order  fit  between  the  other  patches  or  to 
model  the  pylon  as  a  flat  surface.  This  will  probably  require  manual  editing  of  the 
intersecting  patches.  Even  with  a  high  order  fit,  a  wake  cannot  be  added  to  the  trailing 
edge  of  the  pylon.  It  would  impact  the  vertical  tail  surface 

4.    Common  CMARC  Input  File  Errors 

The  patches  created  in  LOFTSMAN  are  assembled  into  a  single  CMARC  input 
file  with  any  text  editor.  The  default  input  file,  which  comes  with  CMARC,  or  any  old 
input  file  may  be  modified.  There  are  many  errors  that  will  cause  CMARC  to  hang  up 
without  an  error  message.  The  two  most  common  errors  are  forgetting  to  designate  the 
last  patch  and  incorrectly  numbering  the  wake  patches. 

The  last  patch  must  be  designated  by  including  a  TNODS=5  setting  in  the  last 
section  of  the  last  patch.  If  it  is  not  included,  CMARC  hangs  up  when  reading  in  the 
geometry.  In  a  similar  manner,  the  last  wake  must  be  designated  with  a  NODEW=5 
setting.  If  the  last  wake  is  not  designated,  CMARC  hangs  up  while  reading  in  the  wake 
information.  Another  common  error  involves  incorrect  wake  to  patch  number 
association.  Patch  numbering  changes  whenever  patches  are  disabled  or  reordered.  The 
KWPACH  field  for  each  wake  definition  must  be  checked  to  make  sure  it  reflects  the 
current  patch  numbering. 

D.         STATIC-PRESSURE  SOURCE  AND  ALPHA  VANE  CORRECTIONS 
THROUGH  OFF-BODY  FLOW  ANALYSIS 

CMARC  is  ideally  suited  for  off-body  flow  analysis.  Off-body  streamlines  may 
be  placed  through  a  point  anywhere  in  the  flow  field.  CMARC  will  then  follow  the 
streamline  up  and  downstream  the  distance  designated  in  the  input  file.  This  is 
particularly  useful  for  flow  visualization.  In  addition,  CMARC  calculates  pressure 
coefficient  and  velocity  at  each  point  along  the  streamline.  For  this  study,  two 
streamlines  are  placed  through  the  locations  of  the  static-pressure  source  and  alpha  probe 
locations.  Pressure  coefficient  is  used  to  quantify  static  source  position  error  and  velocity 
is  used  to  calculate  alpha  probe  position  error  as  a  function  of  FROG  UAV  angle-of- 
attack.  Both  static  pressure  and  AOA  are  digitized  for  down  link  to  the  ground  station 
allowing  the  values  to  be  easily  corrected.  Either  a  look-up  table  or  curve  fit  correction 
can  be  applied  subsequent  to  being  passed  to  the  flight  control  routines.    This  analysis 
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was  previously  published  by  the  author  in  Ref.  [1],  but  is  included  here  to  provide  a 
single  source  of  UAV  modeling  conducted  to  date. 

1.  Description  of  the  FROG  UAV  Pitot-Static  and  AOA  Systems 

The  pitot-static  system  and  angle-of-attack  probe  share  a  common  flight  test 
boom  extending  from  the  nose  of  the  UAV.  The  boom  contains  both  the  total  and  static 
pressure  ports.  Figure  4.3  depicts  the  general  dimensions  of  the  flight  test  boom 
installation  and  the  modeling  coordinate  system. 

2.  Modeling  Off-Body  Streamlines 

Streamlines  are  placed  at  the  two  locations  indicated  in  Figure  4.3,  which 
correspond  to  the  static  source  and  alpha  probe  locations.  Two  off  body  streamlines  were 
activated  in  CMARC  by  setting  NSTLIN=2  in  the  &SLIN1  line.  Only  a  short  distance  of 
2  inches  is  selected  up  and  downstream  in  the  SU  and  SD  fields  to  reduce  the  size  of  the 
output  file.  Figure  4.4  is  a  POSTMARC  rendering  of  the  two  off-body  streamlines  used 
for  sensor  corrections.  With  the  model  at  at=0°,  notice  that  the  streamline  is  curving  up 
at  the  angle-of-attack  vane  location  6.5  inches  in  front  of  the  aircraft  nose. 

3.  Analysis  of  Static  Source  Position  Errors 

The  position  error  pressure  coefficient,  ACp  pc  or  APp/qc,  is  a  function  of  free 
stream  Mach  number  and  angle-of-attack  provided  that  the  static  source  is  located  outside 
of  a  thick  boundary  layer  and  sideslip  is  minimized  [Ref.  11].  In  the  case  of  the  FROG 
UAV  with  incompressible  flow,  APp/qc  becomes  a  function  of  angle-of-attack  only. 
Therefore,  corrections  can  be  simply  defined  as  a  function  of  measured  angle-of-attack. 

A  DOS  batch  file  was  executed  to  step  the  CMARC  model  through  angles-of- 
attack  ranging  from  -8°  to  20°.  The  batch  file  incremented  the  angle-of-attack  using 
CMARC  s  command  line  override  feature.  In  addition,  a  new  output  file  name  was 
designated  for  each  angle-of-attack.  Position  error  pressure  coefficient  is  then  read  from 
the  off-body  streamline  listing  of  the  output  file  at  the  location  corresponding  to  the  static 
source.  Table  4.3  lists  the  values  of  APp/qc  calculated  from  CMARC  data.  Figure  4.5 
displays  ACp  pc  as  a  function  of  indicated  angle-of-attack.  The  second  order  influence  of 
angle-of-attack  is  clear  with  a  second  order  curve  fitting  tightly  through  the  data  points. 
Of  note,  the  error  is  relatively  constant  for  a  ±8°  band  around  trim  angle-of-attack.  For 
incompressible  flow,  position  error  pressure  coefficient  is  independent  of  airspeed  and 
altitude. 
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Figure  4.3  Diagram  of  the  FROG  UAV  Pitot-Static  and  AOA  Systems. 
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Figure  4.4  FROG  Off-Body  Streamline  Visualization  with  POSTMARC  (ct=10°). 

Position  error  pressure  coefficient  can  be  turned  into  position  corrections  for 
airspeed  and  altitude.  The  following  relations  were  developed  which  assume  small  errors 
and  incompressible  flow: 


AVpc  = 


VxACp 


and      AVpc^Vc-Vi 


4.1 


AHpc  = 


Where: 


CTstdSo 


and      AHpc  =HC  -  Hj 


AHpc  is  the  altitude  position  correction. 

AVpc  is  the  velocity  position  correction. 

APD 
ACD  =  — —  or  position  error  pressure  coefficient. 

qc 

cistcj  is  standard  day  density  ratio. 
g0  is  the  gravitational  constant. 


4.2 
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Table  4.3  displays  corrections  calculated  for  both  airspeed  and  altitude  at  the 
FROG  UAV  trim  condition  of  88  ft/s  and  800  ft  MSL.  The  corrections  are  added  to  the 
indicated  value  to  obtain  the  corrected  value.  Figures  4.6  and  4.7  display  the  corrections 
as  a  function  of  indicated  angle-of-attack.  Again,  a  second  order  curve  fits  nicely  through 
the  data  points.  Equations  4. 1  and  4.2  can  be  used  to  implement  a  correction  algorithm 
based  on  airspeed  and  altitude. 

4.  Analysis  of  Alpha  Vane  Position  Error 

Local  flow  field  velocity  is  extracted  from  the  off-body  streamline  listing  to 
obtain  local  angle-of-attack.  The  alpha  vane  is  assumed  to  capture  the  x-z  component  of 
the  local  velocity  field  and  ignore  cross  flow  in  the  y  direction.  Flow  field  velocity  is 
turned  into  indicated  angle-of-attack  and  angle-of-attack  position  correction  with  the 
following  equations: 


a  i  °  =  a  tan 


vz 


* 


180 


TC 


degrees  4.3 


Act  pc  =  ct  t  -  a  j  degrees  4.4 


A  DOS  batch  file  is  executed  to  step  the  CMARC  model,  with  an  off-body 
streamline  located  at  the  vane  position,  through  angles-of-attack  ranging  from  -8°  to  20°. 
Local  velocity  components  are  then  read  from  the  location  corresponding  to  the  alpha 
vane.  Table  4.4  lists  the  values  of  Aapc  calculated  from  CMARC  data.  Figure  4.8 
displays  Aa  pc  as  a  function  of  indicated  angle-of-attack.  Linear  and  second  order  curve 
fit  equations  are  also  indicated  on  Figure  4.8.  Angle-of-attack  correction  is  fairly  linear 
through  the  FROG  operating  envelope,  with  approximately  -1.25  degrees  of  position 
error  at  the  FROG  cruise  trim  condition.  The  corrections  apply  at  all  incompressible 
airspeeds  and  all  altitudes. 

5.  Summary  of  Off-Body  Flow  Field  Analysis 

CMARC  proved  useful  for  both  static-pressure  source  and  alpha  vane  position 
corrections.  Measured  data  may  be  corrected  using  look-up  tables  with  the  values  in 
Table  4.3  and  4.4  or  by  using  the  curve  fits  in  Figures  4.5  through  4.8.  Flight  testing  is 
recommended  for  validation  of  sensor  corrections  obtained  from  this  CMARC  off-body 
flow  field  analysis. 
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2.8 
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18 
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2.3 
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20 
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1.8 

5.1 

Table  4.3       Position  Error  Corrections  for  the  NPS  FROG  UAV  at  V=88  ft/s  and 
H=800  ft  MSL.  Derived  from  CMARC  Panel  Code  Off-Body  Flow 
Field  Analysis. 
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Source:  CMARC  Panel  Code 

Trim  Condition:  V=88  ft/s    H=800  ft  MSI 
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Figure  4.5      Position  Error  Pressure  Coefficient,  ACP  p„  for  the  NPS  FROG  UAV. 
Derived  from  CMARC  Panel  Code  Off-Body  Flow  Field  Analysis. 
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Source:  CMARC  Panel  Code 
Trim  Condition  V=88  ft/s  H=800  ft  MSL 
Flight  Test  Boom  Static  Source  (8.25"  torward  of  nose) 
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AOA   Indicated  (deg) 


Figure  4.6      Altitude  Position  Error,  AHpc,  for  the  NPS  FROG  UAV  at  V=88  ft/s 
and  H=800  ft  MSL.  Derived  from  CMARC  Panel  Code  Off-Body 
Flow  Field  Analysis. 
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Source:  CMARC  Pane!  Code 
Trim  Condition:  V=88  ft/s    H=800  ft  MSL 
Flight  Test  Boom  Static  Source  (8.25"  forward  of  nose) 


0  4  8  12  16 

AOA   Indicated  (deg) 


Figure  4.7      Airspeed  Position  Error,  AVpc,  for  the  NPS  FROG  UAV  at  V=88  ft/s 
and  H=800  ft  MSL.  Derived  from  CMARC  Panel  Code  Off-Body 
Flow  Field  Analysis. 
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UAV  AOA 

aT  (deg) 

Velocity  at  Alpha  Vane 

AUA  Correction 

Aa=aT-otl  (deg) 

AOA  indicated 

a;  (deg) 

Vx  (ft/s) 

Vy  (ft/s) 

Vz  (ft/s) 

-8 

80.92 

1.66 

-12.23 

0.60 

-8.60 

-6 

81.27 

1.65 

-8.65 

0.08 

-6.08 

-3 

81.60 

1.64 

-3.21 

-0.75 

-2.25 

-2 

81.67 

1.63 

-1.47 

-0.97 

-1.03 

-1 

81.71 

1.63 

0.28 

-1.20 

0.20 

0 

81.73 

1.62 

2.13 

-1.49 

1.49 

1 

81.73 

1.61 

3.93 

-1.75 

2.75 

2 

81.70 

1.60 

5.72 

-2.00 

4.00 

3 

81.66 

1.59 

7.51 

-2.25 

5.25 

4 

81.58 

1.58 

9.30 

-2.50 

6.50 

5 

81.48 

1.57 

11.08 

-2.75 

7.75 

6 

81.37 

1.56 

12.88 

-2.99 

8.99 

8 

81.07 

1.53 

16.43 

-3.46 

11.46 

10 

80.67 

1.51 

19.98 

-3.91 

13.91 

12 

80.17 

1.48 

23.50 

-4.34 

16.34 

14 

79.61 

1.46 

26.99 

-4.73 

18.73 

16 

78.93 

1.43 

30.47 

-5.11 

21.11 

18 

78.18 

1.39 

33.90 

-5.44 

23.44 

20 

77.34 

1.36 

37.31 

-5.75 

25.75 

Table  4.4       Angle-of  Attack  Vane  Position  Error  Corrections  for  the  NPS  FROG 
UAV.  Derived  from  CMARC  Panel  Code  Off-Body  Flow  Field 
Analysis. 
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Source:  CMARC  Panel  Code 
Trim  Condition:  V=88  tt/s  H=800  ft  MSL 
Flight  Test  Boom  Alpha  Vane  (6.5"  forward  of  nose) 


4  8  12  16 

AOA  Indicated  (deg) 


Figure  4.8      Angle-of-Attack  Vane  Position  Error,  Aapc,  for  the  NPS  FROG  UAV. 
Derived  from  CMARC  Panel  Code  Off-Body  Flow  Field  Analysis. 
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E.  DEVELOPMENT  OF  STABILITY  DERIVATIVES 

In  this  section,  CMARC  is  used  to  develop  the  longitudinal  and  lateral-directional 
stability  derivatives  for  the  FROG  UAV.  The  development  effort  focuses  on  finding  the 
static,  rate-damping  and  control-power  derivatives.  The  results  obtained  from  CMARC 
are  compared  to  data  sets  obtained  from  empirical  estimation  techniques  and  flight-test 
parameter  estimation.  In  the  next  section,  the  stability  derivatives  are  entered  into  a 
dynamic  model  to  find  modal  frequency,  damping  and  response  to  control  deflections. 

The  potential  flow  analysis  performed  by  CMARC  does  not  provide  accurate 
viscous  drag  values.  Therefore,  total  drag  is  estimated  from  the  flight-test  power-off 
glide  ratio  and  cruise  thrust  required.  The  User  Guide  developed  in  Appendix  A 
describes  in  detail  how  to  gather  stability  derivative  data  with  CMARC.  An  abbreviated 
discussion  is  presented  below. 

CMARC  contains  built-in  functionality  to  integrate  forces  and  moments  in  all 
axes  over  the  surface  of  a  body.  Force  and  moment  coefficients  are  automatically 
nondimensionalized  based  on  the  mean  aerodynamic  chord,  reference  wing  area,  semi- 
span  and  center  of  gravity  location  in  the  CMARC  BINP9  input  line.  Coefficients  are 
presented  in  both  wind  and  body  axes.  Of  note,  CMARC  uses  the  semi-span  to 
nondimensionalize  rolling  and  yawing  moment  coefficients.  Most  texts  on  stability  and 
control,  including  Roskam  [Ref.  12]  and  Etkin  [Ref.  13],  nondimensionalize  rolling  and 
yawing  moments  by  span.  Rolling  and  yawing  moments  are  nondimensionalized  by  span 
in  this  study.  Table  4.5  summarizes  the  factors  for  normalizing  moments  and  angular 
rates.  All  rolling  and  yawing  moment  coefficients  presented  in  this  study  have  been 
normalized  with  span  by  dividing  the  CMARC  output  by  a  factor  of  two.  Table  4.5  also 
indicates  the  characteristic  time,  t*,  employed  for  angle  rate  data  reduction. 

In  addition  to  differences  in  nondimensionalizing  moments,  CMARC  uses  the 
typical  CFD  axis  system  shown  in  Figure  4.9.  For  this  study,  all  work  is  performed  in  the 
stability  axis  system.  Figure  4.9  also  illustrates  the  standard  stability  axes  implemented 
in  this  study.  The  sign  of  CMARC  roll  and  yaw  moments  need  to  be  reversed.  The 
direction  for  positive  control  deflections  is  also  shown  in  Figure  4.9.  All  control  surfaces 
are  patched  with  positive  deflections  using  the  convention  in  Figure  4.9. 
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A  potential  flow  solution  will  not  produce  satisfactory  results  for  bodies  with 
significant  areas  of  flow  separation.  Therefore,  CMARC  models  must  be  analyzed  in  the 
linear  slope  regions  for  valid  results. 


MOMENTS 

NORMALIZING 
PARAMETER' 

RATES 

CHARACTERISTIC 
TIME 

L  -  CtqSb 

b 

pb 

P=  i 

2u0 

/2Uo 

M  =  CmqSc 

c 

re 
r  = 

2u0 

t*=  c/ 

/2uo 

N  =  CrqSb 

b 

„      rb 

r  — 

2u0 

/2uo 

Note:  1)  CMARC  nondimensionalizes  roll  and  yaw  coefficients  with  b/2. 
Table  4.5  Nondimensionalized  Moment  and  Rate  Equations. 

For  the  static  stability  derivatives,  the  CMARC  model  is  run  at  two  different 
angles-of-attack  and  one  sideslip  angle  with  zero  control  surface  deflections.  The  slopes 
of  the  force  and  moment  coefficients  are  then  taken  to  produce  the  Clo  and  Cma 
longitudinal  derivatives  and  the  Cyp,  Qp  and  Cnp  lateral-directional  derivatives. 

For  the  control-power  derivatives,  the  model  is  run  at  the  trim  condition  with 
successive  control  deflections.  The  difference  between  the  results  with  and  without 
control  surface  deflections  yield  the  Cu>t,  Cmsc  and  Co8e  longitudinal  and  the  Cy8r,  Cigr, 
Cn5r,  Ci6a  and  Cn5a  lateral-directional  control-power  derivatives. 

Development  of  the  damping  derivatives  is  not  as  straight  forward.  The  static 
derivatives  are  obtained  with  motion  disabled.  For  the  longitudinal  damping  derivatives, 
the  model  is  run  with  oscillating  vertical  plunging  motion  to  obtain  the  Cl  and  Cm  a-dot 
terms.  The  lift  and  pitching  moment  coefficients  are  broken  into  real  (in  phase  with 
AOA)  and  imaginary  (out  of  phase  with  AOA)  components.  The  imaginary  components 
are  due  to  a-dot  effects.  Next,  the  model  is  run  with  oscillating  pitch  motion  to  obtain 
the  combined  a-dot  and  pitch  rate  terms.  Subtracting  the  a-dot  influence  obtained  from 
the  plunging  motion  isolates  the  pitch  rate-damping  term  from  the  pitch  motion. 
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CMARC  Axis  System 


+$  =  Roll   Angle  ($=P) 
+L  ■  Roll   Moment 


+v  =  Lateral 
Vel. 
+Y  =  Side     /% '/V 

Force  <s  /    +6  =  Pitch  Angle   (8=q) 
+M  =  Pitch  Moment 


\> 


AN,   <  0 


+w  =  Vertical 

Vel. 
+Z  ■  Normal        f     ) 


^    AM.    < 


AM,   <   0 


Force 


2  " 


+*  =  Yaw  Angle  (f=r) 
>N  =  Yaw  Moment 


Stability  Axis  System 


Figure  4.9   CMARC  CFD  Axis  System  Compared  To  Stability  Axis  System.  From 
References  [3]  and  [14]  Respectively. 
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For  the  lateral-directional  analysis,  the  (3-dot  terms  are  assumed  to  be  small.  This 
allows  the  model  to  be  run  with  just  oscillating  roll  and  yaw  motion.  As  with  the 
longitudinal  test  case,  the  imaginary  or  out  of  phase  component  yields  the  combined  (3- 
dot  and  rate-damping  data.  With  the  (3-dot  terms  assumed  to  be  small,  the  oscillating 
motion  yields  the  Cy,  Q,  and  Cn  roll  and  yaw  rate  terms  directly. 

Throughout  this  chapter,  stability  derivatives  obtained  from  CMARC  are 
compared  to  three  other  data  sets.  The  first  comes  from  classical  analysis  through  the 
methods  presented  by  Roskam  in  References  [12]  and  [15].  The  second  set  comes  from  a 
classical  aerodynamic  model  developed  by  Papageorgio  [Ref.  2]  and  further  refined 
through  parameter  estimation  by  Engdahl.  The  resulting  parameter  estimation  model 
provides  a  good  match  to  observed  flight  characteristics  and  is  currently  in  use  by  the 
NPS  Department  for  Aeronautical  Engineering  for  closed-loop  flight  control 
development.  The  third  set  of  data  is  for  the  Cessna  1 72  at  the  cruise  flight  condition. 
The  data  is  only  presented  for  an  order  of  magnitude  comparison. 

1.    Longitudinal  Stability  Derivatives 

a.  Static  Longitudinal  Stability  Derivatives 

Three  basic  longitudinal  stability  derivatives  can  be  measured  with  just 
two  runs  of  the  CMARC  model.  The  model  is  first  analyzed  at  an  angle-of-attack 
corresponding  to  the  estimated  trim  condition.  In  this  case,  at=0°  is  selected  for  the  first 
run.  A  second  CMARC  run  is  conducted  with  angle-of-attack  incremented  by  one  or  two 
degrees.  Cl  and  Cm  are  then  extracted  manually  from  the  data  files.  The  slope  of  Cl  and 
Cm  versus  angle-of-attack  provide  the  Cl0  and  Cma  longitudinal  derivatives.  For  this 
study,  several  angles-of-attack  were  analyzed  to  check  consistency  of  the  slope.  In 
addition,  a^m  is  calculated  from  the  lift  curve  slope  and  trim  lift  coefficient.  For  the 
longitudinal  analysis,  only  half  the  model  is  analyzed.  The  symmetric  calculation  mode 
is  selected  by  setting  both  RSYM=0.0  and  IPATSYM=0  in  the  CMARC  input  file. 
Equations  4.5  through  4.7  are  used  for  these  calculations: 

(cl2 -Cl^iso  ..  ... 

CL    =—, H1    per  radian  4.5 

(a2-a1)        n 
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(cm2  -cmi)    180 


Cm     = — : r — * per  radian  4.6 

(a2-a1)         n 


cLtnm-cL]J     180 
a°trim=a°l  + ~ degrees  4.7 


Two  FROG  UAV  model  configurations  are  analyzed  in  a  build-up 
approach  to  check  results  against  classical  calculations  and  flight-test  data.  Figure  4.10 
shows  the  CMARC  models.  First,  just  the  wing  and  horizontal  tail  are  considered.  The 
patches  for  all  other  surfaces  and  wakes  are  turned  off  and  the  wing  root  is  extended  to 
centerline.  Then,  the  blended  wing/fuselage  and  horizontal  tail  are  analyzed.  Values  of 
Clo  and  Cma  for  these  two  configurations  are  presented  in  Table  4.6. 

Classical  design  calculations  are  also  performed  to  estimate  Cmafor 
comparison  to  CMARC  results.  Equation  4.8  is  used  for  the  calculation  of  Cma : 


^rrir,       aw 


(h-hac)-VH— (l-  — 
v  ac;  awV       da 


4.8 


In  classical  design,  the  horizontal  tail  downwash  derivative,  de/dot,  is  generally  selected 
from  empirical  data.  Using  a  taper  ratio  of  TR=1:1  and  an  aspect  ratio  of  AR=6, 
ds/da=0.4  is  selected  from  empirical  charts  in  Ref.  [16]  for  the  FROG  UAV 
configuration.  Also,  calculations  for  de/da=0.25  are  included  to  illustrate  pitching 
moment  sensitivity  to  the  downwash  derivative.  Table  4.6  lists  Cmafor  each 
configuration. 

The  final  static  longitudinal  parameter  required  is  total  aircraft  drag.  Drag 
coefficient  plays  an  important  role  in  long-period  aircraft  dynamics.  Unfortunately, 
potential  flow  panel  codes  such  as  CMARC  do  not  provide  accurate  total  drag  estimates. 
They  can  provide  good  induced  drag  predictions.  And,  if  equipped  with  a  boundary  layer 
code  like  that  contained  in  CMARC,  they  can  provide  integrated  skin  friction  results. 
However,  a  large  total  drag  contribution  in  the  form  of  pressure  drag  due  to  flow 
separation  is  not  accounted  for.    Total  drag  estimates  are  made  below  using  the  two 
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simple  techniques  shown  in  Equations  4.9  to  4.12.  The  first  method  is  based  on  flight- 
test  glide  ratio.  The  second  is  based  on  cruise  power  required  and  estimated  prop 
efficiency.  Note  that  the  selected  prop  efficiency  is  relatively  low  due  to  the  small 
propeller  diameter,  high  RPM  and  pusher  configuration.  The  two  methods  provide  drag 
predictions  within  10%  of  each  other.  The  results  are  averaged  to  Cd=0.065  and  included 
in  Table  4.6. 

Method  1 :  Lift-to-Drag  Ratio  (L/D=7  from  flight-test) 

L/=1    =>    p.k.E..™   fe=9.67    lbs  4.9 

/D  111 

r       D     9.67   lbs 

CD=  —  = =  0.0611  410 

^      0.5*0.002327    lb-s2 1  ft"  *882    ft2 1  s2  *  17.57    ft2 
Method  2:  Cruise  Power  Setting  (HP=5,  r|p=0.35) 


_D__Tr_,_  HPr*550    ftlbslslHP*rjpIV  4  { { 

qS      qS  qS 


C  5    ///>*550    ft-lbs/s/HP*0.35/SS    ft  I  s         =  Q  ^  4U 


0.5*0.002327    lb-s2/ft4*S82    ft2  /  s2  *  17.57   fr 


Average:    CDave  =(0.06 1 1  +0.069)/2  =  0.065 
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Simplified  Wing  and  Horizontal  Tail 


Blended  Wing  and  Fuselage  with  Horizontal  Tail 


Figure  4.10   CMARC  Models  of  the  FROG  UAV  for  the  Determination  of  Static 
Longitudinal  Stability  Derivatives. 
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Flight-test  data  for  the  short-period  and  phugoid  modes  were  used  for 
longitudinal  parameter  estimation.  Values  for  Clo  and  Cma  based  on  preliminary 
parameter  estimation  work  by  Engdahl  are  presented  in  Table  4.6. 


METHOD 

CONFIGURATION1 

STATIC  LONGITUDINAL  PARAMETERS 

Antrim               CLa               Cma                CD                CDa 

(aegj        (per  raa;      (per  raa;                        (per  raaj 

CMARC 
Panel  Code 

Wing/Horiz  Tail 

-0.87 

4.86 

-0.835 

n/a 

n/a 

Wing/Fuselage/Horiz  Tail 

-0.01 

4.85 

-0.413 

n/a 

0.266 

Classical 
Design3 

Wing/Horiz  Tail  -  5s/6a=0.25 

-0.81 

4.85 

-1.00 

n/a 

n/a 

Wing/Horiz  Tail  -  5s/5a=0.40 

-0.82 

4.82 

-0.70 

n/a 

0.253 

Parameter 
Estimation4 

Flying  Aircraft 

n/a 

4.09 

-0.42 

0.0655 

n/a 

C-1726 

n/a 

0.31 

-0.89 

0.31 

0.13 

NOTES:  1 )  CGX=34.5%  M.A.C.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Zero  lift  wing  incidence  is  +6.5°  from  the  longitudinal  reference  line. 

3)  Classical  design  after  Ref.  [12]. 

4)  Parameter  estimation  from  NPS  flight  test  data  by  Engdahl. 

5)  Average  drag  obtained  from  L/D  and  cruise  power  analysis  from  flight  test  data. 

6)  C-172  data  from  Ref.  [12]. 

Table  4.6  Comparison  of  FROG  UAV  Static  Longitudinal  Stability  Derivatives. 

Good  correlation  of  lift  curve  slope  (Cl<x)  to  classical  design  is 
demonstrated  by  the  two  CMARC  configurations  in  Table  4.6.  Both  techniques  over 
estimated  Clo  compared  to  flight-test  data.  Neither  the  classical  design  technique  nor  the 
CMARC  model  accounts  for  control  surface  leakage  or  separation  effects.  Modeling 
actual  flight  control  gap  geometry  may  provide  a  closer  match  to  experimental  data. 

Excellent  correlation  is  achieved  for  Cma  between  the  blended 
wing/fuselage  CMARC  model  and  flight-test  results.  The  low  value  of  Cma  implies  that 
ds/da  is  greater  than  the  0.4  value  selected  from  empirical  charts  in  Reference  [16].  The 
simplified  wing/tail  model  resulted  in  a  Cma  that  closely  matches  classical  design 
predictions  with  ds/da=0.35.  Clearly,  CMARC  modeling  provides  the  value-added 
capability  to  capture  fuselage  influences  that  are  missed  with  classical  design  methods. 
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An  investigation  was  performed  to  see  if  improved  wake  definition  could 
better  capture  the  ds/da  downwash  derivative.  Initially  a  flexible  wake  was  selected. 
However,  at  the  cruise  angle-of-attack,  the  flexible  wake  impacted  the  horizontal  tail 
causing  inaccurate  results.  A  further  investigation  was  made  into  a  streamline  based 
wake  definition.  A  study  by  Walden  et  al.  [Ref.  17]  studied  wake  turbulence  by 
modeling  an  aircraft  flying  in  trail  of  a  wake-generating  wing.  A  horizontal  tail  trailing 
the  main  wing  is  a  similar  configuration.  The  study  found  that  a  streamline-based  wake 
is  the  best  method  for  modeling  downwash  effects.  In  the  study,  the  model  was  initially 
run  with  a  rigid-wake.  A  streamline  was  also  defined  downstream  of  the  main  wing 
trailing  edge.  After  analyzing  the  first  results,  the  wake  was  then  predefined  to  follow  the 
streamline.  The  result  was  a  first  iteration  on  modeling  a  flexible  wake  with  a  predefined 
rigid-wake.  This  technique  was  tried  during  the  course  of  FROG  UAV  modeling. 
However,  FROG  geometry  proved  unsuited  to  this  technique.  The  streamline  trailing  the 
wing  ended  up  impacting  the  honzontal  tail.  This  indicated  that  a  predefined  wake  would 
also  impact  the  horizontal  tail,  resulting  in  inaccurate  result.  Figure  4. 1 1  shows  the 
streamline  impacting  the  horizontal  tail  at  cruise  angle-of-attack  with  rigid-wakes 
selected.  With  the  already  close  correlation  of  the  rigid-wake  results  to  flight-test  data, 
further  efforts  to  investigate  Cma  were  abandoned. 


Figure  4.11  Rigid-Wake  Streamline  for  Investigation  of  Streamline  Based  Wake 
Definition. 
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In  summary,  for  the  static  longitudinal  stability  derivatives,  CMARC 
produced  accurate  values  for  Cma  and  a  slightly  high  value  of  Cl<x-  Drag  coefficient  is 
obtained  by  averaging  the  values  from  the  lift-to-drag  ratio  and  cruise  thrust  required 
techniques. 


b. 


Longitudinal  Damping  Stability  Derivatives 


The  two  aircraft  motions  illustrated  in  Figure  4.12  are  used  to  develop  the 
longitudinal  rate-damping  derivatives.  A  sinusoidal  vertical  plunging  motion  isolates  the 
a-dot  effects.  Oscillatory  motion  is  controlled  with  the  CMARC  BINP8A  input  file  line. 
The  oscillating  pitch  angle  captures  both  the  a-dot  and  pitch  rate-damping  terms.  It  is 
controlled  with  the  CMARC  BINP8B  input  file  line.  The  a-dot  influence  developed 
from  the  plunging  motion  is  subtracted  from  the  pitching  motion  to  isolate  the  pitch  rate 
effect.  All  motion  is  conducted  at  a  frequency  of  2n  rad/s,  which  equates  to  a  reduced 
frequency  of  k=0.0595  for  this  configuration  and  trim  airspeed. 


Sinusoidal  Plunging  Motion 


Oscillating  Pitch  Motion 


Figure  4.12  Aircraft  Motion  for  the  Determination  of  a-dot  and  Pitch  Rate 
Damping. 
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The  two  models  shown  in  Figure  4.13  are  used  for  the  rate-damping 
investigation.  The  complete  wing/fuselage  horizontal  tail  model  is  used  for  both  the  oc- 
dot  and  pitch  rate-damping  terms,  while  the  simplified  horizontal  tail  only  model  is  used 
to  obtain  a  first  approximation  of  pitch  rate-damping.  The  User  Guide  in  Appendix  A 
provides  a  more  detailed  description  of  the  techniques  used  to  gather  the  dynamic 
longitudinal  stability  derivatives. 
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Blended  Wing  and  Fuselage  with  Horizontal  Tail 


Figure  4.13    CMARC  Models  of  the  FROG  UAV  for  the  Determination  of 
Longitudinal  Rate-Damping. 
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Traditional  design  calculations,  from  Roskam  [Ref.    12],  are  used  for 
comparison  to  CMARC  and  flight-test  results: 

Q0=2*Q„„*,7„*>V—  4.13 


/      ds 
C      =-2*C       *n    *V   -!-*——  414 

Lm°  La"    1h      H  c     da  .      4J4 


C,    =2*CLaH*r,H*VH  4.15 


r>       -—')*C         *  7?     *V    -±- 

^M.    ~      ^      ^LaH       'IH       y  H    — 


4.16 


where  the  dynamic  pressure  ratio  at  the  tail  is  taken  as  r|H=l  -0 

The  dynamic  longitudinal  stability  derivatives  are  presented  in  Table  4.7. 
CMARC-obtained  values  are  compared  to  classical  design  values  and  C-172  stability 
derivatives.  Flight-test  data  is  unavailable  for  the  damping  derivatives.  The  aerodynamic 
model  currently  in  use  by  the  FROG  research  team  implements  longitudinal  damping 
coefficients  obtained  from  classical  design  techniques.  The  CLa-dot  and  Cmq  obtained 
from  CMARC  for  the  complete  model  provided  a  good  match  to  the  classical  design  data. 
However,  both  CLa-dot  and  Cmq  are  over  estimated  by  approximately  50%  as  compared 
to  classical  design  values.  One  possible  explanation  for  the  additional  damping  from 
CMARC  is  that  the  classical  approximation  only  estimates  the  horizontal  tail  influence  on 
damping.  CMARC  also  measures  wing  and  fuselage  damping  influence.  The  horizontal 
tail  by  itself  produces  values  close  to  the  classical  calculations.  This  supports  the 
assertion  that  CMARC  also  captures  the  wing  and  fuselage  influence  for  the  whole 
model. 
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METHOD 

CONFIGURATION1 

DYNAMIC  LONGITUDINAL  PARAMETERS 

C  La-dot 

(aeg) 

cLq 

(per  raa; 

c 

^ma-dot 

(per  raa; 

Cmq 

(per  raaj 

Panel  Code 

Blended  Wing-Fuselage/Horiz  Tail 

1.42 

6.82 

-6.24 

-11.78 

Horizontal  Tail  only 

n/a 

4.37 

n/a 

-11.94 

Classical2 

Wing  /  Horizontal  Tail 

1.56 

3.89 

-4.14 

-11.39 

Par.  Est3 

Flying  Aircraft 

n/a 

n/a 

n/a 

n/a 

C-1724 

1.7 

3.9 

-5.2 

-12.4 

NOTES:  1 )  CGX=34.5%  M.A.C.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Classical  design  after  Ref.  [12]. 

3)  Parameter  estimation  from  flight  test  data  by  Engdahl. 

4)  C-172  data  from  Ref.  [12]. 

Table  4.7  Comparison  of  FROG  UAV  Dynamic  Longitudinal  Stability  Derivatives. 
c.  Longitudinal  Control-Power  Derivatives 


The  elevator  control-power  derivatives  are  obtained  by  substituting  a  0° 
deflection  horizontal-tail  patch  for  one  with  +5°  trailing  edge  down  deflection.  Only  one 
run  is  required.  The  difference  between  the  trim  condition  and  the  deflected  value  is 
divided  by  the  elevator  deflection  as  shown  below: 


[C,     -  Cr     ) 
C,„  -^*57.3    per    rad 


■L5c 


c»*  = 


Se2  -Sex 

V^Ds-  ~  CDSe.  ) 


Jde2 

5e2  -  5ex 


—  *57.3    per    rad 


C, 


\c 


m5e\  I 


lSe2 


lSe 


Sej  -5e\ 


57.3    per    rad 


4.17 


4.18 


4.19 


Elevator  control-power  derivatives  are  presented  in  Table  4.8.  A  classical 
design  estimate  of  elevator  control-power  is  also  shown  for  comparison.  The  following 
relationships  from  Roskam  [Ref.  15]  are  used  for  estimating  elevator  control-power: 
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C,       =CL        *  — 

Cm    =  -C,    *^-    or    Cm    = 

m5e  LSe  pr  mSe 


-CLs*V* 


4.20 
4.21 


where  Clsf  is  the  variation  of  lift  coefficient  with  elevator  deflection  found  from  charts  in 
Roskam[Ref.  15]. 

In  general,  the  elevator  control-power  derivatives  obtained  with  CMARC 
correlate  well  with  classical  design  techniques.  For  both  Cr^e  and  Cm5e,  CMARC  derived 
elevator  power  is  approximately  15  percent  higher  than  obtained  from  classical 
calculations.  However,  CMARC  significantly  under-estimates  elevator  control-power 
compared  to  the  parameter  estimation  model.  The  most  likely  explanation  for  this 
phenomenon  is  that  the  CMARC  solution  does  not  model  the  effects  of  prop  wash.  The 
higher  dynamic  pressure  over  the  horizontal  stabilizer  from  prop  wash  increases  elevator 
effectiveness.  The  parameter  estimation  model  captures  this  phenomenon.  This  same 
trend  can  be  seen  in  the  modeling  of  rudder  control-power.  An  attempt  could  be  made  to 
model  the  prop  disk  effects  using  CMARC's  normal  panel  velocity  definition  capability. 


METHOD 

CONFIGURATION1 

ELEVATOR  CONTROL  POWER 

(per  raa) 

Cd5€ 

(per  raa) 

Cm5e 

(per  raa) 

Panel  Code 

Blended  Wing-Fuselage/Horiz  Tail 

0.438 

0.01 

-1.199 

Classical2 

Wing  /  Horizontal  Tail 

0.39 

n/a 

-1.04 

Par.  Est3 

Flying  Aircraft 

1.13 

n/a 

-1.62 

C-1724 

0.43 

0.06 

-1.28 

NOTES:  1)  CGX=34.5%  M.A.C.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Classical  design  after  Ref.  [12]. 

3)  Parameter  estimation  from  NPS  flight  test  data  by  Engdahl. 

4)  C-172  data  from  Ref.  [12].. 

Table  4.8  Comparison  of  FROG  UAV  Elevator  Control-Power  Derivatives. 
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2.    Lateral  Directional  Stability  Derivatives 

a.    Static  Lateral-Directional  Stability  Derivatives 

Development  of  the  static  lateral-directional  stability  derivatives  is  not  as 
straightforward.  Both  sides  of  the  airframe  must  be  modeled  by  setting  both  RSYM=1.0 
and  IPATSYM=1.  This  creates  symmetric  patches  around  the  y=0  plane  allowing 
CMARC  to  perform  asymmetric  calculations  around  the  entire  body  and  significantly 
increases  processing  times. 

Initially,  a  lateral-directional  solution  was  attempted  with  a  complete 
airframe  model.  The  longitudinal  model  was  modified  by  adding  the  vertical  stabilizer. 
This  required  what  was  perceived  to  be  a  minor  wake  modification.  Figure  4.14 
illustrates  the  configuration.  The  wing  wake  was  terminated  at  the  wing  root  to  prevent 
interference  with  the  vertical  stabilizer.  This  left  a  wake  gap  for  the  vertical  stabilizer. 
The  engine  nacelle  and  pylon  patches  were  turned  off  because  the  wakes  would  also 
impact  the  vertical  stabilizer.  This  configuration  met  with  mixed  success.  The  values  of 
both  Cyb  and  db  were  opposite  the  expected  directions.  Under  closer  inspection,  it  was 
determined  that  the  false  vortex  shed  from  the  wing  root  caused  a  destabilizing  influence 
on  the  vertical  stabilizer.  Figure  4.15  illustrates  the  vortex  using  CMARC's  built  in  off- 
body  streamline  capability. 

In  order  to  obtain  satisfactory  results,  the  model  is  broken  up  into  two 
separate  groups.  The  aircraft  is  modeled  with  the  blended  wing  and  fuselage  as  one 
group  and  the  horizontal  and  vertical  stabilizers  as  another  group.  Separate  solutions  are 
obtained  for  each  group  and  summed  through  superposition  to  obtain  the  whole  aircraft 
solution.  The  two  separate  models  are  shown  in  Figure  4.16.  This  method  should  be 
used  as  a  last  resort.  A  complete  solution  should  be  used  if  the  airframe  geometry  allows 
adequate  wake  separation  from  the  vertical  stabilizer.  A  complete  solution  would  capture 
fuselage,  wing  and  vertical  stabilizer  interactions. 

Once  the  model  is  defined,  it  is  checked  for  lateral-directional  balance  at 
zero  sideslip  angle  (yaw  angle=0°).  The  side  force,  rolling  and  yawing  coefficients 
should  be  zero  when  a  trial  run  is  performed  at  zero  sideslip.  If  lateral-directional  forces 
or  moments  are  present,  the  model  and  wake  geometry  should  be  checked  for  symmetry. 
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Figure  4.14  Unsuccessful  Lateral-Directional  Model  of  the  FROG  UAV.  False 
Wing  Root  Vortex  Caused  Destabilizing  Influence  on  the  Vertical  Stabilizer. 


Figure  4.15  Destabilizing  Wing  Root  Vortex  Visualized  with  Off-body  Streamlines. 
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Figure  4.16  Final  Simplified  Lateral-Directional  Models  of  the  FROG  UAV. 
Separate  Solutions  are  Summed  to  Obtain  Complete  Airframe  Solution. 

Next,  a  single  CMARC  run  is  performed  witha=atrim  and  with  two 
degrees  of  yaw  angle.  The  results  from  both  models  are  summed.  The  lateral-directional 
derivatives,  Cyp,  Qp  and  Cnp,  are  then  obtained  directly  with  equations  4.22  through  4.24: 


Cv.  =  — —  * per  radian 


A/?c 


71 


C]n  = 


!P      AP 

c, 


1  .180 

-  per  radian 


C„„  = 


71 

180 


per  radian 


4.22 


4.23 


4.24 


nP      AP°      n 

It  should  be  noted  that  the  CMARC  wind  axis  is  modeled  with  x-aft  and  z- 
up,  vice  x-forward  and  z-down  for  the  flight  dynamics  stability  axis.  Care  must  be  taken 
to  reverse  the  signs  of  the  appropriate  coefficients  to  convert  from  the  CMARC  wind  axis 
to  the  flight  dynamics  stability  axis  system. 
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Static  lateral-directional  stability  derivatives  obtained  from  CMARC  are 
presented  in  Table  4.9.  For  comparison  three  other  sets  of  data  are  also  presented.  The 
first  comes  from  classical  analysis  using  methods  from  Roskam  [Ref.  12].  The  second 
set  comes  from  estimates  based  on  data  recorded  from  flight-test  sideslip  maneuvers, 
published  by  Papageorgio  in  Ref.  [2].  The  third  set  comes  from  parameter  estimation,  by 
Engdahl,  based  on  dynamic  flight-test  maneuvers. 


METHOD 

CONFIGURATION1 

STATIC  LAT-DIR  PARAMETERS 

Cyp 
(per  raa) 

Cip 
(per  raa) 

Cnp 
(per  raa) 

CMARC 
Panel  Code 

Wing/Fuselage  +  HorzA/ert  Stabs 

-0.249 

-0.063 

0.063 

Classical 
Design2 

Wing/FuselageA/ert  Tail 

-0.511 

-0.055 

0.051 

Flight  Test3 

Flying  Aircraft 

-0.700 

-0.053 

0.057 

Parameter 
Estimation4 

Flying  Aircraft 

-0.987 

-0.094 

0.176 

C-1725 

-0.310 

-0.089 

0.065 

NOTES:  1)  CGX=34.5%  M.A.C.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Classical  design  calculations  by  Roskam's  methods,  after  Ref.  [12]. 

3)  Flight  test  results  from  steady  heading  sideslip,  from  Ref.  [2] 

4)  Parameter  estimation  from  flight  test  data  by  Engdahl. 

5)  C-172  data  from  Ref.  [12]. 

Table  4.9  Comparison  of  FROG  UAV  Static  Lateral-Directional  Stability 
Derivatives. 

CMARC  produced  weak  results  for  the  Cyp  side  force  derivative.  A  large 
component  of  Cyp,  approximately  60%,  comes  from  the  fuselage.  Side  forces  are  not 
modeled  well  by  the  potential  flow  solution.  In  addition,  the  engine  pylon  and  pod  are 
left  off  the  model  due  to  problems  with  their  wakes  impacting  the  vertical  tail.  They 
most  likely  provide  a  significant  contribution  to  side  force. 
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CMARC  results  for  weathercock  stability,  Cnp,  and  dihedral  effect,  C13, 
show  close  correlation  to  the  classical  calculations.  However,  dynamic  parameter 
analysis  indicates  a  considerably  higher  value  for  both  derivatives.  The  parameter 
estimation  most  likely  captures  the  additional  dynamic  pressure  over  the  vertical 
stabilizer  from  prop  wash.  As  with  the  longitudinal  solutions,  an  attempt  should  be  made 
to  model  the  propeller  disk  with  a  normal  velocity  distribution. 

b.  Lateral-Directional  Damping  Derivatives 

As  with  the  static  case,  the  dynamic  solutions  are  run  with  separate 
wing/fuselage  and  vertical/horizontal  stabilizer  models.  The  solutions  are  summed 
through  superposition.  Figure  4.17  shows  the  two  angular  motions  selected  to  develop 
the  lateral-directional  rate-damping  derivatives.  The  (3-dot  terms  are  generally 
considered  to  be  small.  As  a  result,  the  sideways  plunging  motion  test  case  is  not  run. 
Yaw  rate-damping  is  obtained  directly  from  the  oscillating  yawing  motion  without  the 
requirement  to  subtract  P-dot  effects.  Yawing  motion  is  conducted  at  a  frequency  of  2ii 
rad/s,  which  equates  to  a  reduced  frequency  of  k=0.369  for  this  configuration  and  trim 
airspeed. 


Roll  Damping  Motion 


Yaw  Damping  Motion 


Figure  4.17  Aircraft  Motion  for  the  Determination  of  Roll  and  Yaw  damping. 
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The  roll  axis  is  de-coupled  from  both  angle-of-attack  and  sideslip. 
Therefore,  roll  damping  can  be  obtained  directly  from  a  constant  rolling  motion,  as 
illustrated  in  Figure  4.17.  For  the  FROG  study  a  roll  rate  of  20°/sec  is  selected.  Initially 
a  lower  roll  rate  of  5°/sec  was  used.  However,  the  lower  rate  produced  small  values  of 
both  side  force  and  yawing  moment.  A  higher  roll  rate  prevented  machine  resolution 
from  becoming  a  factor  in  data  reduction.  The  User  Guide  in  Appendix  A  provides  a 
more  detailed  description  of  the  techniques  used  to  gather  the  dynamic  lateral-directional 
derivatives. 

The  roll  damping  stability  terms  are  presented  in  Table  4.10.  Traditional 
design  calculations  from  Roskam  [Ref.  12]  were  completed  for  comparison  to  CMARC 
and  flight-test  results.  In  addition,  C-172  data  is  also  presented.  CMARC  roll  damping, 
Op,  is  50%  higher  than  classical  results,  but  right  in  line  with  parameter  estimation. 
CMARC  adverse  yaw  due  to  roll  rate,  Cnp,  is  approximately  one  fourth  that  obtained  from 
either  classical  calculations  or  parameter  estimation.  Of  note,  CMARC  does  capture  side 
force  due  to  roll  rate.  This  term  is  difficult  to  obtain  from  classical  techniques  and  is 
usually  assumed  to  be  negligible. 


METHOD 

CONFIGURATION1 

ROLL  RATE  TERMS 

cyP 

(per  raaj 

(per  raa) 

Cnp 

(per  raa) 

CMARC 
Panel  Code 

Wing/Fuselage  +  Horz/Vert  Stabs 

0.05 

-0.452 

-0.022 

Classical 
Design2 

Wing/Fuselage/Vert  Tail 

0 

-0.300 

-0.072 

Parameter 
Estimation3 

Flying  Aircraft 

0 

-0.448 

-0.108 

C-1724 

-0.037 

-0.47 

-0.03 

NOTES:  1 )  CGX=34.5%  M.A.C.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Classical  design  calculations  by  Roskam's  methods,  after  Ref.  [12]. 

3)  Parameter  estimation  from  flight  test  data  by  Engdahl. 

4)  C-172  data  from  Ref.  [12]. 

Table  4.10  Comparison  of  FROG  UAV  Roll  Rate  Stability  Derivatives. 
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The  yaw  rate-damping  terms  are  presented  in  Table  4.11.  Traditional 
design  calculations  from  Roskam  [Ref.  12]  are  presented  for  comparison  to  CMARC  and 
flight-test  results.  C-172  data  is  also  listed.  CMARC  does  an  excellent  job  of  predicting 
yaw  damping.  Yaw  damping,  Cnr,  directly  matches  the  result  from  parameter  estimation. 
CMARC  underestimates  roll  moment  due  to  yaw  rate,  Qr,  by  25  to  40  percent.  And,  side 
force  due  to  yaw  rate,  Cyr,  is  overestimated  by  a  factor  of  three  compared  to  the 
parameter  estimation  model.  No  explanation  is  readily  available  to  explain  these  errors. 


METHOD 

CONFIGURATION1 

YAW  RATE  TERMS 

cYr 

(per  raa) 

(per  rarj) 

Cnr 

(per  raa) 

CMARC 
Panel  Code 

Wing/Fuselage  +  HorzA/ert  Stabs 

0.337 

0.121 

-0.121 

Classical 
Design2 

Wing/FuselageA/ert  Tail 

0.140 

0.168 

-0.076 

Parameter 
Estimation3 

Flying  Aircraft 

0.110 

0.208 

-0.121 

C-1724 

0.210 

0.096 

-0.099 

NOTES:  1 )  CGX=34.5%  M.A.C.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Classical  design  calculations  by  Roskam's  methods,  after  Ref.  [12]. 

3)  Parameter  estimation  from  flight  test  data  by  Engdahl. 

4)  C-172  data  from  Ref  [12]. 

Table  4.11   Comparison  of  FROG  UAV  Yaw  Rate  Stability  Derivatives. 
c.  Lateral-Directional  Control-Power  Derivatives 


Aileron  and  rudder  control-power  derivatives  are  presented  in  Tables  4.12 
and  4.13.  Classical  design  calculations  are  performed  using  equations  4.25  through  4.30 
from  Roskam  [Ref.  15]: 


'-o 


4.25 
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^Isr  Cy^ 


(Zv  cos  or  -  lv  sin  a) 


4.26 


"Sr  YSr 


(lv  cos  a-  Zv  sin  a) 


4.27 


CY    » small 


C\     =K*C*C, 


4.28 
4.29 


C,    =\aSa\*  — 


C    (C  \ 


4.30 


where  Cysf,  ctg,  k  and  (CWk)  are  empirical  values  from  charts  in  Roskam  [Ref.  15]. 
Additionally,  flight-test  parameter  analysis  and  C-172  data  are  provided  for  comparison. 


METHOD 

CONFIGURATION1 

AILERON  CONTROL  POWER 

Cy5a 

(per  raa) 

C|5a 

(per  raa) 

CnSa 

(per  raa) 

CMARC 
Panel  Code 

Blended  Wing-Fuselage/Horz/Vert  Tails 

-0.021 

0.194 

-0.0121 

Classical 
Design2 

Wing/Fuselage/Vert  Tail 

0 

0.213 

-0.0236 

Parameter 
Estimation3 

Flying  Aircraft 

0 

0.239 

-0.0261 

C-1724 

0 

0.178 

-0.053 

NOTES:  1)  CGX=34.5%  MAC.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Classical  design  calculations  by  Roskam's  methods,  after  Ref.  [12]. 

3)  Parameter  estimation  from  flight  test  data  by  Engdahl. 

4)  C-172  data  from  Ref.  [12]. 

Table  4.12  Comparison  of  FROG  UAV  Aileron  Control-Power  Derivatives. 
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In  general,  CMARC  provides  aileron  control-power  estimates  in  the 
correct  direction  and  same  order  of  magnitude  as  both  the  classical  and  parameter 
estimation  techniques.  Aileron  roll  control-power  shows  excellent  correlation  to  both 
methods.  However,  CMARC  under-estimates  adverse  yaw  due  to  aileron  deflection, 
Cn6a,  by  approximately  50%.  Yawing  moment  due  to  roll  rate,  Cnp,  is  also 
underestimated  by  a  similar  margin.  An  investigation  into  the  source  of  CMARC 
inaccuracies  in  the  modeling  of  adverse  yaw  due  to  aileron  deflection  and  roll  rate  is 
warranted. 

The  rudder  control-power  derivatives  obtained  with  CMARC  show  good 
correlation  to  classical  design  techniques.  Side  force  due  to  rudder  deflection,  Cy8r,  is 
15%  greater  than  classical  calculations  but  closely  matches  flight-test  results.  For  both 
Cn8r  and  Ci6r,  CMARC  estimates  fall  between  classical  calculations  and  parameter 
estimation  results.  The  stronger  rudder  control-power  observed  in  the  parameter 
estimation  model  is  most  likely  due  to  the  capturing  of  prop  wash  effects. 


METHOD 

CONFIGURATION1 

RUDDER  CONTROL  POWER 

^YSr 

(per  raa) 

C|Sr 

(per  raa) 

Cn6r 

(per  raa) 

CMARC 
Panel  Code 

Blended  Wing-Fuselage/Horz/Vert  Tails 

0.0928 

0.0040 

-0.0453 

Classical 
Design2 

Wing/Fuselage/Vert  Tail 

0.081 

0.0056 

-0.0341 

Parameter 
Estimation3 

Flying  Aircraft 

0.093 

0.0004 

-0.0785 

C-1724 

0.187 

0.015 

-0.0657 

NOTES:  1 )  CGX=34.5%  M.A.C.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Classical  design  calculations  by  Roskam's  methods,  after  Ref.  [12]. 

3)  Parameter  estimation  from  flight  test  data  by  Engdahl. 

4)  C-172  data  from  Ref  [12]. 

Table  4.13  Comparison  of  FROG  UAV  Rudder  Control-Power  Derivatives. 
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3.  Summary  of  CMARC  Stability  Derivative  Analysis 

Table  4.14  lists  the  complete  FROG  aerodynamic  model  obtained  from  CMARC 
for  the  trim  condition.  In  summary,  CMARC  produces  reasonably  accurate  stability 
derivatives  for  an  initial  aerodynamic  model.  The  greatest  difficulty  is  encountered 
modeling  side  force  due  to  sideslip  and  roll  rate.  The  potential  flow  solution  from 
CMARC  fails  to  adequately  capture  the  side  force  on  the  fuselage.  Improved  fidelity 
might  be  obtained  by  modeling  fuselage  flow  separation  with  wake  separation  lines. 
CMARC  also  underestimates  yawing  moment  due  to  roll  rate  and  aileron  deflection  by 
large  margins. 


LONGITUDINAL 

LATERAL-DIRECTIONAL 

Derivative 

Value 

Derivative 

Value 

cL 

0.4295 

Gyp 

-0.2493 

cD 

0.065 

Qp 

-0.0630 

cLa 

4.845 

Cnp 

0.6300 

Coa 

0.2664 

Cyp 

0.0488 

Cmo 

-0.4126 

c,P 

-0.4514 

C|_adOt 

1.420 

Cpp 

-0.0220 

^MadOt 

-6.264 

CYr 

0.3370 

Clc, 

6.862 

c,r 

0.1210 

cDq 

0 

Cnr 

-0.1210 

C^q 

-11.78 

Cy5r 

0.0928 

Cl_6e 

0.4378 

C|gr 

0.0040 

Co8e 

0.0092 

CnSr 

-0.0453 

^MSe 

-1.199 

Cyga 

-0.0206 

C|sa 

0.1943 

Cnga 

-0.0121 

Table  4.14  Summary  of  CMARC  Stability  Derivatives  for  the  NPS  FROG  UAV. 
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F.         COMPARISON  OF  DYNAMIC  AERODYNAMIC  MODELS 

This  section  will  discuss  the  dynamic  response  of  the  FROG  UAV  using  the 
CMARC  generated  aerodynamic  model.  Specifically,  modal  frequency  and  damping  are 
obtained  through  eigenvalue  analysis.  Elevator,  aileron  and  rudder  control  response  is 
found  using  a  linearized  state-equation  model.  Results  are  compared  to  the  classical 
design  and  parameter  estimation  models. 

Dimensional  stability  derivatives  are  placed  into  linearized  4x4  models  as 
outlined  by  Schmidt  in  Ref.  [14].  The  short-period  and  long-period  modes  are  obtained 
using  the  MATLAB  eigenvalue  decomposition  routine.  In  a  similar  manner,  the  lateral- 
directional  plant  matrix  yields  the  roll,  spiral  and  Dutch-roll  modes.  Dynamic  response 
due  to  step  and  doublet  control  inputs  is  obtained  by  assembling  the  complete  linear 
system.  This  paper  uses  control  deflection  sign  convention  consistent  with  Figure  4.9 
from  Ref.  [14].  The  MATLAB  "lsim"  command  produces  the  time  based  output  of  the 
linear  system.  MATLAB  is  used  to  plot  the  dynamic  response. 


1.    Longitudinal  Dynamics 

a.  Longitudinal  Dynamic  Modes 

The  longitudinal  response  of  an  aircraft  can  be  reduced  to  a  series  of  four, 
first-order  differential  equations.  They  are  typically  written  in  matrix  form  based  on  the 
state  variables  u/V,  a,  q  and  0.  The  linearized  state-equations,  in  matrix  form  as 
developed  by  Schmidt  in  Ref.  [14],  are  listed  in  equations  4.31  to  4.37: 

{*}-MM+{B)5.  4'31 


{x}=[u/V     a    q    $}T  4.32 


M-Wkl  4-33 
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{B}=[ln}-] 

M 

vxu 

*a 

0 

-gcos©0 

k]= 

vzu 

VMU 

(v  +  zq) 

-gsin0o 
0 

0 

0 

1 

0 

4.34 


4.35 


&}  =  [**    ^    ^&    o]7 


4.36 


[/.]  = 


V  0  0  0 

0  {v-zd)  0  0 

0  -Ma  1  0 

0  0  0  1 


=    Inertial    matrix 


4.37 


The  non-dimensional  stability  derivatives  generated  by  CMARC  are 
converted  to  dimensional  derivatives  through  the  transformations  listed  in  Table  4.15. 
Dimensional  derivatives  are  more  convenient  because  they  will  lead  to  time  histories 
being  expressed  in  seconds  and  frequency  in  rad/s.  The  MATLAB  script  "froguav.m"  in 
Appendix  D  is  used  to  automate  these  calculations. 

The  linearized,  longitudinal,  state-equation  is  given  by  equation  4.31. 
Note  that  the  plant  and  control  matrices  are  obtained  by  pre-multiplying  by  [In]"1.  The 
eigenvalues  of  the  plant  matrix,  [A],  yield  the  longitudinal  dynamic  modes.  The 
eigenvalues  are  obtained  from  the  MATLAB  script  "froguav.m"  listed  in  Appendix  D. 
The  results  are  presented  in  Table  4.16. 
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Term 

Xu 

xa 
xa 
xq 

z* 
zQ 

z, 

Zj 

Mu 

Ma 
M« 

Mq 

M4 


Description 

mV  <9M 

QS,r        dCD, 


m 

m    V2V 

QS  /c_ 

m    V2V 

qs  acD 

m     56 
QS 


dec 


aci 


<9(ac/2V) 

acD 

5(qc/2V) 


QSrr  +  5Cim 
""m"(CD  +  ^T) 

_QSf_cx       dcL 
m    V2V/  d(dc/2V) 

_QS(_c\       dCL 
m    V2V7  d(qc/2V) 
QS  dCL 

m  5<5 
QScdCm 
IyV  dM 
QSc5Cm 

<9ct 

c 


QSC 


dCr 


Iy      \2VJ  d(ac/2V) 
<9Cm 


QSC 

Iy      V2V/  <>(qc/2V) 
Q  S  c  dC 


Iv      5<5 


Units 
s"1 
ft-s"2 
ft-s"1 
ft-s"1 
ft-s"2 
s"1 
ft-s"2 


ft-: 


.-1 


ft-s"1 
ft-s"2 
ft-s"1 


.-1 


.-1 


Table  4.15  Dimensional  Longitudinal  Stability  Derivatives  from  Ref.  [14]. 


63 


METHOD 

CONFIGURATION1 

LONGITUDINAL  DYNAMIC  MODES 

Short- 

Period 

Phugoid 

en 

(raais) 

c 

(raa/s) 

C 

CMARC 
Panel  Code 

Wing/Fuselage/Horiz  Tail 

5.15 

0.909 

0.358 

0.136 

Classical 
Design2 

Wing/Horiz  Tail  -  5s/5a=0.40 

5.90 

0.734 

0.407 

0.110 

Parameter 
Estimation3 

Flying  Aircraft 

4.67 

0.770 

0.397 

0.101 

C-172 

6.027 

0.685 

0.181 

0.116 

NOTES:  1)  CGX=34.5%  M.A.C.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Classical  design  calculations,  after  Ref.  [12]. 

3)  Classical  model  from  Papageorgio  [Ref.  2]  as  modified  through 
parameter  estimation  from  flight  test  data  by  Engdahl. 

4)  C-172  data  from  Ref.  [12]. 

Table  4.16      Comparison  of  FROG  UAV  Dynamic  Longitudinal  Modes. 

Table  4.16  summarizes  the  longitudinal  dynamic  modes.  The  CMARC 
aerodynamic  model  provides  a  better  estimate  of  short-period  frequency  than  the  classical 
design  technique.  The  CMARC  short-period  frequency  is  10%  faster  than  that  observed 
from  parameter  analysis  and  damping  is  21%  higher.  The  high  damping  in  all  three 
models  will  produce  nearly  deadbeat  results.  The  CMARC  phugoid,  or  long  period,  is 
10%  slower  than  observed  from  flight-test  with  40%  more  damping.  Although,  it  must 
be  pointed  out  that  it  is  difficult  to  accurately  capture  the  phugoid  mode  in  flight-test. 
This  is  due  to  difficulty  the  external  pilot  has  in  accurately  maintaining  wings  level  trim 
during  the  extended  period  required  for  capturing  phugoid  data.  Overall,  the  CMARC 
aerodynamic  model  produces  satisfactory  short-period  and  phugoid  modal  data  for  the 
development  of  closed-loop  flight  controls.  Of  note,  CMARC  provides  a  more  accurate 
prediction  of  the  important  short-period  natural  frequency.  The  larger  error  in  phugoid 
frequency  is  less  important  to  the  design  of  a  closed  loop  controller. 
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b.  Longitudinal  Dynamic  Response  to  Control  Input 

Following  dynamic  mode  analysis,  the  response  to  an  elevator  step-input 
and  a  doublet-input  is  modeled  using  MATLAB.  Dynamic  response  for  the  CMARC, 
classical  and  parameter  analysis  models  are  overlaid  for  comparison.  The  linear  system 
is  set  up  using  the  MATLAB  script  "froguav.m"  in  Appendix  D.  The  "lsim"  command 
outputs  the  time-based  response  from  a  time-based  control  input  vector.  A  -2°  (TEU) 
elevator  step-input  is  selected  to  keep  the  response  in  the  linear  region.  A  5°  elevator 
doublet  with  a  1 .2  second  period  is  used  to  excite  the  short-period  mode.  The  selected 
doublet  period  closely  matches  the  short-period  mode.  Dynamic  response  is  displayed  in 
Figures  4.18  and  4.19. 

The  FROG  UAV  response  to  a  -2°  (TEU)  step  elevator  input  is  displayed 
in  Figure  4.18.  All  three  aerodynamic  models  are  presented  for  comparison.  The  model 
adjusted  by  parameter  analysis  shows  a  much  larger  response  than  either  the  CMARC  or 
classical  design  models.  The  CMARC  model  produces  a  final  pitch  angle  change  after  4 
seconds  of  35°  versus  45°  for  the  flight-tested  model.  On  closer  review,  it  is  noted  the 
Cm8e  =-1.199  from  CMARC  is  considerably  less  than  the  Cm5e  =  -1.621  from  flight-test. 
It  should  be  noted  that  the  parameter  estimation  model  compensates  for  real  world  factors 
including;  air  load  distortions,  sensor  measurement  errors  and  prop  wash.  In  other  words, 
the  CMARC  model  assumes  uniform  displacement  and  perfect  sensors  while  the  flight- 
test  model  is  empirically  fit  to  a  measured  response. 

The  response  to  a  5°  elevator  doublet  is  displayed  in  Figure  4.19.  All 
three  models  show  a  similar  frequency  and  a  high  degree  of  damping.  As  expected,  the 
parameter  estimation  model  produces  a  larger  response.  The  magnitude  of  the  CMARC 
response  is  approximately  two  thirds  of  that  observed  from  the  flight-test  model.  Again, 
this  is  due  to  the  much  larger  value  of  Cm6e  obtained  from  empirically  fitting  the 
parameter  estimation  model  to  observed  aircraft  response.  Modeling  prop  wash  with 
CMARC  should  produce  improved  results. 

In  summary,  the  CMARC  aerodynamic  model  showed  satisfactory 
longitudinal  control  response  in  frequency  and  damping.  It  is  recommended  that  Cm8e  be 
adjusted  to  increase  elevator  response. 
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Figure  4.18  FROG  UAV  Dynamic  Response  to  a  -2°  (TEU)  Elevator  Step  Input 
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LONGITUDINAL  RESPONSE  TO  A  5  DEGREE  ELEVATOR  DOUBLET 
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Figure  4.19  FROG  UAV  Dynamic  Response  to  a  5°  Elevator  Doublet 
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2.  Lateral-Directional  Dynamics 

a.  Lateral-Directional  Dynamic  Modes 

The  lateral-directional  response  of  an  aircraft  can  also  be  reduced  to  a 
series  of  four  first-order  differential  equations.  They  are  typically  written  in  matrix  form 
based  on  the  state-variables  p\  p,  <)>,  and  r.  The  linearized  state-equations  in  matrix  form 
as  developed  by  Schmidt  in  Ref.  [14]  are  listed  in  equations  4.38  to  4.44: 


{x}=[A}{x}+{B}5alr 
{x}=\p     p     <{>     r]7 

M=[/„rk] 


4.38 
4.39 
4.40 

4.41 


k]= 


YB      F       gcos0o     (Yr+V) 


h 

LP 

0 

Lr 

0 

1 

0 

0 

N? 

N, 

0 

A^ 

4.42 


{*.}  = 


8a 


L*     0    N^ 


>' 


Or 


-\T 


K    o    NSr 


4.43 


[/„]  = 


V          0  0    0 

0  (V-Zd)  0    0 

0  -Ma  10 

0           0  0     1 


=    Inertial    matrix 


4.44 
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As  with  the  longitudinal  axis,  the  stability  derivatives  generated  by 
CMARC  are  dimensionless.  However,  it  is  more  convenient  to  use  dimensional 
derivatives.  They  will  lead  to  time  histories  being  expressed  in  seconds  and  frequency  in 
rad/s.  The  lateral-directional  stability  derivatives  are  turned  into  dimensional  derivatives 
through  the  transformations  listed  in  Table  4.17  from  Ref.  [14].  The  MATLAB  script 
"froguav.m"  in  Appendix  D  is  used  to  automate  these  calculations. 


Term 


Yi3 
YP 

Yr 

Ys 

u 

Lp 

Lr 

U 

>"p 
Nr 


Description 


d(3 


QS 

m 

QS 

m 


Q^acy 

m 
Jb_ 
2  V 

_b_ 
2V 


acv 


<9(pb/2V) 

dCy 
d(rb/2V) 

QS     5Cy 

m     dS 
QSb  dCt 

Ix        5/? 


QSb 

Ix 

QSb 

Ix 


<9Q 


2v;  a(pb/2V) 

b  \        <9Q 


QSb 

I, 

QSb 


2V;   d(rb/2V) 
QSb  dCt 

Ix      <?<5 
QSb  dCn 

Iz      5/? 
_b_\        <3Cn 
2VJ  d(pb/2V) 
_b_\        aCn 
2Vy   d(rb/2V) 
QSb  dCn 


I, 


5.5 


Units 


ft-s"2 
ft-s"1 


ft-s"1 


ft-s"2 


.-i 


.-i 


.-i 


,-i 


Table  4.17  Dimensional  Lateral-Directional  Stability  Derivatives  from  Ref.  [14] 
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The  linearized,  lateral-directional  state-equation  is  given  by  equation  4.38. 
Note  that  the  plant  and  control  matrices  are  obtained  by  pre-multiplying  by  [In]"1.  The 
dimensional  derivatives  from  Table  4.16  populate  the  linearized  4x4  plant  matrix  and 
control  matrices.  The  plant  matrix,  [A],  is  used  to  obtain  the  lateral-directional  dynamic 
modes.  The  eigenvalues  are  obtained  from  MATLAB  using  the  "froguav.m"  script  in 
Appendix  D.  The  results  are  presented  in  Table  4.18. 


METHOD 

CONFIGURATION1 

LATERAL-DIRECTIONAL  DYNAMIC  MODES 

Dutc 

h-Roll 

Roll 

Spiral 

(raa/s) 

c 

©n 

(raa/s) 

(raa/s; 

CMARC 
Panel  Code 

Wing/Fuselage  Group 
plus  HorizA/ert  Stab  Group 

2.54 

0.13 

-3.80 

0.000 

Classical 
Design2 

Wing/Fuselage/Horiz  Tail 

2.48 

0.09 

-2.83 

0.065 

Parameter 
Estimation3 

Flying  Aircraft 

4.27 

0.14 

-3.97 

0.090 

C-1724 

3.38 

0.20 

12.43 

0.01 

NOTES:  1)  CGX=34.5%  M.A.C.  /  CGZ=8.6"  from  bottom  of  fuselage. 

2)  Classical  design  calculations  after  Ref.  [12]. 

3)  Classical  model  from  Papageorgio  [Ref.  2]  as  modified  through 
parameter  estimation  from  flight  test  data  by  Engdahl. 

4)  C172  data  from  Ref.  [12]. 

Table  4.18      Comparison  of  FROG  UAV  Lateral-Directional  Dynamic  Modes. 

The  lateral-directional  dynamic  modes  are  summarized  in  Table  4.18.  The 
CMARC  aerodynamic  model  provides  significantly  better  estimates  of  frequency  and 
damping  than  the  classical  design  technique  for  both  the  Dutch-roll  and  roll  modes.  The 
CMARC  Dutch-roll  natural  frequency,  although  better  than  the  classical  response,  is  40% 
lower  than  predicted  by  parameter  estimation.  Adjustments  should  be  made  to  the 
CMARC  aerodynamic  model  to  improve  Dutch-roll  response. 
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The  primary  contributor  to  the  Dutch-roll  frequency  is  weathercock 
stability,  Cnp.  A  review  of  Table  4.9  shows  that  Cnp  obtained  from  CMARC  is 
considerably  less  than  the  value  obtained  through  parameter  estimation.  One  potential 
source  of  error  is  the  lack  of  modeling  the  increased  dynamic  pressure  due  to  prop  wash. 
The  FROG  model  should  be  re-worked  to  include  a  prop  disk  to  investigate  the  ability  of 
CMARC  to  capture  prop  wash  effects. 

The  CMARC  model  predicts  a  neutral  spiral  mode.  This  is  most  likely 
due  to  a  combination  of  equal  and  opposite  Qr  and  Cnr  ratios.  When  the  Cir/Cnr  ratio  is 
changed  to  -2  (approximately  that  of  the  analytical  and  parameter  estimation  models),  the 
spiral  mode  comes  out  close  to  the  analytical  model.  The  weak  roll  due  to  yaw-rate,  Qr, 
in  the  CMARC  solution  seems  to  be  the  main  source  of  error  in  the  spiral  mode.  The 
value  will  need  to  be  modified  to  produce  an  acceptable  spiral  response. 

b.  Lateral-Directional  Dynamic  Response  to  Control  Input 

Following  dynamic-mode  analysis,  the  response  to  aileron  and  rudder 
inputs  is  modeled  using  MATLAB.  Dynamic  response  for  the  CMARC,  classical  and 
parameter-analysis  models  are  overlaid  for  comparison.  The  linear  system  is  set  up  using 
the  MATLAB  script  "froguav.m"  in  Appendix  D.  The  "lsim"  command  outputs  the  time- 
based  response  from  a  time-based  input  control  vector. 

The  lateral  response  to  a  +2°  (right  wing  down)  step  aileron  deflection  is 
shown  in  Figure  4.20.  CMARC  steady-state  roll-rate  is  less  than  that  predicted  by  the 
other  two  models.  This  is  expected  considering  the  CMARC  model  has  the  lowest  ratio 
of  aileron  control-power  to  roll-damping.  In  addition,  the  side- force  terms,  Cyp  and  Cy8a, 
modeled  by  CMARC  both  fight  roll-rate.  They  are  not  included  in  the  other  models. 
Perhaps  they  should  be  assumed  small  and  set  to  zero.  CMARC  does  a  significantly 
better  job  of  modeling  the  Dutch-roll  response.  The  roll-rate  and  sideslip  traces  from  the 
classical  design  model  clearly  show  excess  excitation  of  the  Dutch-roll  mode.  The 
CMARC  model  shows  a  better  correlation  to  the  parameter  estimation  model  for  Dutch- 
roll  amplitude  and  damping.  As  predicted  by  the  eigenvalues,  Dutch-roll  frequency  is 
slower  than  the  parameter  estimation  model. 

The  lateral-directional  response  to  a  5°  aileron  doublet  is  shown  in  Figure 
4.21.  The  1.5  second  period  clearly  excites  the  Dutch-roll  mode  in  all  three  models.  As 
with  the  step-input  response,  the  CMARC  model  provides  a  significantly  better  match  to 
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Dutch-roll  excitation  than  the  classical  design  model.  Amplitude  and  damping  show  a 
close  correlation  to  the  parameter  estimation  model.  A  lower  Dutch-roll  frequency  is 
evident  in  the  CMARC  sideslip  trace. 

The  lateral-directional  response  to  a  +2°  (nose  left)  rudder  deflection  is 
shown  in  Figure  4.22.  The  CMARC  model  demonstrates  a  close  match  to  the  classical 
design  model  for  Dutch-roll  frequency  and  damping.  CMARC  shows  higher  sideslip 
excitation  than  the  classical  design  model.  A  slower  Dutch-roll  frequency  is  evident  in 
both  the  CMARC  and  classical  design  models. 

The  lateral-directional  response  to  a  5°  rudder  doublet  deflection  is  shown 
in  Figure  4.23.  Again,  the  1.5  second  doublet  excites  the  Dutch-roll  mode  in  all  three 
models.  CMARC  provides  a  similar  response  to  the  classical  design  model.  As 
expected,  Dutch-roll  natural  frequency  is  about  40%  slower  than  the  parameter  estimation 
model  with  similar  damping. 

In  summary,  the  CMARC  lateral-directional  model  provides  FROG 
dynamics  that  are  similar  to  the  classical  design  calculations.  However,  current  CMARC 
lateral-directional  dynamics  are  not  adequate  for  closed-loop  controller  or  autopilot 
design.  Minor  adjustments  will  be  required  to  better  match  observed  flight 
characteristics.  The  lateral-directional  model  should  be  modified  to  include  a  propeller 
disk.  Higher  dynamic  pressure  over  the  tail  surfaces  should  provide  a  stronger  Cnp 
derivative  and  result  in  faster  Dutch-roll  frequency.  In  addition,  the  FROG  model  should 
be  adjusted  to  provide  a  higher  aileron  control-power  to  roll-damping  ratio.  This 
modification  will  improve  roll-rate  response. 
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Figure  4.20  FROG  UAV  Dynamic  Response  to  a  +2°  (RWD)  Aileron  Step  Input 
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Figure  4.21  FROG  UAV  Dynamic  Response  to  a  5°  Aileron  Doublet 
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LATERAL-DIRECTIONAL  RESPONSE  TO  A  2  DEGREE  RUDDER  STEP  INPUT 
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Figure  4.22  FROG  UAV  Dynamic  Response  to  a  +2°  (TEL)  Rudder  Step  Input 
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LATERAL-DIRECTIONAL  RESPONSE  TO  A  5  DEGREE  RUDDER  DOUBLET 
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Figure  4.23  FROG  UAV  Dynamic  Response  to  a  5°  Rudder  Doublet 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

CMARC  is  a  DOS  personal-computer  hosted  panel-code  adopted  from  the  NASA 
Ames  PMARC  code.  AeroLogic,  Inc.,  created  CMARC  by  converting  PMARC 
FORTRAN  77  source  code  into  the  C  language.  Significant  memory  management  and 
command-line  enhancements  were  also  added.  CMARC  solves  for  inviscid, 
incompressible  flow  over  complex  three-dimensional  bodies.  Emphasis  in  this  study  is 
placed  on  expanding  the  use  of  CMARC  to  develop  a  full  aerodynamic  model  of  the 
Naval  Postgraduate  School  FROG  UAV.  The  CMARC  stability  derivatives  are 
compared  to  models  derived  from  classical  design  calculations  and  parameter  estimation. 
In  addition,  pitot-static  and  angle-of-attack  sensor  position  corrections  are  obtained 
through  CMARC  analysis. 

The  LOFTSMAN  and  POSTMARC  portions  of  the  Personal  Simulation  Works 
software  suite  are  used  exclusively  for  the  pre-process  modeling  and  post-process 
visualization  of  CMARC  files.  The  LOFTSMAN  capability  to  automatically  format  and 
generate  CMARC  input  patches  is  an  enhancing  characteristic.  Functionality  should  be 
added  to  allow  the  modeling  of  wing  tip  ribs  that  are  not  parallel  to  the  aircraft  centerline. 

POSTMARC  is  an  excellent  tool  for  visualizing  CMARC  output  files.  The 
capability  to  create  streamlines  and  perform  boundary-layer  calculations  external  to 
CMARC  is  extremely  useful.  However,  much  time  could  be  saved  if  POSTMARC 
maintained  previous  settings  and  selections  following  translations,  rotations  and  re- 
scaling.  Additionally,  a  capability  to  overlay  multiple  data  types  is  desired. 

CMARC  off-body  flow-field  analysis  is  useful  for  both  static  source  and  angle-of- 
attack  sensor  position  corrections.  In-flight  measurements  may  be  corrected  using  look- 
up tables  or  through  curve-fits  of  CMARC-derived  data.  Flight  testing  is  recommended 
for  validation  of  sensor  corrections  obtained  from  the  CMARC  off-body  analysis. 

For  the  static  longitudinal  analysis,  CMARC  produces  accurate  values  for  a^ 
and  CMa  and  a  slightly  high  value  of  CLa.  Elevator  control-power  (CM6e)  from  CMARC  is 
considerably  weaker  than  the  value  obtained  from  parameter  estimation.  CMARC- 
derived  pitch-damping  is  stronger  than  the  pitch-damping  from  both  classical  design 
calculations  and  parameter  estimation. 

Longitudinal  dynamic-mode  analysis  shows  an  acceptable  match  for  short-period 
and  long-period  frequency  and  damping.    As  expected,  a  lower  dynamic  response  to 
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elevator  control  inputs  is  observed  when  compared  to  the  parameter  estimation  model. 
However,  in  all  cases,  the  CMARC  aerodynamic  model  demonstrates  better  dynamic 
response  than  the  classical  design  model.  A  propeller  disk  should  be  added  to  the 
CMARC  model  in  an  attempt  to  capture  prop-wash  effects  on  the  horizontal  stabilizer 
and  elevator.  Increased  elevator  control-power  would  provide  a  better  match  with 
observed  flight  characteristics. 

In  general,  CMARC  provides  satisfactory  lateral-directional  derivatives.  Roll- 
damping  closely  matches  the  value  obtained  through  parameter  estimation.  The  greatest 
difficulty  encountered  is  in  modeling  side-force  due  to  both  sideslip  and  roll-rate.  The 
potential  flow  solution  from  CMARC  fails  to  adequately  capture  the  side  force  on  a 
slender-body  fuselage.  Additional  wakes,  placed  along  separation  lines,  may  improve 
fuselage  side-force  prediction. 

Dynamic  response  to  rudder  control  input  shows  close  correlation  in  amplitude 
and  damping  to  the  parameter  estimation  model.  However,  both  the  CMARC  and 
classical  models  of  the  FROG  UAV  produced  a  slower  Dutch  roll  frequency.  Steady 
state  roll-rate  obtained  from  the  CMARC  model  is  somewhat  slower  than  either  the 
classical  design  or  parameter  estimation  models. 

Overall,  the  CMARC  panel  code  is  found  to  be  suitable  for  aerodynamic 
modeling  of  the  Naval  Postgraduate  School  FROG  UAV.  CMARC-derived  stability 
derivatives  are  sufficiently  accurate  for  incorporation  into  an  initial  aerodynamic  model. 
The  CMARC  aerodynamic  model  demonstrated  better  longitudinal  dynamic  response 
than  the  classical  design  model.  Lateral-directional  response  is  similar  to  that  obtained 
from  classical  design  techniques.  Adjustment  through  comparison  with  flight-test  data  is 
still  required  to  optimize  the  CMARC  model.  Future  studies  should  concentrate  on 
improving  CMARC  modeling  of  fuselage  side  force  through  the  addition  of  separation 
wake  lines.  Additionally,  the  propeller  disk  should  be  modeled  in  an  attempt  to  capture 
the  effects  of  increased  dynamic  pressure  over  the  horizontal  and  vertical  tail  surfaces. 
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A.         INTRODUCTION 

1.  CM  ARC  Panel  Code  and  Limitations 

CMARC  is  a  PC  hosted  panel  code.  The  code  is  useful  for  gathering  aircraft 
stability  data  subject  to  the  limitation  of  inviscid,  potential  flow  solutions.  Panel  codes 
assume  fully  attached  flow.  Therefore,  care  must  be  taken  to  gather  stability  data  in  the 
linear  region  of  the  lift  curve  slope  where  separation  effects  are  minimized.  The  results 
obtained  from  CMARC  should  always  be  spot  checked  with  classical  design  calculations. 

CMARC  and  PMARC  have  been  shown  to  produce  functionally  equivalent 
results.  The  guide  describes  the  use  of  DOS  batch  files  for  automating  CMARC  data 
runs  on  a  PC.  PMARC,  hosted  on  an  SGI  machine,  will  require  the  use  of  UNIX  batch 
commands. 

2.  Overview  of  Obtaining  Stability  Derivatives  with  CMARC 

An  overview  of  the  methods  to  obtain  stability  derivatives  follows.  The  detailed 
description  and  data  reduction  methods  can  be  found  in  the  appropriate  section  of  the  user 
guide. 

For  the  static  stability  derivatives,  the  CMARC  model  is  run  at  two  different 
angles-of-attack  and  one  sideslip  angle.  The  static  solutions  are  obtained  without  control 
surface  deflections.  The  slopes  of  the  force  and  moment  coefficients  are  then  taken  to 
produce  the  CLa  and  Cma  longitudinal  derivatives  and  the  CYp,  Cip  and  CnP  lateral- 
directional  derivatives. 

For  the  control  power  derivatives,  the  model  is  run  at  the  trim  condition  with 
successive  elevator,  aileron  and  rudder  control  deflections.  The  difference  between  the 
results  with  and  without  control  surface  deflections  yield  the  CL5e,  CM5e  and  CD5e 
longitudinal  and  the  CY5r,  C16r>  Cn5r,  C15a  and  Cn5a  lateral-directional  control  power 
derivatives. 

Development  of  the  damping  derivatives  is  not  as  straight  forward.  For  the 
dynamic  derivatives,  motion  is  enabled  around  the  center  of  gravity.  With  motion  turned 
on,  it  is  extremely  important  to  ensure  that  the  global  origin  and  C.G.  are  co-located.  If 
they  are  not,  a  rotation  around  a  specific  axis  will  also  create  a  translation  at  the  C.G. 
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For  the  longitudinal  damping  derivatives,  the  model  is  run  with  oscillating  vertical 
plunging  motion  to  obtain  the  CL  and  CM  a-dot  terms.  The  lift  and  pitching  moment 
coefficients  are  broken  into  real  (in  phase  with  AOA)  and  imaginary  (out  of  phase  with 
AOA)  components.  The  imaginary  components  are  due  to  a-dot  effects.  Next,  the  model 
is  run  with  oscillating  pitch  motion  to  obtain  the  combined  a-dot  and  pitch  rate  terms. 
Subtracting  the  a-dot  influence  obtained  from  the  plunging  motion  isolates  the  pitch  rate 
damping  term  from  the  pitch  motion. 

For  the  lateral-directional  analysis,  the  P-dot  terms  are  generally  negligible.  This 
allows  the  model  to  be  run  with  just  oscillating  roll  and  yaw  motion.  As  with  the 
longitudinal  test  case,  the  imaginary  or  out  of  phase  component  yields  the  combined  (3- 
dot  and  rate  damping  data.  With  the  p-dot  terms  assumed  negligible,  the  oscillating 
motion  yields  the  CY,  C„  and  Cn  roll  and  yaw  rate  terms  directly. 

3.  Coefficient  Normalization  and  Stability  Axis  System 

CMARC  contains  built-in  functionality  to  integrate  forces  and  moments  over  the 
surface  of  a  body.  Forces  and  moments  are  automatically  normalized  into  non- 
dimensional  coefficients  based  on  the  mean  aerodynamic  chord,  reference  wing  area, 
semi-span  and  center  of  gravity  location  in  the  CMARC  BINP9  input  line.  CMARC 
outputs  coefficients  in  both  the  wind  and  body  axes.  Of  note,  CMARC  uses  the  semi- 
span  to  normalize  rolling  and  yawing  moment  coefficients.  Most  texts  on  stability  and 
control,  including  Roskam's  "Aircraft  Flight  Dynamics"  and  Etkin's  "Dynamics  of 
Flight,"  normalize  rolling  and  yawing  moments  by  span.  This  user  guide  will  normalize 
roll  and  yaw  moments  by  span.  Table  A.l  summarizes  the  factors  for  normalizing 
moments  and  angular  rates.  All  rolling  and  yawing  moment  coefficients  presented  in  this 
guide  have  been  normalized  with  span  by  dividing  the  CMARC  output  by  a  factor  of  two. 
Table  A.l  also  indicates  the  characteristic  time,  t*,  employed  for  angle  rate  data 
reduction. 

In  addition  to  differences  in  normalizing  moments,  CMARC  uses  the  typical  CFD 
axes  system  shown  in  Figure  A.l.  This  user  guide  will  perform  all  calculations  in  the 
stability  axes  system  as  illustrated  in  Figure  A.l.  The  sign  of  CMARC  roll  and  yaw 
moments  need  to  be  reversed.  The  direction  for  positive  control  deflections  is  also  shown 
in  Figure  A.l.  All  control  surfaces  are  patched  with  positive  deflections  using  the 
convention  in  Figure  A.  1 . 
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Note:  1 )  CMARC  normalizes  roll  and  yaw  coefficients  with  b/2. 

Table  A.l  Normalized  Moment  and  Rate  Equations. 
B.         MODEL  VALIDATION 

1.         Initial  Solution  at  the  Trim  Condition 

Once  a  CMARC  model  is  built  up,  a  reference  solution  is  obtained  at  the  desired 
trim  condition.  Initial  hand  calculations  or  flight  test  data  can  be  used  to  select  the 
estimated  cruise  angle-of-attack.  Use  POSTMARC  to  visualize  the  reference  solution. 
Check  to  make  sure  the  pressure  distribution  is  consistent  with  expected  results. 
Stagnation  zones  (Cp=l)  should  be  evident  on  the  leading  and  trailing  edges  of  the  flight 
surfaces  and  aerodynamic  bodies.  Lower  pressure  should  be  observed  on  the  top  of 
lifting  surfaces  and  tip  loss  should  be  observed  near  the  wing  tips.  Wakes  should  trail 
from  the  expected  surfaces. 

If  abnormal  pressure  distribution  is  observed,  check  for  inverted  panel  orientation. 
This  is  accomplished  by  changing  the  color  for  positive  and  negative  panel  orientations 
using  the  "View  »  Color"  pull-down  menu.  The  default  color  is  white  for  both 
orientations.  Change  the  inside  to  a  contrasting  color  such  as  red.  All  panels  should  face 
in  the  positive  out  direction.  If  a  patch  has  panels  oriented  inward,  change  the  orientation 
of  the  patch  with  the  "IREV"  term  in  the  initial  patch  definition. 

Check  that  the  model  is  symmetric  around  the  lateral  and  directional  axes  with 
sideslip  set  to  zero.  Residual  rolling  or  yawing  moment  without  a  control  deflection  is  an 
indication  of  an  asymmetric  model. 


85 


CMARC  Axis  System 
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Vel. 
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+6  =  Pitch  Angle  (8=q) 
+H  =•  Pi  ten  Moment 


+w  =  Vertical 

Ve 
+Z   =■  Normal   r*  ) 


^  1MX  <  0 


Force 


■t-t  =  Yaw  Angle  (tar) 
♦N  =  Yaw  Moment 

z  r 

Stability  Axis  System 


Figure  A.l    CMARC  CFD  Axis  System  Compared  To  Stability  Axis  System.   From 
References  [3]  and  [14]  Respectively. 
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Pitching  moment  should  be  checked  for  a  reasonable  value.  Don't  forget  to  verify 
that  the  origin  of  the  global  coordinate  system  is  set  to  the  center-of-gravity  (C.G.).  It  is 
important  to  remember  that  the  damping  derivatives  will  be  obtained  by  inducing  motion 
around  the  center-of-gravity.  If  the  C.G.  and  origin  are  not  collocated,  angular  motion 
will  generate  undesired  translation  at  the  C.G. 

Next,  the  lift  coefficient  from  the  initial  solution  is  checked  to  make  sure  it  is  in 
the  "ballpark."  If  a  major  discrepancy  is  noted,  the  input  file  will  have  to  be  checked 
closely  for  error.  Be  sure  to  check  the  normalizing  factors  such  as  wing  area,  reference 
chord  and  semi-span.  If  every  thing  looks  good,  try  increasing  or  decreasing  the  angle- 
of-attack  to  see  if  a  slight  adjustment  will  bring  the  lift  coefficient  in  line  with  the  hand 
calculation.  Once  the  model  exhibits  symmetry  with  a  reasonable  lift  distribution,  you're 
ready  to  start  gathering  data. 
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C.         LONGITUDINAL  STABILITY  DERIVATIVES 

This  section  will  describe  CMARC  methods  for  the  development  of  longitudinal 
stability  derivatives  to  include  the  static,  rate  damping  and  control  power  derivatives.  As 
stated  earlier,  the  results  obtained  from  CMARC  should  be  checked  against  classical 
design  calculations.  The  potential  flow  analysis  performed  by  CMARC  does  not  provide 
accurate  viscous  drag  values.  Total  drag  can  be  estimated  from  the  flight  test  power-off 
glide  ratio  or  cruise  thrust  required.  If  this  data  is  unavailable,  then  empirical  data  from 
publications  such  as  Hoemer's  "Fluid  Dynamic  Drag"  can  be  used  to  estimate  total  drag. 

For  the  longitudinal  analysis,  only  half  the  model  is  analyzed.  This  cuts  the 
model  size  in  half,  resulting  in  much  quicker  solutions.  The  symmetric  calculation  mode 
is  selected  by  setting  both  RSYM=0.0  and  IPATSYMO  in  the  CMARC  input  file. 
Remove  the  vertical  tail  patch  if  it  interferes  with  wing  or  fuselage  wakes.  Figure  A.2 
shows  the  FROG  model  configuration  used  to  find  the  longitudinal  derivatives.  A  rigid 
wake  is  selected  that  is  continuous  from  wingtip  to  wingtip.  As  will  be  seen  in  a  later 
section,  the  vertical  tail  is  added  for  the  lateral-directional  solutions.  This  requires  the  use 
of  a  modified  wing  wake  to  avoid  the  vertical  tail. 


7 

!_ 
z      x 


Blended  Wing  and  Fuselage  with  Horizontal  Tail 


Figure  A.2  FROG  UAV  Model  for  Obtaining  Longitudinal  Stability  Derivatives. 
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1.  Static  Longitudinal  Stability  Derivatives 

a.  Static  Longitudinal  Stability  Equations 

Three  basic  longitudinal  stability  derivatives  can  be  measured  with  just 
two  runs  of  the  CMARC  model.  The  model  is  first  analyzed  at  an  angle-of-attack 
corresponding  to  the  estimated  trim  condition.  For  the  FROG  case,  a  =0°  is  selected  for 
the  first  run  based  on  hand  calculations.  A  second  CMARC  run  is  conducted  with  angle- 
of-attack  incremented  by  one  or  two  degrees.  CL,  CD  and  Cm  are  then  extracted  manually 
from  the  data  files.  The  slope  of  CL,  CD  and  Cm  versus  angle-of-attack  provides  the  CLa, 
CDa  and  Cma  longitudinal  derivatives.  Additional  solutions  should  be  obtained  to  check 
for  lift  slope  linearity.  In  addition,  aaim  is  calculated  from  the  lift  curve  slope  and  trim  lift 
coefficient.  Equations  A.  1  through  A.4  are  used  for  these  calculations. 

CL    =—7-^ \-* per  radian  A.l 

a        V*2-a\)       n 


[CDl  -C^j^igp 
>a        {a2-ax)        n 


Cd    =    t    -J-* per  radian  A.2 


Cma=[Cp'Cm^—  Per  radian  A.3 

a°trim  =  a°i  +^5 —*—  degrees  A.4 


b.  Sample  Static  Longitudinal  Stability  Data  Reduction 

Table  A.2  presents  CMARC  solutions  for  the  FROG  UAV  at  two  angles- 
of-attack.  Sample  calculations  for  obtaining  the  static  derivatives  are  demonstrated 
below.  The  calculations  are  easily  implemented  in  a  spreadsheet  or  with  a  MATLAB 
script. 
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RUN# 

AOA 

CL 

CD 

CY 

C  m 

C  n 

C  1 

1 

0° 

0.4300 

0.0277 

0.0000 

0.0444 

0.0000 

0.0000 

2 

2° 

0.5991 

0.0370 

0.0000 

0.0300 

0.0000 

0.0000 

Table  A.2  FROG  UAV  Static  Longitudinal  Stability  Data. 


Sample  Calculations: 


C 


L 


(0.5991-0.4300).  180 


a 


(2-0) 


K 


=  4.8444  per  radian 


A.5 


C, 


_  (0.0370 -0.0277)^  180  _ 


'a 


(2-0) 


K 


=  0.2664  per  radian 


A.6 


r 


(0.0300-0.0444).  180 


m 


a 


(2-0) 


71 


=  -0.4125  per  radian 


A.7 


<*°trim=0°  + 


(0^295-0.4300).180 


4.8444 


71 


A.8 


c.  Sample  Total  Drag  Coefficient  Calculation 

The  final  static  longitudinal  parameter  required  is  total  aircraft  drag.  CD 
plays  an  important  role  in  long  period  aircraft  dynamics.  Unfortunately,  potential  flow 
panel  codes  such  as  CMARC  do  not  provide  accurate  total  drag  estimates.  They  can 
provide  good  induced  drag  predictions.  In  addition,  if  equipped  with  a  boundary  layer 
code  like  that  contained  in  CMARC,  they  can  provide  integrated  skin  friction  results. 
However,  a  large  total  drag  contribution  in  the  form  of  separation  drag  is  not  accounted 
for.  Total  drag  estimates  are  made  below  using  the  two  simple  techniques  shown  in 
Equations  A.9  to  A.  12.  The  first  method  is  based  on  the  flight-tested  glide  ratio.  The 
second  is  based  on  cruise  power  required  and  estimated  prop  efficiency.  Note  that  the 
selected  prop  efficiency  is  relatively  low  due  to  the  FROG  UAV  small  propeller  diameter, 


90 


high  RPM  and  pusher  configuration.   The  two  methods  provide  drag  predictions  within 
10%  of  each  other.  The  results  are  averaged  to  CD=0.065. 

Method  1 :  Lift-to-Drag  Ratio  (L/D=7  from  flight  test) 

L/=1  =>  0  =  £  =  £  =  ^ZJfe=9.67  lbs  A.9 

/D  111 

D                                         9.67    lbs 
CD=  —  = =  0.0611  A>10 

q*      0.5*0.002327    lb  ■  s2 1  ft"  *  882    ft2/s2*\7.57    ft2 


Method  2:  Cruise  Power  Setting  (HP=5,  r)P=0.35) 
D      TR  _HPR* 550   ftlbslslHP*T]pIV 

Lr)  — — — — — ■ —  A.  1 1 

qS      qS  qS 

5    HP*550   ft-lbs/s/HP*035/SS    ftls 
CD  = =  0.069  A.12 

0.5*0.002327   lb  ■  s2 1  ft"  *  882    fi2ls2*\1.51    ft2 

2.         Longitudinal  Control  Power  Stability  Derivatives 

cu         Longitudinal  Control  Power  Equations 

The  elevator  control  power  derivatives  are  obtained  by  substituting  a  0° 
deflection  horizontal  tail  patch  for  one  with  positive  elevator  deflection.  Only  one  run  is 
required.  For  the  FROG  UAV  study,  +5°  (TED)  deflection  is  used.  The  difference 
between  the  trim  condition  and  the  deflected  value  is  divided  by  the  elevator  deflection  as 
shown  below.  Note  that  CD5e  from  CMARC  only  includes  induced  drag  due  to  elevator 
deflection: 

lC45e2   "C45eJ     180  J  A.13 

C/o   = per    raa 

**  Se2  -5e\  n 


lCAfe2  ~CA5eJJk180  .  A-14 
= ~  * per    raa 

8e-i  -  5e\  n 

91 


c 


m5e 


r  w&2      Cffl&!  /„  180 


&?c 


&2  ~~  ^1  K 

F^rriiH  (&=<))/,,,  180 


per    rad 


trim 


c 


m6e 


71 


deg 


A.15 


A.16 


b.  Sample  Longitudinal  Control  Power  Data  Reduction 


Table  A.3  presents  CMARC  solutions  for  two  elevator  deflections. 
Sample  calculations  for  obtaining  elevator  control  power  are  demonstrated  below.  The 
calculations  are  easily  implemented  in  a  spreadsheet  or  with  a  MATLAB  script. 


C 


L 


(0.4641-0.4259).  180 


6e 


(5-0) 


7Z 


=  0.4378  per  radian 


A.17 


C 


DSe 


(0.0178-0.0170).  180 


(5-0) 


K 


=  0.0092  per  radian 


A.18 


C 


(-0.0942-0.0104).  180 


m5e 


(5-0) 


71 


=  -1.199  per  radian 


A.19 


Se< 


trim 


_  (0.0104)^  180  _ 


-1.199       n 


=  -0.50c 


A.20 


RUN# 

Se 

CL 

CD 

CY 

C  m 

C  n 

C  1 

1 

0° 

0.4259 

0.0170 

0.0000 

0.0104 

0.0000 

0.0001 

2 

5° 

0.4641 

0.0178 

0.0000 

-0.0942 

0.0000 

0.0001 

Table  A.3  FROG  UAV  Elevator  Control  Power  Data. 
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3.  Longitudinal  a-dot  Damping  Derivatives 

a.  Longitudinal  a-dot  Derivative  Methods  and  Equations 

Two  motions  selected  are  to  develop  the  longitudinal  rate  damping 
derivatives.  A  sinusoidal  plunging  motion  in  the  z-axis  is  used  to  extract  the  a-dot  terms. 
Then  an  oscillatory  pitching  motion  is  used  to  obtain  the  combined  a-dot  and  pitch  rate 
terms.  The  a-dot  are  then  subtracted  to  yield  the  pitch  rate  damping.  All  motion  for 
FROG  data  gathering  is  conducted  at  a  frequency  of  2n  rad/s,  which  equates  to  a  reduced 
frequency  of  k=0.0595  for  this  configuration  and  trim  airspeed. 

The  sinusoidal  plunging  motion  illustrated  in  Figure  A. 3  is  used  to  isolate 
the  a-dot  derivatives.  Z-axis  plunging  motion  is  controlled  with  the  CMARC  BINP8B 
input  file  line.  A  frequency  of  In  rad/s  and  an  amplitude  of  C/2=10  inches  are  selected 
for  the  FROG  study.  Greater  amplitude  was  initially  selected,  but  the  induced  angle  of 
attack  caused  the  rigid  wake  to  impact  the  horizontal  tail.  An  example  BINP8B  input  line 
is  shown  below: 

&BINP8B        DXMAX=0.0,  DYMAXO.O,  DZMAX=10.0 

WTX=0.0,  WTY=0.0,  WTZ=6.283,  &END 

The  number  of  time  steps  and  time  step  interval  is  chosen  to  create  a  nice 
sinusoidal  output  through  at  least  two  cycles  of  motion.  In  the  FROG  study,  a  plunging 
motion  frequency  of  co=27i  rad/s  or  one  cycle/sec  is  utilized.  Fifty  time  steps  are  chosen 
with  an  interval  of  0.05  seconds,  which  creates  2.5  cycles  of  plunging  motion.  After  the 
solution  is  obtained,  the  "total  coefficient"  data  is  extracted  for  plotting.  The  data  can  be 
picked  out  manually,  or  a  data  retrieval  program  can  be  created  for  the  task.  For  the 
FROG  study,  the  data  is  extracted  manually  and  pasted  into  an  "Excel"  spreadsheet  for 
plotting.  In  addition  to  the  CMARC  output,  plunging  motion  phase  angle  and  induced 
angle-of-attack  are  calculated  using  Equations  A.21  and  A.22.  Induced  angle-of-attack  is 
used  to  find  the  phase  angle  of  the  response  with  respect  to  angle-of-attack.  Table  A.4 
shows  representative  FROG  data  for  20  time  steps. 
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180 


$  plunge  -# timestep*^1 ' w plunge  *  degrees 


71 


A.21 


a induced  ~  cos\r  plunge  ) 


■™  plunge 


u, 


ru 


plunge 


180 


7T 


degrees 


A.22 


Where: 


plunge 


u 


induced 


plunging  motion  phase  angle 

induced  angle-of-attack  from  the  plunging  motion 


#timestep  -  time  step 

dt         -  time  step  interval 

^plunge    -  plunging  frequency 

Aplunge    -  plunging  amplitude 

U0         -  reference  free  stream  velocity 


Step 

<t>  (deg) 

O  induced 

CL 

CD 

CY 

C_m 

C_n 

c_i 

0 

0 

-3.41 

0.0001 

0.0327 

0.0000 

0.0173 

0.0000 

0.0000 

1 

18 

-3.24 

0.1318 

0.0300 

0.0000 

0.0744 

0.0000 

0.0000 

2 

36 

-2.76 

0.1928 

0.0299 

0.0000 

0.0485 

0.0000 

0.0000 

3 

54 

-2.00 

0.2631 

0.0316 

0.0000 

0.0294 

0.0000 

0.0000 

4 

72 

-1.05 

0.3451 

0.0318 

0.0000 

0.0218 

0.0000 

0.0000 

5 

90 

0.00 

0.4336 

0.0297 

0.0000 

0.0169 

0.0000 

0.0000 

6 

108 

1.05 

0.5216 

0.0252 

0.0000 

0.0126 

0.0000 

0.0000 

7 

126 

2.00 

0.6009 

0.0189 

0.0000 

0.0087 

0.0000 

0.0000 

8 

144 

2.76 

0.6636 

0.0123 

0.0000 

0.0067 

0.0000 

0.0000 

9 

162 

3.24 

0.7029 

0.0069 

0.0000 

0.0089 

0.0000 

0.0000 

10 

180 

3.41 

0.7160 

0.0042 

0.0000 

0.0139 

0.0000 

0.0000 

11 

198 

3.24 

0.7009 

0.0048 

0.0000 

0.0232 

0.0000 

0.0000 

12 

216 

2.76 

0.6594 

0.0085 

0.0000 

0.0344 

0.0000 

0.0000 

13 

234 

2.00 

0.5952 

0.0142 

0.0000 

0.0474 

0.0000 

0.0000 

14 

252 

1.05 

0.5141 

0.0204 

0.0000 

0.0596 

0.0000 

0.0000 

15 

270 

0.00 

0.4254 

0.0255 

0.0000 

0.0704 

0.0000 

0.0000 

16 

288 

-1.05 

0.3370 

0.0286 

0.0000 

0.0762 

0.0000 

0.0000 

17 

306 

-2.00 

0.2582 

0.0297 

0.0000 

0.0776 

0.0000 

0.0000 

18 

324 

-2.76 

0.1964 

0.0292 

0.0000 

0.0734 

0.0000 

0.0000 

19 

342 

-3.24 

0.1561 

0.0283 

0.0000 

0.0713 

0.0000 

0.0000 

20 

360 

-3.41 

0.1398 

0.0277 

0.0000 

0.0727 

0.0000 

0.0000 

Table  A.4  Sample  Plunging  Motion  Data  for  20  Time  Steps 
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Figure  A.3  Vertical  Plunging  Motion  for  Obtaining  the  a-dot  Derivatives. 

The  a-dot  stability  derivatives  are  extracted  using  the  methods  outlined  in 
Etkin  [Ref.  13].  First,  the  CL  and  Cm  responses  are  plotted  as  a  function  of  plunging 
phase  angle.  Induced  angle-of-attack  is  also  plotted  on  the  right  hand  vertical  axis. 
Figure  A.4  is  a  representative  plot  of  CMARC  data  for  the  FROG  UAV  study.  The 
response  to  plunging  motion  can  be  broken  up  into  real  and  imaginary  components. 
Figure  A.5  is  a  graphical  representation  of  the  real  and  imaginary  components.  The  out- 
of-phase  (imaginary)  portion  of  CL  or  Cm  is  the  a-dot  contribution.  It  is  normalized  by 
dividing  by  the  amplitude  of  ajnduced  and  the  reduced  frequency.  The  phase  angle  is 
measured  between  the  lift  or  pitching  moment  response  and  the  induced  angle-of-attack. 
Equations  A.23  through  A.27  are  used  to  solve  for  CLad0,  and  C^^. 
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NPS  FROG  UAV 
Vertical  Plunging  Motion:  coplunge=2pi  rad/s,  A-cbar/2 
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Figure  A.4  Representative  Vertical  Plunging  Motion  Data  for  the  FROG  UAV. 


Figure  A.5  Complex  Response  to  Vertical  Plunging  Motion. 
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The  phase  angle  between  the  CL  (or  C  J  response  and  the  induced  angle- 
of-attack  is  measured  graphically  or  through  numerical  curve  fit  techniques.  If  the 
graphical  method  is  used,  the  graph  scales  may  be  narrowed  around  the  area  of  interest  to 
facilitate  the  phase  angle  measurement.  For  instance,  Figure  A.4  shows  that  lift 
coefficient  and  angle-of  attack  are  nearly  in-phase.  When  the  x-axis  scale  is  zoomed  into 
the  peak  of  the  curves,  one  can  see  that  lift  coefficient  leads  angle-of-attack  by  one 
degree.  Alternatively,  the  phase  angle  may  be  measured  from  where  the  parameter 
crosses  the  trim  value. 


>\Cl0\ 


Ct     =         °  A.23 

^a        a0k 

=  _L_^£J  A.24 


m 


a        a0k 


a0=—0>  plunge  A-25 

U  o 

4CI0  }=  CLo  sinfo)  and    /[c^  ]=  CmQ  sin^)  A.26 

k  =  jf&  plunge  A-27 

U  o 


Where:  (f>  -  Phase  angle  between  CL  (or  C  J  response  and  ainduced 

conduced  _  Induced  angle-of-attack  from  the  plunging  motion 

CLo  -  Amplitude  of  lift  coefficient  response 

Cmo  -  Amplitude  of  pitching  moment  coefficient  response 

a0  -  Amplitude  of  the  induced  angle-of-attack 

k  -  Reduced  frequency 

U0  -  Reference  free  stream  velocity 

b.  Sample  Longitudinal  a-dot  Data  Reduction 
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Sample  data  reduction  is  presented  below  in  Equations  A.28  through  A.32. 
The  phase  angle  between  the  response  and  the  induced  angle-of-attack  was  measured 
graphically  by  zooming  the  graph  axes  to  expand  the  area  of  interest. 


c/  20   in, 

l[cLo ]=  CLo  rinGO-fr-7160-0-'398* * sm(t°)  =  0.005028 


CmQ  ]=  Cmo  sm{(f>)  =  (0.0776  -  0.0087)*  sin(-  140°)=  -0.0222 
A  20  in/  2 


a  plunge  =  —  * lTC  =  0.0595    rad* =  3.41° 


U, 


\056in/s 


180 

71 


ACL0\  0.005028  , 

Lr       — — —  1. 

La,  ~    ;,         nAcr»c*AA«oc 


'or        a0k       0.0595*0.0595 


=  1.420 


A.28 


A.29 


A.30 


A.31 


=  /jCj=      -0.0222 


w 


or        a0£        0.0595*00595 


A.32 
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4.  Longitudinal  Pitch  Rate  Damping  Derivatives 

a.  Longitudinal  Pitch  Rate  Damping  Methods  and  Equations 

A  sinusoidal  plunging  motion  in  the  z-axis  was  used  to  extract  the  a-dot 
terms.  Now,  an  oscillatory  pitching  motion  is  used  to  obtain  the  combined  a-dot  and 
pitch  rate  terms.  The  a-dot  influence  is  then  subtracted  to  yield  pitch  rate  damping.  All 
motion  for  FROG  data  gathering  is  conducted  at  a  frequency  of  2tc  rad/s,  which  equates 
to  a  reduced  frequency  of  k=0.0595  for  this  configuration  and  trim  airspeed. 

The  sinusoidal  pitch  rate  motion  illustrated  in  Figure  A.6  is  used  to  isolate 
the  combined  a-dot  and  pitch  rate  influence.  Oscillating  pitch  motion  is  controlled  with 
the  CMARC  BINP8A  input  file  line.  A  frequency  of  2n  rad/s  and  an  amplitude  of  ±1.5° 
are  selected  for  the  FROG  study.  Any  larger  pitch  amplitude  would  cause  the  wake  to 
impact  the  horizontal  tail.  An  example  BINP8A  input  line  is  shown  below.  Note  that 
pitch  amplitude  is  in  degrees  and  frequency  is  in  rad/sec: 

&BINP8A       PHIMAX=0.0,  THEMAX=1.5,  PSIMAX=0.0 

WTX=0.0,  WTY=6.283,  WTZ=0.0,  &END 


Figure  A.6  Pitching  Motion  for  Obtaining  the  Pitch  Rate  Damping  Derivatives. 

The  number  of  time  steps  and  time  step  interval  is  chosen  to  create  a  nice 
sinusoidal  output  through  at  least  two  cycles  of  motion.   In  the  FROG  study,  a  pitching 
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motion  frequency  of  (o=2n  rad/s  or  one  cycle/sec  is  utilized.  Fifty  time  steps  are  chosen 
with  an  interval  of  0.05  seconds.  This  combination  yields  2.5  cycles  of  pitching  motion. 
After  the  solution  is  obtained,  the  "total  coefficient"  data  is  extracted  for  plotting.  The 
data  can  be  picked  out  manually,  or  a  data  retrieval  program  can  be  created  for  the  task. 
For  the  FROG  study,  the  data  is  extracted  manually  and  pasted  into  an  "Excel" 
spreadsheet  for  plotting.  In  addition  to  the  CMARC  output,  pitch  motion  phase  angle  and 
angle-of-attack  are  calculated  using  Equations  A.33  and  A.34.  Lift  and  pitching  moment 
phase  angle  is  obtained  with  respect  to  angle-of-attack.  Usually,  the  second  motion  cycle 
is  used  to  ensure  that  all  start  up  transients  have  subsided.  Table  A. 5  shows 
representative  FROG  data  for  20  time  steps. 


180 
t  pitch  =#timesteP  *dt  *  &  pitch  * degrees  A.33 

a  =  6  =  sm{<f> Pitch  )Apitch  degrees  A.34 

Where:  (j>pnCh     -  pitching  motion  phase  angle 

a  -  Angle-of-attack  from  the  pitching  motion 

#timesteP  -  time  step 

dt         -  time  step  interval 

©pitch    -  pitching  frequency 

Apitch    -  pitch  amplitude 


100 


Step 

<t»  (deg) 

a 

CL 

CD 

CY 

C_m 

C_n 

c_i 

0 

0 

0.00 

0.0540 

0.0273 

0.0000 

0.0828 

0.0000 

0.0000 

1 

18 

0.46 

0.4545 

0.0364 

0.0000 

-0.0462 

0.0000 

0.0000 

2 

36 

0.88 

0.4961 

0.0344 

0.0000 

0.0303 

0.0000 

0.0000 

3 

54 

1.21 

0.5312 

0.0347 

0.0000 

0.0226 

0.0000 

0.0000 

4 

72 

1.43 

0.5491 

0.0351 

0.0000 

0.0270 

0.0000 

0.0000 

5 

90 

1.50 

0.5541 

0.0348 

0.0000 

0.0310 

0.0000 

0.0000 

6 

108 

1.43 

0.5435 

0.0338 

0.0000 

0.0447 

0.0000 

0.0000 

7 

126 

1.21 

0.5227 

0.0321 

0.0000 

0.0535 

0.0000 

0.0000 

8 

144 

0.88 

0.4925 

0.0302 

0.0000 

0.0618 

0.0000 

0.0000 

9 

162 

0.46 

0.4557 

0.0281 

0.0000 

0.0692 

0.0000 

0.0000 

10 

180 

0.00 

0.4167 

0.0262 

0.0000 

0.0727 

0.0000 

0.0000 

11 

198 

-0.46 

0.3787 

0.0247 

0.0000 

0.0741 

0.0000 

0.0000 

12 

216 

-0.88 

0.3457 

0.0236 

0.0000 

0.0720 

0.0000 

0.0000 

13 

234 

-1.21 

0.3218 

0.0230 

0.0000 

0.0648 

0.0000 

0.0000 

14 

252 

-1.43 

0.3075 

0.0227 

0.0000 

0.0585 

0.0000 

0.0000 

15 

270 

-1.50 

0.3055 

0.0228 

0.0000 

0.0499 

0.0000 

0.0000 

16 

288 

-1.43 

0.3158 

0.0233 

0.0000 

0.0402 

0.0000 

0.0000 

17 

306 

-1.21 

0.3369 

0.0241 

0.0000 

0.0323 

0.0000 

0.0000 

18 

324 

-0.88 

0.3674 

0.0254 

0.0000 

0.0246 

0.0000 

0.0000 

19 

342 

-0.46 

0.4039 

0.0271 

0.0000 

0.0191 

0.0000 

0.0000 

20 

360 

0.00 

0.4429 

0.0290 

0.0000 

0.0153 

0.0000 

0.0000 

Table  A.5  Sample  Pitch  Motion  Data  for  20  Time  Steps 

The  Cl  and  Cm  responses  are  plotted  as  a  function  of  pitch  motion  phase 
angle.  Angle-of-attack  is  also  plotted  on  the  right  hand  vertical  axis.  Figure  A.7  is  a 
representative  plot  of  CMARC  data  for  the  FROG  UAV  study.  The  out-of-phase 
(imaginary)  portion  of  Cl  or  Cm  is  the  combined  a-dot  and  pitch  damping  contribution. 
As  seen  in  Equations  A.35  and  A.36,  the  a-dot  contribution  is  subtracted  from  the  total  to 
yield  the  pitch  damping  influence.  The  pitch  rate  contribution  is  normalized  by  pitch  rate 
and  t*.  The  phase  angle  is  measured  between  the  lift  or  pitching  moment  response  and 
angle-of-attack.  Equations  A.35  through  A.44  are  used  to  solve  for  CLq  and  Cmq. 
Equations  A.35  and  A.36  assume  that  the  a-dot  and  pitch  rate  contributions  are  in  phase 
with  each  other.  The  a-dot  contribution  is  small,  so  there  isn't  harm  done  if  the  two 
contributions  are  somewhat  out-of-phase.  Consequently,  the  a-dot  contribution  is 
calculated  based  on  the  maximum  a-dot  rate  observed. 


l[CLd+\=l[CLJ+l[Cr\ 


A.35 


4C^+J=4^M^ 


A.36 
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4cwJ=sinWk,J 


A.37 


4^HnWk,J 


A.38 


i[cLd]=cLd*dmax*t' 


J[cmd]=Cmdl*d 


\CLa\=CLc 


max 


*  n  * 


-^J=CV 


A.39 


A.40 


A.41 


C 


m. 


~  ^rtiq      *7max     *■ 


A.42 


4* 


^max       ^max 


*  C     TT 


'deg 


57.3 


deg. 


**VcA  rad/s 


rad 


A.43 


A.44 


Where: 


[Cl  a-dot+q] 
L^-rn  a-dot+qj 
I[Cl  a-dot+q] 
H^nia-dot+qJ 
I[CLq] 

I[Cmq] 

I[Cl  a-dot] 
l\S^  ma-dot  J 


-  Amplitude  of  Cl  response  from  pitching  motion 

-  Amplitude  of  Cm  response  from  pitching  motion 

-  Out-of-phase  Cl  due  to  a-dot  and  pitch  rate  damping 

-  Out-of-phase  Cm  due  to  a-dot  and  pitch  rate  damping 

-  Cl  coefficient  contribution  from  pitch  rate  damping 

-  Cm  coefficient  contribution  from  pitch  rate  damping 

-  Cl  coefficient  contribution  from  a-dot  damping 

-  Cm  coefficient  contribution  from  a-dot  damping 

-  Characteristic  time 
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Figure  A.7  Representative  Pitch  Motion  Data  for  the  FROG  UAV. 

b.  Sample  Pitch  Rate  Damping  Data  Reduction 

Sample  pitch  rate  damping  data  reduction  is  presented  below  in  Equations 
A.45  through  A. 54.  The  phase  angle  between  angle-of-attack  and  the  Cl  or  Cm  response 
was  measured  graphically  by  zooming  the  graph  axes  to  expand  the  area  of  interest. 


1056^ 

sec 


=  0.009470  sec 


A.45 


W  =  «max  =      ^eg/     *****  =  ^T*lK   ra^  =  01645rad/s        A'46 

/  rad 


7[C/.d]=CIa  +drmax  *'*  =l-420*01645*0-009470  =  a002212 


A.47 
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l[cmd  ]=Cmd  *drmax  *t*  =  -6.264*0.1645*0.009470  =  -0.009758  A.48 

^Ld+q  \=M<fi)[CLd+q  J=sin(6°)*  [0.1234]  =  0.0129  A.49 

l]frnd+q  ]=  sin(<P)[Cmd+q  ]=  sin(-100°)  *  [0.02853]  =  -0.02810  A.50 

l[cL   J=4CZ»+    J- /[Q^  0.0129 -0.002212  =  0.01069  A.51 

![Cmq  ]=  4C"^+<7  J"  7lC««  J=  "0-02810  -  (-0.009758)  =  -0.01834  A.52 

Cl    =    I[Cl<[=         0Q1Q69  =  6.862  A.53 

1      a        *t        0.1645*0.009470 
'/max      ' 

_    4fi»J    _        -0.01834        _    n7g  A54 

w<7      ^        *,*      0.1645*0.009470 

'/max     ' 
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C.         LATERAL-DIRECTIONAL  STABILITY  DERIVATIVES 

This  section  will  describe  CMARC  methods  for  the  development  of  lateral- 
directional  stability  derivatives  to  include  the  static,  rate  damping  and  control  power 
derivatives.  The  results  obtained  from  CMARC  should  be  spot  checked  against  classical 
design  calculations. 

Of  note,  CMARC  uses  the  semi-span  to  normalize  rolling  and  yawing  moment 
coefficients.  Most  texts  on  stability  and  control,  including  Roskam  [Ref.  12]  and  Etkin 
[Ref.  13],  normalize  rolling  and  yawing  moments  by  span.  This  study  will  normalize  roll 
and  yaw  moments  by  span.  Table  A.  1  summarizes  the  factors  for  normalizing  moments 
and  angular  rates.  All  rolling  and  yawing  moment  coefficients  presented  in  this  study 
have  been  normalized  with  span  by  dividing  the  CMARC  output  by  a  factor  of  two. 
Table  A.l  also  indicates  the  characteristic  time,  t*,  employed  for  angle  rate  data 
reduction. 

For  the  lateral-directional  analysis,  both  sides  of  the  body  must  be  modeled.  The 
easiest  way  to  do  this  is  to  activate  the  symmetric  patch  toggle  for  each  patch 
(IPATSYM=1).  Then,  turn  off  symmetric  calculations  (RSYM=1.0).  This  creates 
symmetric  patches  around  the  y=0  plane,  allowing  CMARC  to  perform  asymmetric 
calculations  around  the  entire  body.  Processing  times  are  significantly  increased 
compared  to  the  symmetric  case.  Wakes  must  be  defined  for  each  side  of  the  model 
based  on  the  new  patch  numbers.  In  addition,  don't  forget  to  update  the  adjacent  patch 
numbers  of  the  wingtips. 

The  vertical  tail  is  activated  for  the  lateral-directional  study.  This  requires  that  the 
wing  and  fuselage  group  be  run  separately  from  the  horizontal  and  vertical  stabilizer 
group.  This  de-coupling  is  required  to  keep  the  wing  wake  from  hitting  the  vertical 
stabilizer.  The  two  sets  of  results  are  summed  to  get  the  total  response.  Typically,  the 
sidewash  derivative,  ds/dp  is  small,  making  separate  solutions  feasible.  The  main  draw 
back  is  that  separate  solutions  will  not  capture  the  finer  interactions  that  a  complete 
model  would  capture.  However,  many  aircraft  configurations  will  not  require  the  de- 
coupling of  the  empennage  surfaces.  If  the  configuration  permits,  run  the  lateral- 
directional  test  cases  as  a  complete  airframe  model.  Figure  A.8  shows  the  FROG  model 
configurations  used  to  find  the  lateral-directional  derivatives. 
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It  should  be  noted  that  the  CMARC  wind  axes  are  modeled  with  X-aft/Z-up,  vice 
X-forward/Z-down  for  the  typical  stability  axes  system.  Positive  yaw  angle  in  CMARC 
creates  positive  sideslip  in  the  stability  axes  system.  Care  must  be  taken  to  reverse  the 
signs  of  the  appropriate  coefficients  to  convert  from  the  CMARC  wind  axes  to  the 
stability  axis  system  shown  in  Figure  A.l. 


Figure  A.7  FROG  UAV  Lateral-Directional  Models. 

1.  Static  Lateral  Directional  Stability  Derivatives 

a.  Static  Lateral-Directional  Stability  Equations 

Three  lateral-directional  stability  derivatives  are  measured  with  the  first 
two  runs  of  the  CMARC  model.  The  model  is  first  analyzed  at  the  estimated  trim 
condition  with  zero  sideslip.  The  baseline  solution  should  be  compared  to  the 
longitudinal  baseline.  If  the  lift  coefficient  differs  significantly,  the  wake  modification 
assumptions  may  need  to  be  revisited.  The  model  is  then  checked  for  lateral-directional 
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balance  at  zero  sideslip  angle  (yaw  angle  =0).  The  side  force,  rolling  and  yawing 
coefficients  should  be  zero  when  a  trial  run  is  performed  at  zero  sideslip.  If  lateral- 
directional  forces  or  moments  are  present,  the  model  and  wake  geometry  should  be 
checked  for  symmetry.  Visualization  with  LOFTSMAN  will  assist  in  spotting  the 
problem. 

For  the  FROG  case,  0^=0°  and  (3=0°  is  selected  for  the  baseline.  A 
second  CMARC  run  is  then  conducted  with  yaw  angle  incremented  by  one  or  two 
degrees.  CY,  C,  and  Cn  are  then  extracted  manually  from  the  data  files.  The  slope  of  CY, 
C,  and  Cn  versus  sideslip  provides  the  CYp,  C,p  and  CnP  longitudinal  derivatives. 
Additional  solutions  can  be  obtained  to  check  for  slope  linearity.  Equations  A.55  through 
A. 5 7  are  used  for  these  calculations: 


CYo  =  Cr  *180 

n 

C,   ,180 


■''-& 


71 


per  radian 
per  radian 


C[P  =  A/? 

C„a  =  — —  *  -    -  per  radian 


'nP      A/?°      n 


A.55 
A.56 

A.57 


b. 


Sample  Static  Lateral-Directional  Data  Reduction 


Table  A.6  presents  CMARC  solutions  for  the  FROG  UAV  at  two  angles- 
of-attack.  Data  is  for  the  summed  contribution  of  the  wing/fuselage  and  tail  group 
models.  Roll  and  yaw  moments  are  already  divided  by  a  factor  of  two.  This  normalizes 
by  span  to  compensate  for  the  CMARC  output,  which  is  normalized  by  semi-span. 
Sample  calculations  for  obtaining  the  static  derivatives  are  demonstrated  below.  The 
calculations  are  easily  implemented  in  a  spreadsheet  or  with  a  MATLAB  script. 


RUN# 

3 

CL 

CD 

CY 

C  m 

C  n 

C  1 

1 

0° 

0.4432 

0.0145 

0.0000 

0.0048 

0.0000 

0.0000 

2 

2° 

0.4452 

0.0148 

-0.0087 

0.0038 

0.0022 

-0.0022 

Table  A.6  FROG  UAV  Static  Lateral-Directional  Stability  Data  at  ct^^O0. 
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Sample  Calculations: 

_  CY  *180     -0.0087,,  180  deg        .  _  ...  ,.  .  __ 

Cyn  = — —  * = z * —  =  -0.2493  per  radian  A.55 

YP      AJ3°      n  2  n    rad 

_  Ci   ^  180     -0.0022  ^180  deg      .  __ft  ,. 

Cj  0  =  — —  * = 5 * —  =  0.0630  per  radian 

lP      A/?°      7t  2  n    rad 

A.56 

=  ^,180  =  M022t  ISOdeg 
"^      A/?°      n  £         n   rad 


3.  Lateral-Directional  Control  Power  Stability  Derivatives 

a.  Lateral-Directional  Control  Power  Equations 

The  rudder  control  power  derivatives  are  obtained  by  substituting  a  0° 
deflection  vertical  stabilizer  patch  for  one  with  positive  rudder  deflection.  Only  one  run 
is  required.  For  the  FROG  UAV  study,  +5°  (TEL)  deflection  is  used.  The  difference 
between  the  trim  condition  and  the  deflected  value  is  divided  by  the  rudder  deflection  as 
shown  below: 

r       [cYsa2  -Cy&t)mm  ,  A„ 

CVo    = —  * per    rad  A.58 

00  Sa2  —5a\  7i 

r  [Cl5a2  ~C'&i^l80 

C/o    = —  * per    rad  A.59 

00  Sa2  -&i  n 


lC*&2       CnSai  J ^  180 
Sa2  —  5a-[  7t 


cnSa  =  ^2      ;™i'*i^  per  r^  A.60 


r        lc^  ~Cy*J«i80  .  AM 

Cyo    = —  * per    rad  A.ol 

**  8r2  -  Sr^  7i 
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Sr2  —  Sri  n 


C\~   = ^ p         * per    rad  A.62 


C„  r.    = — —  * per    rad  A.63 

dr  Sr2  -  S^  k 


b.         Sample  Lateral-Directional  Control  Power  Data  Reduction 

Tables  A.7  and  A.8  present  CMARC  solutions  for  aileron  and  rudder 
deflections.  Sample  calculations  for  obtaining  aileron  and  rudder  control  power  are 
demonstrated  below.  The  calculations  are  easily  implemented  in  a  spreadsheet  or  with  a 
MATLAB  script. 

Q*.  =[Cy&2  -Cr*l™J-°-™\-°).m  =  -o.o103  A.64 

rSa  Sa2-Sai  7i  (5-0)  it 

=  (Q&2  -Cl&1  J,  180  JO.01695-0),  180  =  Q  1943 
&  &2-&1  *  (5-0)  n 

Cn*.ynSn  "C"&1^  J-0-O0105-O)tlgg  =  -0,Q120  A.66 

*  &2-&1  f  (5-0)  JT 

=  frfr  -C^  J.  180  =  (0.0081-0).180  =  0Q928  A6? 

^r  Jr2-^i  *  (5-0)  ;r 

=  ic^  -c^^M  =  (QQQQ35-o),i8o  =  00040  A68 

/Jr  Sr2-Sri  n  (5-0)  ;r 

=  lC^  -C^J,180  =  (-0.00395-0), 180_    ^^  A  69 

"*  £r2-&i  n  (5"°)  * 
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RUN# 

5a 

CL 

CD 

CY 

C  m 

C  n 

C  1 

1 

0° 

0.4259 

0.0170 

0.0000 

0.0104 

0.0000 

0.0000 

2 

5° 

0.4005 

0.0173 

-0.0018 

0.0199 

-0.00105 

0.01695 

Table  A.7  FROG  UAV  Aileron  Control  Power  Data  at  a^^O0. 


RUN# 

5r 

CL 

CD 

CY 

C  m 

C  n 

C  1 

1 

0° 

0.4259 

0.0170 

0.0000 

0.0104 

0.0000 

0.0000 

2 

5° 

0.4258 

0.0170 

0.0081 

0.0107 

-0.00395 

0.00035 

Table  A.8  FROG  UAV  Rudder  Control  Power  Data  at  a^^O0. 


3.  Yaw  Rate  Damping  Derivatives 

a.  Yaw  Rate  Derivative  Methods  and  Equations 

Only  one  motion  is  required  for  the  yaw  rate  terms.  The  P-dot  terms  are 
generally  considered  negligible.  Therefore,  a  sideways  plunging  motion  is  not  required. 
The  yaw  rate  terms  are  yielded  directly  from  an  oscillating  yawing  motion  as  depicted  in 
Figure  A.8.  Yawing  motion  data  for  the  FROG  is  gathered  at  a  frequency  of  2n  rad/s, 
which  equates  to  a  reduced  frequency  of  k=0.369  for  this  configuration  and  trim  airspeed. 
Oscillating  yaw  motion  is  controlled  with  the  CMARC  BINP8A  input  file  line.  An 
amplitude  of  ±2°  is  selected  for  the  FROG  study.  An  example  BENP8A  input  line  is 
shown  below.  Note  that  yaw  amplitude  is  in  degrees  and  frequency  is  in  rad/sec: 


&BINP8A 


PHIMAX=0.0, 
WTX=0.0, 


THEMAX=1.5, 
WTY=6.283, 


PSIMAX=0.0 
WTZ=0.0, 


&END 
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Figure  A.8  Yawing  Motion  for  Obtaining  the  Yaw  Rate  Derivatives. 

The  number  of  time  steps  and  time  step  interval  is  chosen  to  create  a  nice 
sinusoidal  output  through  at  least  two  cycles  of  motion.  In  the  FROG  study,  a  yawing 
motion  frequency  of  (o=2n  rad/s  or  one  cycle/sec  is  utilized.  Fifty  time  steps  are  chosen 
with  an  interval  of  0.05  seconds,  which  creates  2.5  cycles  of  yawing  motion.  After  the 
solution  is  obtained,  the  "total  coefficient"  data  is  extracted  for  plotting.  The  data  can  be 
picked  out  manually,  or  a  data  retrieval  program  can  be  created  for  the  task.  For  the 
FROG  study,  the  data  is  extracted  manually  and  pasted  into  an  "Excel"  spreadsheet  for 
plotting.  As  discussed  earlier,  the  FROG  model  was  split  into  wing/fuselage  and  tail 
groups  for  the  lateral-directional  study.  Both  sets  of  data  are  pasted  into  the  Excel 
spreadsheet  and  summed.  The  total  value  is  then  plotted.  In  addition,  yawing  motion 
phase  and  sideslip  angles  are  also  calculated  using  Equations  A.70  and  A.71.  Yawing 
motion  phase  angle  is  plotted  on  the  x-axis  and  sideslip  is  plotted  on  the  right  hand  axis. 
Table  A.9  shows  representative  FROG  yaw  rate  data  for  20  time  steps. 


180 
0vaw  =#timesteD  *dt  *mvaw* degrees 


yaw     ^timestep 


yaw 


A.70 


n 
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p  =  sin(^vavv)  *  Ayaw  degrees 


A.71 


Where:  <|)yaw  -  yawing  motion  phase  angle 

P  -  Sideslip  angle  from  yawing  motion 

#t,mestep  -  time  step 

dt  -  time  step  interval 

coyaw  -  Yawing  frequency 

Ayaw  -  Yaw  angle  amplitude 


Step 

<t>  (deg) 

P 

CL 

CD 

CY 

C_m 

C_n 

c_i 

0 

0 

0.00 

0.0534 

0.0111 

-0.0005 

0.085 

0.0003 

-0.0002 

1 

18 

0.62 

0.4004 

0.0195 

-0.007 

-0.0107 

0.0023 

-0.0023 

2 

36 

1.18 

0.4329 

0.0161 

-0.0087 

-0.0009 

0.0026 

-0.0026 

3 

54 

1.62 

0.4388 

0.0152 

-0.0098 

0.0019 

0.0027 

-0.0027 

4 

72 

1.90 

0.4416 

0.0149 

-0.0098 

0.0027 

0.0026 

-0.0026 

5 

90 

2.00 

0.4428 

0.0149 

-0.0089 

0.0033 

0.0022 

-0.0022 

6 

108 

1.90 

0.4435 

0.0147 

-0.0071 

0.0039 

0.0016 

-0.0016 

7 

126 

1.62 

0.444 

0.0146 

-0.0046 

0.0044 

0.0008 

-0.0008 

8 

144 

1.18 

0.4444 

0.0146 

-0.0017 

0.0049 

0.0000 

0.0000 

9 

162 

0.62 

0.4447 

0.0144 

0.0014 

0.0051 

-0.0009 

0.0009 

10 

180 

0.00 

0.445 

0.0144 

0.0043 

0.0048 

-0.0016 

0.0016 

11 

198 

-0.62 

0.4452 

0.0145 

0.0069 

0.0044 

-0.0022 

0.0022 

12 

216 

-1.18 

0.4453 

0.0146 

0.0087 

0.004 

-0.0026 

0.0026 

13 

234 

-1.62 

0.4452 

0.0147 

0.0098 

0.0035 

-0.0027 

0.0027 

14 

252 

-1.90 

0.4451 

0.0147 

0.0098 

0.0034 

-0.0026 

0.0025 

15 

270 

-2.00 

0.4449 

0.0148 

0.0089 

0.0036 

-0.0022 

0.0021 

16 

288 

-1.90 

0.4448 

0.0146 

0.0071 

0.0041 

-0.0016 

0.0015 

17 

306 

-1.62 

0.4449 

0.0146 

0.0046 

0.0045 

-0.0008 

0.0008 

18 

324 

-1.17 

0.445 

0.0145 

0.0017 

0.005 

0.0000 

-0.0001 

19 

342 

-0.62 

0.4452 

0.0144 

-0.0014 

0.0051 

0.0009 

-0.0009 

20 

360 

0.00 

0.4454 

0.0144 

-0.0044 

0.0049 

0.0016 

-0.0016 

Table  A.9  Sample  FROG  UAV  Yawing  Motion  Data  for  20  Time  Steps 

The  CY,  C,  and  Cm  responses  are  plotted  as  a  function  of  yawing  phase 
angle.  Sideslip  is  also  plotted  on  the  right  hand  vertical  axis.  Figure  A.9  is  a 
representative  plot  of  CMARC  data  for  the  FROG  UAV  study.  The  out-of-phase 
(imaginary)  portion  of  CY,  C,  or  Cm  is  due  to  the  yaw  rate  damping  contribution.  The  |3- 
dot  contribution  is  considered  small  and  is  ignored.  The  phase  angle  is  measured 
between  the  coefficient  response  and  sideslip  angle.     In  this  case,  the  zero  crossing 
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method  of  measuring  phase  angle  proves  to  the  easiest  and  most  accurate.  The  phase 
angle  is  simply  measured  from  the  parameter  to  sideslip  angle  where  they  cross  the  x-axis 
in  the  same  direction.  The  yaw  rate  damping  contribution  is  normalized  by  maximum 
yaw  rate  and  t*.  Equations  A. 72  through  A. 76  are  used  to  solve  for  CYr,  Clr  and  Cm. 

CYr  =    7        *  .    where    l[CY  ]  =  sin(^flM, \cY  ]  A.72 


r        *  t 
'max     * 


C/r  = — [—L^r,    where   l[Ct  ]  =  sin(^aw  |Q  ]  A.73 


'max 


Cnr  =     7^l  ,    where    l[Cn ]  =  sin(^avv \cn ]  A.74 

rmax     ^ 

>"max  = A^7—*0)yaw    rad/s  A-75 


57.3' 

/  rad 

b/ 
t*  =  ££-  sec  A.76 


C/i 


o 


Where: 


[CY]  -  Amplitude  of  CY  response  from  yawing  motion 

[C,]  -  Amplitude  of  C,  response  from  yawing  motion 

[Cn]  -  Amplitude  of  Cn  response  from  yawing  motion 

I[CY]  -  CY  coefficient  contribution  from  yaw  rate  damping 

I[C,]  -  C,  coefficient  contribution  from  yaw  rate  damping 

I[CJ  -  Cn  coefficient  contribution  from  yaw  rate  damping 

t*  -  Characteristic  time 
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YAW  OSCILLATIONS  -  NPS  FROG  UAV 
co=2ji  rad/s,  A=±2° 
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Figure  A.9  Representative  Yawing  Motion  Plot  for  the  FROG  UAV. 


b.  Sample  Yaw  Rate  Damping  Data  Reduction 

Sample  data  reduction  is  presented  below  in  Equations  A.28  through  A.32. 
The  phase  angle  between  the  response  and  the  sideslip  angle  is  measured  graphically  by 
zooming  in  on  the  x-axis  zero  crossing. 
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124m, 


U0      1056»V 
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max 
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rad 
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rad 
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/sec 


A.78 
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l[CY]    _  sm(# yaw\Cy  ]  _  sin(l54°X0.0099]  _ 


CY   =      l   JJ,= — v  ,.«-»,  J  =  -...v.-.  A J  =  0.3371  A.79 

r      r       */  r       *t  0.2193*0.0587 

'max     *  'max     ' 


/[Q]     _  sin^vaH,|Q]  _  sin(l44°)[0.00265] 


Q    = — ^-4-  =  — v"     \      = — * /lJ  =  0.1210  A.80 

r      r       *t  r       *t  0.2193*0.0587 

'max     '  'max     l 


l[Cn]     _  sin(^flWJ[Cw]  _  Sjn(- 36°Xo.Q0265] 


Cn    =       Ln\= — "     '    -J  =  — v         *— J  =  -o.i210  A.81 

r      r       *t  r       *t  0.2193*0.0587 
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4.  Roll  Rate  Damping  Derivatives 

cl  Roll  Rate  Derivative  Methods  and  Equations 

Gathering  roll  damping  data  is  straightforward.  Unlike  the  pitch  or  yaw 
rate  terms,  there  is  no  change  in  angle-of-attack  or  sideslip  with  the  rolling  motion. 
Therefore,  roll  damping  effects  can  be  measured  with  pure  rolling  motion  around  the  x- 
axis.  As  with  the  yaw  damping  case,  the  wing/fuselage  and  tail  models  are  run 
separately.  The  results  are  then  summed.  A  rigid  wake  seems  to  work  well,  which 
provides  a  solution  after  just  a  few  time  steps. 

Roll  damping  data  for  the  FROG  is  gathered  at  a  20  deg/sec  roll  rate. 
Initially,  a  5  deg/sec  roll  rate  was  selected,  but  the  yawing  moment  coefficient  was  low 
enough  that  numerical  resolution  would  effect  the  results.  The  higher  rate  provided  a 
sufficiently  large  response  for  all  terms.  Pure  rolling  motion  is  controlled  by  the 
PHIDOT  term  in  the  CMARC  BINP8  input  file  line.  An  example  BINP8  input  line  is 
shown  below.  Note  that  roll  rate  is  in  deg/sec: 

&BINP8  ALDEG=0.0,  YAWDEG=0.0,  PHIDOT=20.0 

THEDOT=0.0,  PSIDOT=0.0,  &END 

The  roll  rate  damping  terms  are  obtained  by  normalizing  the  CMARC 
output  by  roll  rate,  p,  and  characteristic  time,  t*.  Equations  A.82  through  A.85  are  used 
for  these  calculations. 
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Cy 

Cy 

p*t* 

'p 

Ci 

* 
p*t 

nP 

Cn 

* 
p*t 

* 
t    = 

b/ 

jLJz-  sec 
^0 

nere: 

[Cy] 

[CJ 

[CJ 

P 

t* 

A.82 


A.83 


A.84 


A.85 


-  Amplitude  of  CY  response  from  roll  rate 

-  Amplitude  of  C,  response  from  roll  rate 

-  Amplitude  of  Cn  response  from  roll  rate 

-  Roll  rate 

-  Characteristic  time 
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d.  Sample  Roll  Rate  Damping  Data  Reduction 

Tables  A.  10  presents  the  CMARC  roll  damping  solution  for  the  FROG 
UAV.  The  data  is  the  summed  contributions  from  the  wing/fuselage  and  tail  surface 
models.  Sample  roll  damping  calculations  presented  below  in  Equations  A. 86  through 
A.89.  The  calculations  are  easily  implemented  in  a  spreadsheet  or  with  a  MATLAB 
script. 


*      /n         62in 
t   =  ^-  =  .     .  =  0.0587  sec 
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A.86 
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0.0010 
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sec     180  rad 


=  0.0488 


*0.0587  sec 


A.87 
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0.00925 


*      p*t*      20^*^^*0.0587sec 
sec     180  rad 


=  0.4514 


A.88 


-0.00045 


C 

"'      P *>"  "  20 M*JLf!f£* 0.0587 sec 
sec     180  ra<i 


=  -0.0220 


A.89 


RUN# 

P 

CL 

CD 

CY 

C  m 

C  n 

C  1 

1 

0  deg/s 

0.4432 

0.0145 

0.0000 

0.0048 

0.0000 

0.0000 

2 

20deg/s 

0.4435 

0.0143 

0.0010 

0.0047 

-0.00045 

-0.00925 

Table  A.10  FROG  UAV  Roll  Damping  Data  at  0^=0°  and  (3=0°. 
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APPENDIX  B 
LOFTSMAN  INPUT  FILES  FOR  FROG  UAV  MODELING 
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APPENDIX  D 
MATLAB  SCRIPT  FOR  DYNAMIC  MODE  ANALYSIS 
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%  froguav.m:  This  MATLAB  script  calculates  FROG  UAV  dynamics.   Eigenvalue  modal  analysis 
%  and  response  to  control  input  is  obtained  for  three  aerodynamic  models. 

% 
clear 

%   Calculate  trim  lift  and  drag  coefficients 

%  FROG  Trim  condition 


%  lbs 

%  slugs 

%  ftA2,  2530  in*2 

%  ft/s  =  1056  in/s 

%  trim  mach  number 

%  ft  =  20  in 

%  ft  =  124  in 

%  slugs/ft A3  800  ft  MSL  std  day 


W  =  67.73; 

m  =  W/32.2; 

S  =   2530/144; 

V  =  88; 

M  =  V/1118; 

cbar  =  20/12; 

b  =  126.5/12; 

rho  =  0.002327; 

cruise_HP  =  5; 

etaP  =  0.35; 

L_D=7; 

CL0  =  0.4295; 

Ixx  =  12.52;  Iyy  =  8.43;  Izz  =  18.55;  1X2=0;  %slugs-ftA2 

q  =  0.5*rho* (V) A2; 

CL  =  W/(q*S) ; 

%   Input  the  non-dimensional  stab  derivatives: 


r********i 


:************! 


r  *****  i 


r**********************. 


r  ******* i 


r*********** 


for  i  =  1 : 3 

if  i==l 

source= ' CMARC ' 

%    Longitudinal  Stability  Derivatives 

D_L_D=W/L_D;CD_L_D  =  D_L_D/ (q*S) ; 

Tr=cruise_HP*etaP*550/V;  CD_Tr=Tr/ (q*S) ; 

CD= (CD_L_D  +CD_Tr) /2 ; 

CLalpha  =  4.845;  CDalpha  =  0.2664;  CMalpha 

CLalphadot  =  1.420;  CMalphadot  =  -6.264;  CDalphadot 

alphatrim  =0;  %  radians 

thetaO  =  alphatrim; 

CLmach  =  0 ;  CDmach  =  0 ;  CMmach  =  0 ; 

CLq  =  6.82;  CMq  =  -11.78;  CDq  =  0; 

CLde  =  0.4378;  CMde  =  -1.199;  CDde  =  0.0092; 

%    Lateral  -  Directional  Stability  Derivatives 

CYbeta  =  -0.2493;  Clbeta  =  -0  .0630,-Cnbeta  =  0.0630; 
Clp  =  -0.4514;  Cnp  =  -0.0220;  CYp  =  0.0488; 
Clr  =  0.1210;  Cnr  =  -0.1210;  CYr  =  0.3370; 
Cldr  =  0.0040;  Cndr  =  -0.0453;  CYdr  =  0.0928; 
Clda  =  0.1943;  Cnda  =  -0.0121;  CYda  =  -0.0206; 

elseif  i==2 

source= ' Classical_Analysis ' 


%  Drag  from  Lift  to  Drag  Ratio 
%  Drag  from  thrust  required  calculation 
%  CD  averaged  from  both  methods 
0.4126; 

I;         %  Assumed  to  be  small 


%  Assumed  to  be  small 

%  Induced  drag  contribution  only 


%    Longitudinal  Stability  Derivatives 
D_L_D=W/L_D;CD_L_D  =  D_L_D/ (q*S) ; 
Tr=cruise_HP*etaP*550/V;  CD_Tr=Tr/ (q*S) 
CD=(CD  L  D  +CD  Tr) /2; 


%  Drag  from  Lift  to  Drag  Ratio 

%  Drag  from  thrust  required  calculation 

%  CD  averaged  from  both  methods 


CLalpha  =  4.82;  CDalpha  =  0.2  53;  CMalpha  =  -0.70  %  d_epsilon/d_alpha  =  0.4  0 
1.56;  CMalphadot  =  -4.14;  CDalphadot  =  0; 


CLalphadot 

alphatrim  =  0; 

thetaO  =  alphatrim; 

CLmach  =  0 ;  CDmach  =  0 ;  CMmach  =  0 ; 

CLq  =  3.89;  CMq  =  -11.39;  CDq  =  0; 

CLde  =  0.39;  CMde  =  -1.04;  CDde  =  0.00, 


%  Assumed  to  be  small 
%  radians 


%  Assumed  to  be  small 
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%    Lateral  -  Directional  Stability  Derivatives 

CYbeta  =  -0.511;  Clbeta  =  -0 . 055 ;Cnbeta  =  0.051; 
Clp  =  -0.300;  Cnp  =  -0.072;  CYp  =  0; 
Clr  =  0.168;  Cnr  =  -0.0762;  CYr  =  0.140; 
Cldr  =  0.0056;  Cndr  =  -0.0341;  CYdr  =  0.081; 
Clda  =  0.213;  Cnda  =  -0.0236;  CYda  =  -0.00; 


%  Assumed  to  be  small 


else 


source= ' Parameter_Analysis ' 

%     Longitudinal  Stability  Derivatives 
CL  =  0.4295;  CD  =  0.0614; 

CLalpha  =  4.0907;  CDalpha  =  0.23;  CMalpha  =  -0.4174; 
CLalphadot  =  1.3877;  CMalphadot  =  -3.7115;  CDalphadot  =  0; 
alphatrim  =  0 ;  %  radians 

thetaO  =  alphatrim; 

CLmach  =  0 ;  CDmach  =  0 ;  CMmach  =  0 ; 
CLq  =  3.35;  CMq  =  -8.8818;  CDq  =  0; 
CLde  =  1.1249;  CMde  =  -1.6208;  CDde  =  0.0676; 

%    Lateral  -  Directional  Stability  Derivatives 

CYbeta  =  -0.9867;  Clbeta  =  -0 . 0942 ; Cnbeta  =  0.1755; 
Clp  =  -0.44  83;  Cnp  =  -0.1077;  CYp  =  0; 
Clr  =  0.2078;  Cnr  =  -0.1212;  CYr  =  0.1096; 
Cldr  =  0.0004;  Cndr  =  -0.0785;  CYdr  =  0.0926; 
Clda  =  0.2387;  Cnda  =  -0.0261;  CYda  =  0; 


%  Assumed  to  be  small 


%  Assumed  to  be  small 


%  Assumed  to  be  small 


end 


%  Calculate  the  Dimensional  Stability  Derivatives 


Xu  =  -q*S/ (m*V) * (2*CD+M*CDmach) ; 

Xalpha  =  q*S/m* (CL-CDalpha) ; 

Xalphadot  =  -q*S/m*cbar/ (2*V) *CDalphadot ; 

Xq  =  -q*S/m*cbar/ (2*V) *CDq; 

Xde  =  -q*S/m*CDde; 


%  Xu  -  1/sec  ,  assumes  dCD/dM=0  (ie  no  Mach  effect) 
%  Xalpha  -  ft/sec"2 
%  Xalphadot  -  ft /sec 
%  Xq  -  ft/sec 
%  Xde  -  ft/sec~2 


Zu  =  -q*S/  (m*V) * (2*CL+M* CLmach)  ;  %  Zu  -  1/sec  ,  assumes  dCL/dM=0  (ie  no  Mach  effect) 

Zalpha  =  -q*S/m* (CD+CLalpha) ;         %  Zalpha  -  ft/sec~2 

Zalphadot  =  -q*S/m*cbar/ (2*V) *CLalphadot ;     %  Zalphadot  -  ft/sec 
Zq  =  -q*S/m*cbar/ (2*V) *CLq  ;  %  Zq  -  ft/sec 

Zde  =  -q*S/m*CLde;  %  Zde  -  ft/sec A2 


Mu  =  q*S*cbar/ (Iyy*V) *CMmach; 

Malpha  =  q*S*cbar/Iyy*CMalpha; 

Malphadot  =  q*S*cbar/Iyy*cbar/ (2*V) *CMalphadot ; 

Mq  =  q*S*cbar/Iyy*cbar/ (2*V) *CMq; 

Mde  =  q*S*cbar/Iyy*CMde; 


%  Mu  -  ft/sec 

%  Malpha  -  l/secA2 

%  Malphadot  -  1/sec 

%  Mq  -  1/sec 

%  Mde  -  l/sec*2 


%  Linearized  Longitudinal  4x4  Plant  Matrix 

%  Form  the  An  Bn  and  In  plant  matrices 

An  =  [V*Xu  Xalpha  0  -32 . 17*cos (thetaO) ;  V*Zu  Zalpha  (V+Zq)  -32  .  17*sin (thetaO) 

V*Mu  Malpha  Mq  0;  0  0  1  0] ; 

Bn   =    [Xde   Zde  Mde   0]  '  ; 

In   =    [V   0    0    0;    0     (V-Zalphadot)    0    0;    0    -Malphadot    1    0;    0    0    0    1]  ; 

A   =    inv (In) *An; 
B   =    inv(In)*Bn; 

%  Find  Short  and  Long  Period  Natural  Frequency  and  damping 

out=source 

out= ' Longitudinal   mode   E-values1 

P   =   poly (A) ; 
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R  =  roots (P) 

out=' Short  Period' 

phi  =  angle (R) ; 

Zl(i)  =  (l/(l+tan(phi(l) )*2) )*.5 

Wnl(i)  =  real(R(l) )/(-Zl(i) ) 

Wdl(i)  =  (l-Zl(i)*2)*.5*Wnl(i) ; 

Tdl(i)  =  2*pi/Wdl(i) ; 

out='Long  Period' 

Z3(i)  =  (l/(l+tan(phi(3) )*2))*.5 

Wn3(i)  =  real(R(3) )/(-Z3(i) ) 

Wd3(i)  =  (1-Z3(i) "2) A.5*Wn3 (i) ; 

Td3(i)  =  2*pi/Wd3 (i) ; 

out='Step  response' 

Cl=  [1  000; 0100, -0010, -0001];  D1  =  0; 

U=zeros(1001, 1) ; 

U(101:1001) =-2/57.3; 

t=  (0:  .01:10)  ; 

SYS=SS(A,B,C1,D1) ; 

[Y,T]=lsim(SYS,U,t) ; 

tl=t; 

de_s=U*57.3; 

alpha_s( : , i)=Y( : , 2) *57.3; 

pr_S(:,i)=Y(: ,3) *57.3; 

theta_s(: ,i)=Y(: ,4) *57.3; 

Cl=[l  0  0  0;0  1  0  0;0  0  1  0;  0  0  0  1] ;  D1=0; 

U=zeros (1001, 1)  ; 

U(101:160) =-5/57.3; 

U(161:220)=5/57.3; 

t=(0:  .01:10)  ; 

SYS=ss(A,B,Cl,Dl)  ; 

[Y,T]=lsim(SYS,U,t) ; 

tl=t; 

de_d=U*57.3; 

alpha  d( : ,i)=Y( : ,2)*57.3; 

pr_d(: , i)=Y( : ,3) *57.3; 

theta_d( :  , i) =Y( :  ,4) *57  .  3  ; 

%  Lateral -Directional  4x4  Plant  Matrix 

%  Calculate  the  Dimensional  Stability  Derivatives 

Ybeta  =  q*S/m* (CYbeta) ; 

%  Ybeta  -  ft/secA2 

Yp  =  q*S/m*b/ (2*V) *CYp; 

%  Yp  -  ft/sec 

Yr  =  q*S/m*b/ (2*V) *CYr; 

%  Yr  -  ft/sec 

Yda  =  q*S/ra*CYda; 

%  Yda  -  ft/secA2 

Ydr  =  q*S/m*CYdr; 

%  Ydr  -  ft/secA2 

Lbeta  =  q*S*b/Ixx* (Clbeta) ; 

%  Lbeta  -  l/sec~2 

Lp  =  q*S*b/Ixx*b/ (2*V) *Clp; 

%  Lp  -  1/sec 

Lr  =  q*S*b/Ixx*b/(2*V)*Clr; 

%  Lr  -  1/sec 

Lda  =  q*S*b/Ixx*Clda; 

%  Lda  -  l/sec*2 

Ldr  =  q*S*b/Ixx*Cldr; 

%  Ldr  -  l/sec"2 

Nbeta  =  q*S*b/Izz* (Cnbeta) ; 

%  Nbeta  -  l/secA2 

Np  =  q*S*b/Izz*b/ (2*V) *Cnp; 

%  Np  -  1/sec 

Nr  =  q*S*b/Izz*b/ (2*V) *Cnr; 

%  Nr  -  1/sec 

Nda  =  q*S*b/Izz*Cnda; 

%  Nda  -  l/sec*2 

Ndr  =  q*S*b/Izz*Cndr; 

%  Ndr  -  l/sec*2 

%  Linearized  Longitudinal  4x4  Plant  Matrix 

source 
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longmodes= ' Lateral  Directional  4x4  Plant  Matrix1 

%  Form  the  An  Bn  and  In  plant  matrices 

Cn  =  [Ybeta  Yp  32 . 17*cos (thetaO)  (Yr-V) ;  Lbeta  Lp  0  Lr;  0  10  0;  Nbeta  Np  0  Nr] 

Dn  =  [Yda  Lda  0  Nda;  Ydr  Ldr  0  Ndr] ' ; 

Iln  =  [V  0  0  0;  0  1  0  -Ixz/Ixx;  0  0  1  0;  0  -Ixz/Izz  0  1] ; 

C  =  inv(IIn) *Cn; 
D  =  inv(IIn) *Dn; 

%  Find  Dutch  roll,  roll  and  spiral  Natural  Frequency  and  damping 

out=' Lateral  Directional  mode  E-values' 

PI  =  poly (C) ; 

Rl  =  roots (PI) 

if  i==l 

out=' Dutch  Roll  Mode1 

phi  =  angle (Rl) ; 

Z2(i)  =  (l/(l+tan(phi(2) ) A2) )".5 
Wn2(i)  =  real(Rl(2) )/ (-Z2 (i) ) 
Wd2(i)  =  (1-Z2 (i) A2) A.5*Wn2(i) ■ 
Td2(i)  =  2*pi/Wd2(i) ; 

out='Roll  Mode' 
roll(i)  =  (RKD) 
out= ' Spiral  Mode ' 
Spiral (i)  =  (Rl (4) ) ; 

elseif  i==2 

out=' Dutch  Roll  Mode' 

phi  =  angle (Rl) ; 

Z2(i)  =  (l/(l+tan(phi(l) ) "2) ) A.5 

Wn2(i)  =  real  (Rl(D) /(-Z2(i)  ) 

Wd2(i)  =  (1-Z2(i) "2) A.5*Wn2(i) ; 

Td2(i)  =  2*pi/Wd2 (i) ; 

out='Roll  Mode1 

roll(i)  =  (Rl(3)) 

out=' Spiral  Mode' 

Spiral (i)  =  (Rl(4) ) ; 

else 

out=' Dutch  Roll  Mode' 

phi  =  angle (Rl ) ; 

Z2(i)  =  (l/(l+tan(phi(D ) A2) ) A.5 

Wn2(i)  =  real  (RKD  )/(-Z2(i)  ) 

Wd2(i)  =  (1-Z2(i) A2) A.5*Wn2(i) ; 

Td2 (i)  =  2*pi/Wd2 (i) ; 

out='Roll  Mode1 

rolKi)  =  (RK3)) 

out=' Spiral  Mode' 

Spiral (i)  =  (Rl(4)) 

end 

%  Aileron  step 

Cl=[l    0    0    0;0    1    0    0;0    0    1    0;     0    0    0    1] ;    D1=0; 
U=zeros (1001,  2)  ; 
U(101:1001,l) =2/57.3; 
t= (0: .01:10) ; 
SYS  =  ss(C,D,Cl,DD  ; 
[Y,T] =lsim(SYS,U,t) ; 
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tld=t; 

da_sa=U*57.3; 

rr_sa( : , i)  =Y( : ,2) *57.3; 

phi_sa(: ,i)=Y( : ,3) *57.3; 

beta_sa( : , i) =Y ( : , 1) *57 . 3 ; 

%  Aileron  doublet 

Cl=[l  0  0  0;0  1  0  0;0  0  1  0;  0  0  0  1] ;  D1=0; 

U=zeros(1001, 2)  ; 

U(101:175, 1) =-5/57.3; 

U(176 :250, 1) =5/57.3; 

t= (0: .01:10) ; 

SYS=SS(C,D,C1,D1) ; 

[Y,T] =lsim(SYS,U,t) ; 

tld=t; 

da_da=U*57.3; 

rr_da( :,i)=Y(:  ,2) *57.3; 

phi_da(: , i)=Y(: ,3) *57.3; 

beta_da(: ,i)=Y(: ,1)*57.3; 

%  Rudder  step 

Cl=[l  0  0  0;0  1  0  0;0  0  1  0;  0  0  0  1] ;  D1=0; 

U=zeros(101,2) ; 

U(101:1001,2)=2/57.3; 

t=(0: .01:10) ; 

SYS=ss(C,D,Cl,Dl) ; 

[Y,T]=lsim(SYS,U,t)  ; 

tld=t; 

dr_sr=U*57.3; 

beta_sr ( : , i) =Y( : , 1) *57 .3 ; 

lr_sr ( : , i) =Y( : ,2) *57.3; 

rr_sr ( : , i) =Y( : , 4) *57 . 3 ; 

%  Rudder  doublet 

Cl=[l  0  0  0;0  1  0  0;0  0  1  0;  0  0  0  1] ;  D1=0; 

U=zeros(1001,2) ; 

U(101:175,2) =-5/57.3; 

U(176:250,2)=5/57.3; 

t=(0: .01:10) ; 

SYS=SS(C,D,C1,D1) ; 

[Y,T] =lsim(SYS,U,t)  ; 

tld=t; 

dr_dr=U*57.3; 

beta  dr( : , i) =Y( : ,1) *57.3; 

lr_dr(:,i)=Y( : ,2) *57.3; 

rr_dr(:,i)=Y( :  ,4)  *57.3; 

end 

^****#****+*************************+************************************************************** 

%  Plot  Responses 

%  Plot  Elevator  Step  Response 

figure 

subplot (4, 1, 1) ,  plot (tl,de_s) , title ( 'LONGITUDINAL  RESPONSE  TO  A  2  DEGREE  ELEVATOR  STEP 

INPUT' ) ,ylabel ( 'Elevator  (deg) ' ) , axis ( [0  5  -4  4]), grid 

subplot (4,1,2) ,  plot (tl,alpha_s ( : , 1) , ' - . ' , tl,alpha_s ( : ,2) , ' -- ' , tl,alpha_s ( : , 3) , ' - ' ) ,ylabel ( 'A.O.A 

(deg) ' ) ,axis( [0  5  -5  10]), grid 

subplot (4,1,3) ,  plot (tl,pr_s ( : , 1) ,  ' -  .  ' , tl,pr_s ( : , 2) , ' -- ' , tl,pr_s( : ,3) ,  ' - ' ) ,ylabel  (  'q 

(deg/s)  ' ) ,axis(  [0  5  -20  20]), grid 

subplot (4,1,4) ,  plot (tl,theta_s  ( : , 1) , ' - . ' , tl, theta_s( : ,2)  ,  ' -- ' , tl,theta_s ( : , 3) ,  ' - ' ) ,ylabel ( 'Theta 

(deg) ' ) ,grid 

xlabel ( 'Time  (sec)1)  ,axis([0  5  -10  50]),  pause, 

legend ( 'CMARC , 'Classical ',' Parameter  Est      ',0) 

%  Plot  Elevator  Doublet  Response 

figure 
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subplot (4,1,  1)  ,  plot (tl,de_d)  , title ( 'LONGITUDINAL  RESPONSE  TO  A  5  DEGREE  ELEVATOR 

DOUBLET '), ylabel (' Elevator  (deg)  ' ) , axis (  [0  5-10  10]), grid 

subplot (4,1,2)  ,  plot (tl,alpha_d( : ,1) ,  ■ - . • , tl, alpha  d( : ,2) ,  '  -- ' , tl, alpha  d( : ,3) 

, ' - ' ) , ylabel ( 'A.O.A 

(deg) ' ) ,axis( [0  5  -12  12]), grid 

subplot (4,1,3)  ,  plot  (tl,pr_d( :,l),'-.',tl,pr  d(:,2),'--',tl,pr  d(:,3) ,  '-') , ylabel ( 'q 

(deg/s) ' ) ,axis( [0  5  -60  50]), grid 

subplot (4,1,4)  ,  plot (tl, theta_d( :,l),'-.',tl, theta_d ( : , 2) ,  ' -- ' , tl.theta  d(: ,3) 

, ' - ' ) , ylabel ( 'Theta 

(deg) ' ) ,grid 

xlabel('Time  (sec) '),, axis ( [0  5  -10  30]), pause 

legend ( 'CMARC' , 'Classical ',' Parameter  Est       ',0) 

%  Plot  Aileron  Step  Response 

figure 

subplot (4,1,  1)  ,  plot (tld,da_sa( :  ,1)  ),  title ( 'LATERAL-DIRECTIONAL  RESPONSE  TO  A 

2  DEGREE  AILERON 

STEP  INPUT') , ylabel ( 'Aileron  (deg)  '), axis (  [0  5  -4  4]), grid 

subplot (4,1,2)  ,  plot  (tld,rr_sa( : ,1) ,  ' - . ' ,tld,rr  sa (:,  2) ,'--', tld, rr  sa(:,3),'- 

' ) , ylabel ( 'p 

(deg/s) ') ,axis( [0  5  -10  25]), grid 

subplot (4,1,3)  ,  plot (tld,phi_sa( : , 1) ,  ' - .  ' ,tld,phi_sa( : , 2) ,  ' -- ' ,tld,phi_sa( : ,3) 

,'-'), ylabel ( 'Bank 

Angle  (deg)  '), axis (  [0  5  -10  70]), grid 

subplot (4, 1,4)  ,  plot (tld,beta_sa( : ,1) ,  ' - . ' , tld,beta_sa ( : ,2) ,  ' -- ' , tld,beta_sa ( : 

,3)  ,  '- 

' ) , ylabel ( ' Beta  (deg) ' ) , grid 

xlabel('Time  (sec)')  ,axis([0  5  -5  10]),  pause, 

legend ( 'CMARC  ,' Classical ',' Parameter  Est       ',0) 

%  Plot  Aileron  Doublet  Response 

figure 

subplot (4,1,1)  ,  plot (tld, da_da( : ,1) ), title ( 'LATERAL-DIRECTIONAL  RESPONSE  TO  A 

5  DEGREE  AILERON 

DOUBLET' ) , ylabel ( 'Aileron  (deg)  '), axis (  [0  8  -10  10]), grid 

subplot (4,1,2)  ,  plot (tld,rr_da( :, 1) ,'-.', tld, rr_da ( : ,2) ,  ' --' , tld,rr_da( : , 3) ,  ' - 

' ) , ylabel ( 'p 

(deg/s) ') ,axis( [0  8  -50  60]), grid 

subplot (4,1,3)  ,  plot (tld,phi_da( : , 1) ,  ' - .  ' , tld,phi_da( : , 2)  ,  ' -- ' , tld,phi_da( : , 3) 

, ' - ' ) , ylabel ( 'Bank 

Angle  (deg) '), axis ( [0  8  -30  10]), grid 

subplot (4,1,4)  ,  plot (tld,beta_da( : , 1) ,  ' - . ' , tld,beta_da( : , 2) ,  ' -- ' , tld,beta_da( : 

,3) ,  '- 

' ) , ylabel ( 'Beta  (deg)'), grid 

xlabel('Time  (sec) '), axis ( [0  8  -10  15]),  pause, 

legend ( 'CMARC , 'Classical ', 'Parameter  Est       ',0) 

%  Plot  Rudder  Step  Response 

figure 

subplot (4,  1,  1)  ,  plot (tld, dr_sr( :, 2)  ), title ( 'LATERAL-DIRECTIONAL  RESPONSE  TO  A 

2  DEGREE  RUDDER  STEP 

INPUT' ), ylabel ( 'Rudder  (deg) '), axis ( [0  8  -4  4]), grid 

subplot (4,1,2),  plot (tld,beta_sr ( : , 1) , ' - . ' , tld,beta_sr ( : , 2) , ' -- ' , tld,beta_sr ( : 

,3)  ,  '- 

' ) , ylabel ( 'Beta  (deg) '), axis ( [0  8  -1  3]), grid 

subplot (4,1,3),  plot (tld,rr_sr ( : , 1) , ' - . ' ,tld,rr_sr ( : ,2) , ' -- ' ,tld,rr_sr( : , 3) , ' - 

' ) , ylabel ( 'r 

(deg/s) ') .axis ( [0  8  -10  5]), grid 

subplot (4,1,4)  ,  plot (tld, lr_sr ( :, 1) ,'-.', tld, lr_sr ( : , 2) ,  ' -- ' , tld,lr_sr ( : , 3) ,  ' - 

' ) , ylabel ( 'p 

(deg/s) ' ) ,grid 

xlabel('Time  (sec) '), axis ( [0  8  -5  10]),  pause, 

legend ( 'CMARC , 'Classical ',' Parameter  Est       ',0) 

%  Plot  Rudder  Doublet  Response 

figure 

subplot (4,1,1)  ,  plot (tld,dr_dr( : ,2)  )  , title ( 'LATERAL-DIRECTIONAL  RESPONSE  TO  A 

5  DEGREE  RUDDER 

DOUBLET' ), ylabel ( 'Rudder  (deg) '), axis ( [0  8  -10  10]), grid 

subplot (4,1,2)  ,  plot (tld,beta_dr( : ,1) ,  ' - . ' , tld,beta_dr ( : ,2) ,  ' -- ' , tld,beta_dr ( : 

,3)  ,  '- 

' ) , ylabel ( 'Beta  (deg) '), axis ( [0  8  -10  10]), grid 

subplot (4, 1,3)  ,  plot (tld, rr_dr( : ,1) , '-. ' ,tld,rr_dr(: ,2) ,  •--• , tld, rr_dr ( : , 3) ,  ' - 

' ) , ylabel ( 'r 

(deg/s) ') ,axis( [0  8  -30  30]), grid 

subplot (4,1,4)  ,  plot (tld, lr_dr( : ,1) , '-. ' ,tld,lr_dr( : ,2) ,  •--' ,tld,lr_dr( : ,3) ,  ' - 

' ) , ylabel ( 'p 

(deg) ' ) ,grid 

xlabeK'Time  (sec)'),  ,axis([0  8  -20  20]), pause. 

legend! 'CMARC , 'Classical' ,' Parameter  Est       ',0) 
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